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DUAL CONTROL THEORY. 0 
A, A, Fel'dbaum 


Moscow 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 11, 


pp. 1453-1464 , November, 1960 
Original article submitted March 23, 1960 


Basic formulas are derived and the optimum control algorithm is determined in the general case, first for an open- 
loop and then for a closed-loop nonlinear system of dual control [1]. 


Similarities and differences between the solutions for open and closed systems are indicated, 


1. Derivation of the Risk Formula in 
Open-Loop Systems 








We start first by deriving the formulas for an open 
system shown in Fig, 1, This is done to demonstrate 
certain methods which have led to the obtained results; 
also to compare certain characteristics of the open and 
closed systems, The open system is basically more 
simple than the closed one, and this makes the deriva- 
tion of the formulas easier in the former case, 

The problem is stated as follows [1]: all quantities 
are functions of discrete time at the time moments 
0, 1, see» S» ove, M, Where n is fixed, The input is 


rs = 2% (s, A), (1) 


where A is the parameter vector 
A= (Aj;4.., Aq) (2) 


with an a priori probability density Py(A) = P(A). The 
input x* becomes mixed with noise h* in the channel 
or system H*, at whose output y* is obtained, Statisti- 
cal properties of the noise, and the method of combin- 
ing the signal and the noise in H* are known, Conse- 
quently the conditional probability P(y*,/ x*,) is also 
known, being identical for alls , as the probability 
density P(h*,) of the noise is assumed as not varying 
with s*, 

The characteristic of the controlled object B is 
given by the formula 









































Ip = Fo(Zs, Us), (3) 
where F, is a known function, and 
Zz; => Zs (s, p), (4) 
1h" z 
| ee 
xy LY A u B v C “= 
Fig. 1 


in which p is the parameter vector 
= (Hi, Ha, - ++» Hm) (5) 


with given a priori probability density Po(z) = P(g. 
The quantity v, is obtained at the output of the channel 
G, in which quantity u, is mixed in a known manner 
with the noise g,, whose probability density P(g.) does 
not vary with s, Thus, the conditional probability den- 
sity P(v,| us) is also known, It is now required to find 
such a sequence of probability densities T,(u,|y§-1), 
that the average risk, that is,the mathematical expec- 
tation of the quantity 


s6=—n 


W = > W, (S, 2, 2%), (6) 
s=0 


will be minimum, 
Therefore, the minimum of the quantity 
s=n ‘e=n 
R=MW)=>MW)j= DR ( 


8=0 s=0 


is required, The function I) is the required algorithm 
of the controlling member A, 

Further notation is introduced, Let P(x,| u,) be the 
conditional probability density of x, with given u,. 
This function can be computed from the formulas (3) 
and (4), when the probability density P(g) is known, 
In addition, let P(y$_,| x$-1) be the conditional proba- 
bility density of the vector 


oe “ (y5» yy > meng ¥,_,) (8) 
when the vector 
x, = (x, aie «as x) (9) 


is known, 


It follows from the properties of the channel H* 
that 
i=s—1 


Pry, |jx_)= [] P@lai). ay 
i-o 





* All external noises are regarded as independent, 
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The vector x$-, depends on s and A, and therefore 


P(y$-1| X$-1) depends also on s and A, 

We denote by (Xs, Vg, Ug, Y$-1) the region of va- 
riation of the parameters xg, Vs, Ug, Y$_1- An infinitely 
small element of this region is denoted by 


dQ (x5, Vs, Us, Y;_,) = da, dv, du, dy, ...dy'_, (11) 


At first we shall write down the expression for the 
conditional partial risk r, at the sth moment, under- 
standing by the latter the magnitude of the risk R, when 
the vector xs is kept fixed or, what amounts tothe same, 
when the vector A remains fixed, Then 


re= M{W,|xs} = 
x5 .y Mg ¥,_ 4) 


x P (Us lus) Ts (us | Ey 
P am | x,_,) dQ (Zs, Vs, Us, y,_,): 


W,(s, 23 , %s) P(x%s| vs) X 


(12) 


Let 2(A) be the region of variation of the vector 
X and dQ(A) an infinitely small element of it, The 
partial risk R, is given by the formula 


R, = \ r.P(A)dQ(Q) = 
Pel 


mee 


O(N, Ug. Ug a 


W,(s, x3, &) P(xs| 0s) * 
a) 


. 
s—1 


x P(vs\us) Ts (us| y,_,) Piy, |x 


s—1 
P (A) dQ (x, Us, Ye Y,_y, 4) = 

. \ 
A(X y My Vs ¥y—y) 


{ \ W,[s, 2: (s, 4), 2s] X 


QQ) 


P (a,\ v5) P (vs| us) Ts (Us| ¥3_,) 


x P (yy, |x)? (Ad dQ(A)} dQ (as, ws, vs, ¥%,)- 
(13) 
The expression in the bracesrepresents an integral 
over the region 2(A), as W, and P(y%-_,| x$-1) general- 
ly depend on A . Having performed the integration, we 
obtain a function p, in the braces depending on xg, s, 


¥%-1. The dependence on s is shown in the index 8, and 
therefore one can write 


Ps = Pa (Zs, Y;_,.) = \ W,[s, 23 (8, A), 2s]° 
ai) 


° P(y;_,|x,_,) P (A) dQ (a). (14) 


R= | P (a4 0.) P (vs | ws) T's (tts | y2-1) 


Ax 5,05 stty.¥, 4) 


X Ps(Xs, Y;_,) AQ (Le, Ve, Us, Yz_yle (1) 


From (7) the total risk R is determined from the 
expression 


R=DR=D J — Plalv) Plu) x 
s=0 8-20 A(x ,4,v o%e¥,_) 
xT, (us| Yi) Ps (2s, Yo~s) AQ (as, Vs, Us, Yi) 


(16) 


The functions I, are to be selected in such a way 
that the value of R is a minimum, 


2. Determination of Optimum Strategy 
for Open-Loop Systems 

It can be seen from the formula (16) that the se- 
lection of functions I, for fixed s only affects the com- 
ponent R, corresponding to the sth time moment, In 
this way the total risk is identical with action risk and 
one is allowed to select I, such that it minimizes a 
single R, in (15), As the function I, represents probabi- 
lity density, we have 








| T.(w,)dQ(w,) = 4. (17) 


Q(ug) 


The expression for R(s) is rewritten as follows: 


R= Patualyd{ § Pe|v.) P(ro|us) x 
Aus,,_) QA(Xg Pg) 
X Ps (Xe, ¥5_,) dQ (as, vs) dQ (u,, y;_,)- 
(18) 


The integral in the braces represents a function of 
us and y$_, which we shall denote by & ,(Us, y$-1)» 
that is, 

Es (Us, y;_,) — 


= \ P (%|v,) P (ve | Us) Pe(%s, Y;_,) F2 (Ls, Ve). 
2QA(xXs Vg) (19) 


Then 
Bin \ I (y;_,)dQ(y%_,)s 


aty,_,) 


where 
T(y2_,)= | Te(usfyt_,) 8 (as, yj_,) 42 (ws) = 


Q(u,) 


=(Eday § Ve(tslys_,) €Q (ts) = Edgy > (ednin’ 


2Q(ug) ( 20) 


One can take the average value (&,) , , outside the 
integration in accordance with the integral mean value 








theorem, The function I) represents probability density 
and hence its integral over the region 2(u,) equals 
unity, 

It follows from the expression (20) that a minimum 
value of I(y$-,) is (€s) mine It is not difficult to show 
that I attains this value if the functionT, is chosen ac- 
cordingly, In fact, let u§ be the value of u, correspond- 
ing to the minimum (for instance, the least of all local 
minima) of the function €,(ug) in the region Q(u,), 

We assume that such a value exists either within or on 
the boundary of the region, 

Consider the function 


I’, (us) = d(u.— 4; ), (21) 
where 5 is the unit impulse function (otherwise Dirac's 
5 - function), Formula (21) provides the optimum algo- 
rithm of the member A, as substituting (21) into the 
left-hand side of (20) and using a well-known property 
of 5- function, we fine 


[= \ ly (us | <3 Es (Us, y;_.) dQ (u,) = 


(ug) 


zt i 3 5 (us — ui) Es (us, ys) dQ (u,) = 


= §.(u;)= min €,(u,)=(/) min® (22) 
us EQ(u,) 

The minimum values of I for each y$_; provide 
also the minimum of R,. 

Thus, the optimum strategy, as seen from the for- 
mula (21), does not prove random but regular, The val- 
ue of ug should be taken as equal to a u§. As seen from 
the formula (20), the value u§ which minimizes the 
function & .(us, ¥$-1) depends ony $_;. In other words, 
as 


Es (us, Ys1) = min, (23) 
the quantity u§ is a function of y$_;: 
uy = Us (ys-1). (24) 


This makes also the function r, in the formula (21) 
depend on y$-;: 
e= 5(u,—u;,) = d[u,— uj (ys-,)}. (25) 
Thus, the optimum strategy proves regular; it is 
determined in the general case from all the preceding 
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values of y$ (i=0, ., 81), It should be shown that the 
regularity of the optimum decisions occurs in a wide 
group of problems in Wald's general theory of statistical 
decisions, 

The integral £.(ug, y$_;) can be evaluated with 
the aid of a computer CO (Fig, 2). The block diagram 
of the optimum controlling member can then be rep- 
resented as shown in Fig, 2, The automatic optimizer 
AO selects a value u, such that €, becomes least, This 
value is now sent to the output of the block A, Were it 
possible to determine the value of u§ analytically as 
a function of y$-,, the block A could be constructed 
differently, However, the formula for &.(us, y$-;) is as 
a rule most difficult to obtain explicitly and that is 
why the computer CO must perform automatically the 
integrations with respect to x, and vg, as shown in the 
formula (18), In this way the minimization of the inte- 
gral €, can be achieved directly, 


3, Derivation of Risk Formula for Closed- 
Loop Systems 
Consider now a block diagram of a closed system 
shown in Fig, 3, In order to simplify the exposition we 
neglect the interference h* in the channel H* (see, 
for example, Fig, 1 or Fig. 3 of [1]); x* is thus assumed 
known, It was shown previously how the interference h* 
should be taken into consideration, and should such a 
need arise, this can be done without any real difficul- 
ties, Our notation remains the same, In addition, there 
is also the feedback channel H with interference 
which has a constant in time density distribution P(h,). 
The way in which the signal and noise combine 
in the channel H is given and the conditional probabi- 
lity density ‘P(y | X,) can therefore be found,and con- 
sequently the conditional probability density 


t=s 


P (y. |x.) = [] P (ye|2)- (26) 
i=0 








In a closed system it is necessary to take into ac- 
count also the vectors 


Xs = (20, T1>-++) Ze), 
Yo = (Yor Y1r +++» Ys)» 
U, = (Uo, Uy,~~ +» Us)» (27) 


Ve = (Vo, Vy, ++ +9 Ve)s 
X_ = (Zo, Ty +++) Ts): 


z 
emer r 


h 


Les 











= |* 
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The determination of the optimum controlling 
member A is reducible to the determination of its stra- 
tegy, random in the general case, ensuring the minimi- 
zing of the risk R. The function T,(u, | x", ys-1» Ws-1)s 
is assumed to represent random strategy, that is, the 
conditional probability density of u, when the vectors 
Xs", Ys-1» aNd Us., remain fixed, Now I; is a function 
of 2°, Ys-1» and Us. 

First we shall write down the expression for the 
conditional partial risk rz and we shall now understand 
by it the risk R, when all the preceding motions of the 
system A are known, that is, when the vectors x,*, ys-1, 
and @,., are fixed, Thus, 


r= M {W,|x; » Ys—1, Us_,} <3 


= \ W, (s, x; ° Xs) Pz, | Vs) x 
2 (x goP g,%g) 


X P (v,| Us) Ts (Us| Xs » Yor» Us) dQ (xy, Vs, Us). 
(28) 

In this tormula P,(xs | v,) is understood to repre- 
sent the conditional probability density of the output 
x from the object B when its input y is known; as shown 
below, this probability need not be the same for dif- 
ferent s, 

If one considers the values of r, corresponding to 
different trials of the system, then generally speaking 
the vectors Us_; and ys-;, not known beforehand, can 
for different trials assume different values, Therefore, 
in computing the average partial risk R, it is necessary, 
when averaging f,, to take into account the joint prob- 
ability distributions P(u,_;,. ¥s-1) of the vectors u,-; 
and ys-;, which in the general case are not statistically 
independent of one another. Thus, 


R= M {r,} = 
= \ W.(s, 2% , Xs) Ps (x5| vs) P (vs | Us) X 


O(x as u s Ys— ) 


KT, (us| X51, Us—1, Yo—1) P (us_i, Ys—1) 








dQ (2s, Us, Us, Ys—1)- (29) 
The values of the quantities uj and y; (i=0,..., 
s-1) are memorized in the controlling member A. 
Thus the input and the output of B are known — admit- 
tedly not exactly, but with some error due to the pre- 
sence of g and h, the input and the output interference 
/ respectively, This knowledge enables one to obtain a 
more accurate probability density of the parameter 
vector g on which in turn depends the interference z, 
acting on the object B, The better g is known, the more 


fully deterrnined is the method of regulating the object, 


One can derive a Bayesian formula in order to deter- 
mine the a posteriori probability density of vector g 
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in the sth beat, Consider the probability density t 

P(p, Ws_ 1. Ys-1) Of a joint event consisting in a simul- 
taneous occurrence of the vectors y, U,_; and ys-;. Let 
P(ts_ is Ys-1 |e) be the conditional joint probability 
density of the vectors t,_; and y,_, when p is kept fixed 
and P(p | uy_. Ys-1) is the conditional probability 
density of g when the vectors u,_; and y,_; are fixed, 
The latter is the required a posteriori probability den- 
sity of the vector y. Then, of course, 


P (pt, Us—y, Ysa) = P (Us—1, Yo—s |p) P (p) = 
P (p| Us—y, Yo) P (uy_y, Ys—1)- (30) 


P(g) denotes the a priori probability density of 
vector gp. Hence the formula is obtained giving the a 
posteriori probability density of vector p: 


Po) Puy) Yea |M) _ 
P (u,_,, Ys—,) 





P (p| Usa, Ysa) = 


P (uw) P (u,_,» ¥p_, | #) : 
J P (u,_,. yo_, 1B) P(e) dQ (m) 
2p) 

In the latter formula Q(y) denotes the range of 
variability of vector g and dQ(jp its infinitely small 
element, 

In formula (31) a special role is played by the 
likelihood function P(ts_1. Ys~1! @ ). The probability 
densities of quantities are called likelihood functions 
when the vector g, whose values are subject to hypothe- 
sis testing, is kept fixed, We obtain the likelihood 
function for the closed-loop system shown in Fig. 3, 
Obviously, we have 


(31) 





P (Us1, Ys—1 |) = P (Uo, yolp) P (ui, "1 |", Uo, Yo) X 
x P (Ug, Yo | p, U1, ¥:) eee P (Us—1, Ys—1 |p, Us_», Ys—2): 
(32) 


A typical factor on the right-hand side is a proba- 
bility density 


P (uj, ys|, Wi-1, Yi-1) = 
=P (yi lp, ui, Ui-y, Yi-1) P (ui | Pp, Ui—1, Yi-1) _ 
= P(yi|p, 7, wi) P (us| pp, Wi, Yi) = 


= P(y|p, i, wa) Ts (ui | u—1, Yi—-a)- (33) 


Thus P(y; | p. Uj. Bj-1, ¥j-1) depends neither on 
@j-; nor on y;_;, but only ong, on the beat index i and 
on the value of uj. In addition, P(u; | p, tj-1» Fi-1) 
does not depend on g but only on uj-,; and y;., and this 
is precisely the stochastic algorithm I; of the control- 
ling member A. 


tIn a number of subsequent formulas the dependence 
of the probability density on x? is not stated explicitly 
but it is always implied. 














By substituting the expression (33) into (32) we 
find that 


s—l 


P (ua, Yea| o) = [TJ P@| m4 ad) | 


s—1 


[ I] YP; (uj | Ui-1, Yi—1) ] Py (Up). (34) 
t=] 
Here Po(uo) is the initial probability density of u 
which also should be given. 
Now one is able to substitute the likelihood func- 
tion obtained in (34) into the a posteriori probability 
density formula (31), We obtain 


P, (p) = P (p| usa, Ya) = 


s—1 


8—] 
[[] Pmleiud [FL 
i==0 


P (yp) —— 





P (U,_4» Ys_,) : (39) 


In the above formula the following notation has 
been introduced; 
Do = Po (Uo). (36) 


Knowing the a posteriori probability density P.(p) 
one is able to determine more precisely the expression 
of the probability density P,( x, | v,) which appears in 
the formula (29), From the formulas (3) and (4) we 
find that 

t= Fy (Z., Vs) = F (s, p, Vs) (37) 


The quantities s and v, which appear in the func- 
tion F are parameters, This formula enables one to 
find the conditional probability density P,( xs v,) when 
P(g) is known, 

Since 


P, (Zs | Vs) P (v, | us) dQ (v5) = P, (25 | Us) = 
(Q0g) 
= | P(eslm, ms) Po(p)dQ(p), (38) 
Xp) 
the expression (29) giving R, can be re-written in the 
following way: 


s— \ W, (s, 3, Z,) 
QUxg Us -¥g—1) 
[ \ P, (25| 04) P (04| ws) dO (v4) | x 
QA(r,) 


xP, (us | Ties Us—1, Ys—1) 
P (Ug-1, Ys—1) GQ (Xs, Us, Yo—1) = 


= \ W,,(S, 2, 2») P (xe |p, ws) 
Q(x, Ps Us Yg—y) 


P, (p)Ts (us| Ss: Us—;, Ys—1) X 
x P (us—1, Yo—1) dQ (2, fp, Us, Ye—1)- (39) 


Substituting the formula (35) for P.(g) and can- 
celling P(ts.15 ¥s-1) we obtain 


R, = W.(s, £5 , 2s) P (x,|p, us) P (p) X 
A(Xs Ms .Y¥s—y) 
s—1 5 
x I[ Py; p, i,j) Il TydQ (z,, p, Us, Vs—,). (40) 
i=0 i=0 


The total risk is given by 
s=n s6=n 
aT at a 
s=0 8=0 Q(x, .p,Us .¥g—)) 
W o(8, 2s , Zs) P (5 |p, Us) P (p) x 
s—1 8 
* I] P (yi | p, 7, wi) I] PdQ (25, 1, Us, Ys—s)- 
i=o i=0 (41) 
Let us compare the expressions (16) and (41) of 
risks in open and closed systems, In the formula ( 16) 
only the sth component R, was influenced by the func- 
tion I; the total risk associated with I; was exclusive- 
ly an action risk, In the formula (41), however, the 
function T° influenced not only the term R, (action 
risk) but also all the 9 = components R; where i> s, 


As regards I), the sum = 1 R; represents the investiga- 


tion risk, In fact, the selection of I, influences the 
character of the process in the subsequent beats, as it 
causes either a better or a worse investigation of the 
characteristics of the object B; the latter amounts to 
the determination of the probability density P,(m). 
Thus the total risk regarded as a function of T, (0=s= 
= n) represents the sum of the action risk and the in- 
vestigation risk, As regards Tp, only a single action 
risk occurs . 

When the object B has memory the selection of 
I, influences not only R, but also the terms Rj when 
i>s even when the additional investigation of the ob- 
ject B has not taken place (see Paper IV in this series). 
Consequently, in this case the action risk is expressed 
by a more involved formwvla, and the interconnection 
between the action and the investigation risks proves 
also to be more involved, 
4, Determination of Optimum Strategy 
in Closed Systems 

It is required to select the functions r in the for- 
mula (41) such that the risk R be minimum, As T; rep- 
resents probability densities they must satisfy the con- 
dition (17), 

To be able to select I, we consider the last com- 
ponent R, of R: 

R, = 
AXn# Un. Yn—,) 
Wa(n, Zn Zn) P (en |, Un) P (pe) x 


n—l n—t1 


x ]] Pile. tu) [] Pn dQ Gn, f, Uns Yn—a)-(42) 
i=u i=0 








1037 














We assume that T», ..., ',_; are given and only 
I, is to be selected such that R, becomes minimum 
provided a constraint as in (17) is satisfied: 


Pn (Un, Un—1, Yn—1)dQ (un) = 1. (43) 
A(uy) 
In order to simplify the expression (42), put 


On = On (Un, Yn—1) = 


s= \ W,(n, Tn; Tn) P (Zn |p, Un) 
Q(Xp) 


n—l1 
P(p) []? &!# i, uj) dQ (ap, pr). (44) 
imo 


Further, put ‘ 
8. = IT r;; (45) 
i=o 


R,= 
Quy Un—1-Yn—) 
Xn (Un, Un-1, Yn-1) Bn—s rn dQ (Un, Yn—1) = 


= { Bn—1%n (Un—1, Yn—1) dQ (Un—1, Yn—1); 
QX(Uy 1. ¥n—1) (46) 
where 
%n (Un—1, Yn—1) = 


On (Un, Un—i» Yn—1) rn (Un, Un—1, Yn—1) dQ (Un). 
Q(uy) (47) 


If the T), is selected such that for any set of vec- 
tors Bp.; and yp-; the quantity « , remains least, then 
the integral (46) will also be minimum (bearing in 
mind that the quantity 6,_; is assumed known. By re- 
garding the vectors Bp.1 and Yo-1 in the formulas (44) 
and (47) as parameters, one is able to find a value uy 
such that the function a, is minimum, We put 


Yn = Om (Un, Un—1, Yn—1) = MIN Oy (Un, Un—1, Yn—1)- 
UnEQ(uy) 
(48) 
The quantity u,f, of course, depends on @,-1 and 
Yn-1: 
un = Uy (Un—i, Yn—1)- (49) 


Then the optimum function Ip satisfying the con- 
dition (43) is 


lr, (Un; Un-1, Yn-1) = b (un — Un » (50) 
because in this case 
Hn = On (Un » Un—t» Yn—1) = (On) = (nding? (2) 


This follows from the considerations similar to 
those given above when the conditions (19), (20), and 
(21) were derived. 
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Now the optimum function I,.; will be deterr 
mined, Consider the sum of the last two terms R,_, and 
Rp. First of all, we introduce the function 


Ly = Ay (Uy, UZ-1, Yr) = 


W , (ke, Xe , Zp) P (Ty |p, Ue) x 


Q(x,B) 
k—1 
x P(p) [[ Pilati w)dQ(r%e,p) O<k<n). 


k-1 
There the expression I, P(y; | ws. i, uy) equals 
unity when k=0, 
Now the sum of the last two terms becomes 


Sn—1 — Ras + Ry = 
= AOn—j Bn—1 dQ (Un—1, Yn—2) + 
A(Un—1-Yn—g) 
+ J anBn AQ (ttn, Yn) (53) 
Qn. Yn—)) 


When i=n-2, the Gi are assumed to be known andl, 
is always selected in accordance with the formula (50), 
Then 


: Cs 
Se; = ( Bn—s { \ l n—1An—] dQ (U»—,) > ws 
QAUy—2 ¥n—z)— Mry_)) 
+ \ Pn—1hn (u’, Un—), Yn—)) dQ (Wn—, Yn—1) } . 


Ain —.¥n—y) (54) 


This quantity is least when the expression within 
the braces is least; we write the latter in full as 


%n—1 (Un—2, Yn—2) = 


mis 


QA(uy_1) 


/ De—jtn—1-+- \ 


Avy) 


Pp—sOn (Un » Wn—1, Yn—y) dQ (yn) } dQ (Un) = 


== \ Pa—: (Un—1, Un—2, Yn—2) X 
QX(uy _,) 


x~ { An—| -|- 


On (Un, Un—1, Yn—1) AQ (y,, -1) }aQ (Un—1). (55) 


QX(Vn—4) 


Now consider the function 
Yn—1 = Tn—1 (Un—1, Yn—2) = On-, + 


+ \ On (Un, Un—1, Yn—1) dQ (Yn). 
XAVy—1) ( 56) 
We find the value up*, which is the minimum 
value of this function, This value obviously depends on 
Bn-2 and Yn-2: 


Un—1 = Un—; (Un—2, Yn—2)- (57) 














Then the most suitable function T'p.; such that 
it ensures the minimum of Kp., and satisfies the con- 
dition (17) will be 


Pp—1 (Un—1) Un—2» Yn—2) = 6 (Un—, — Wns). (58) 


By applying analogous arguments one can deter- 
mine the whole sequence of the functions I; using the 
following procedure, 

We introduce the function 


TYn—k = S—*K + 
2X(¥n—x) 
Yn—K+1 (Un—p+1 Un—x » Yn—x) dQ (Yn—k)s (59) 


where yp = %p. Let the minimum of yp., over Up-, 
be denoted by yj}: 


Yn—k = (Yn—k)tn_,p=ur_,. 


= min Yn—k (Un—x; Un—k—1» Yn—k—-1 ); 
Upn—KEQ(Uy_;,) 
Un—k = Un—, (Un—K—13 Yn—K—1)- (60) 


Then the optimum function T*;, is given by the 
formula 
Dan = 8(Un—« — Un_-x), (61) 


where 6 is the unit impulse function, Then 
i= 


(Sn min ( > Rn; ) nisi 


i=0 


= Bn—r—1T3_, dQ (Un—K—1) Yn—x—1), 
O(Up—p—1-Yn—k—1) (62) 


In this way Ij are determined when i=n, n-1, ..., 
eee, 1. AS far as the determination of Ty = Po(ug) is 
concerned, this function may be given beforehand and 
then no selection is required, If the best Ty can be cho- 
sen, the selection takes place in accordance with the 
procedure shown above, When k=n-=1 we arrive at the 
formula 


(St}nin = \ 


12( Ue, Ve) 
The total risk therefore when all Tj (i= 1, s+. 0) 
are optimal, is expressed by the formula 
‘ R= R,+(Si)min = 
- \ W (0, 2 , Xo) P (xo |p, Uo) x 
2( ko, Ue) 


x P (p) Po (uo) dQ (Xo, #, Uo) + 
+ \ Po (Uo) Vi (Vor Yo) dQ (Uo, Yo) = 


Po (Uo) Ti (Yor Yo) G2 (Uo, Yo) 63) 


Q2Q4(Us, Vo) PY 
= | Po(we)| a0(to) + 
2(ue) 
+ | Ti (Hoyo) 42 (yo) ]dQ(u), (64 


2 Yo) 





where 
Ay = Aq (Uo) = 
= | WO, 25 , 0) P (20|p, wo) P (p) dQ (xo, p). 
O(a.) (65) 
We put 


vi (Uo, Yo) dQ (Yo). 
2( ve) (66) 


¥5 = Yo(Uo) = 9 (Uo) + 


It can be seen from the formula (64) that the 
quantity R becomes minimum when 


Po - Po (Uo) =I, = d(u, — U5), (67) 


where the value uf makes the function y (ug) minimum. 

The value uf is dependent on the a priori proba- 
bilities, also generally on the given data, but not on 
observations which are not available at the initial mo- 
ment of time, 

Thus in the case of a closed system as well, the 
optimum strategy proves to be not random but a regu- 
lar one, This conclusion is naturally valid within the 
bounds of the limitations accepted when formulating 
this theory, 

The determination of sequence of functions y p., 
may prove rather tedious in actual examples if elec- 
tronic computers are not used, Even so, it has only 
been possible so far to solve only relatively simple 
problems in which the memory capacity necessary to 
memorize the functions is not too great, 

The obtained system with the optimum control 
algorithm If is, generally speaking, a dual control 
system in which the u; serve not only as directing but 
also as investigating operations, In the absence of in- 
vestigation, the probability density P,( yg) would remain 
equal to Po(g) and the investigation risk would be null. 

It may happen that in a particular case of a dual 
control system the process of investigating the object 
proceeds in the same way, no matter what the values 
of the regulating operations within the admissible 
bandwidth; what is only required is that the regulating 
operations take place, Such a system (an example is 
given in the next paper) is called neutral, In a neutral 
system the investigation risk is independent of the 
quantity u;, but Ps(g) is not identical with Po(p) but 
varies with s, 
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The problem is considered of estimating the accuracy of self-adaptive systems with extremal continuous ad- 


THE DYNAMICS OF SELF-ADAPTIVE SYSTEMS 
WITH EXTREMAL CONTINUOUS ADJUSTMENT 
OF THE COMPENSATING NETWORK IN THE PRESENCE 





justment of the compensating network parameters, where this adjustment is made by using the gradient method, 
and variable controlling signals, variable object parameters, and random disturbances are considered. 


1, Formulation of the Problem 

The development and application of self-adaptive 
systems gives rise to problems in the probability-theo- 
retic analysis of these systems under the influence of 
random disturbances, The solution of these problems 
can be used as a basis for the synthesis and design of 
the systems, The statistical theory of self-adaptive sys- 
tems is in its initial stage of development, At the pre- 
sent time, known results have been obtained in the 
theory of discrete (step) systems of extremal control. 
The mean error of control and the optimal size of the 
steps have been determined for such systems in the case 
of random disturbances of particular types. This work 
has appeared in several papers, The main directions of 
these investigations are well described in [1-4]. 

The statistical dynamics of self-adaptive systems 
with continuous operation have not yet been sufficient- 
ly studied. Moreover, the design of such systems is of 
practical interest because of the simplicity of their 
construction, From the probability-theoretic investiga - 
tions devoted to this type of self-adaptive system, we 
particularly note the works [1, 5, 6], in which the prob- 
lem we are considering has been formulated in connec- 
tion with systems of continuous operation. In [7] the 
dynamics are investigated of systems with continuous 
operation that use the gradient method, Also, in [8], 
the analysis is carried out of a paiticular case of a 
system of extremal continuous control, operating in the 
presence of random disturbances. 

In the present work, we consider the problem of 
analyzing the accuracy of a self-adaptive system with 
continuous extremal adjustment of the compensating 
network parameters, where adjustment is obtained by 
the gradient method in relation to the variation of the 
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external controlling signals and the parameters of the 
controlled object, 

An automatic self-adaptive system with extremal 
adjustment of the compensating network parameters is 
essentially a system of automatic synthesis of a com- 
pensating network with a given structure, Such a system 
must ensure the adjustment of the compensating net- 
work parameters, such that an extremum is maintained 
in some index of the performance of the control pro- 
cess, The block diagram of such a linear system is 
shown in Fig, 1, It consists of a main control network 
and an auxiliary network for the automatic adjustment 
of the compénsating apparatus, The structure of the 
compensating network is given, and the parameters of 
this network depend on the adjustment. 

The main control network consists of the object 
of the control together with the control apparatus and 
the compensating network. The object and the control 
apparatus have m variable parameters y;, ..., y,, and 
are characterized by the linear operator ®p(yj, .... Ym: 
p), where p=d/dt, The compensating device includes 
adjustable parameters x,, ..., x, and can also be 
characterized by a linear operator P(X, voey Xp P)- 
To the control system is fed the useful control signal 
Z(t) of an arbitrary unknown form, and the random dis- 
turbance X(t). The function Z(t) is smooth, and varies 
slowly in comparison with the processes in the main 
control circuit, 

The equations describing the operation of the main 
control network shown in Fig, 1 are 





(1) 


Y = Op (yy, --- + Ym P)V, 
V = ®, (%,.-.-. 5 tm, p) (A+ XI, 
A=Z—Y 
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Fig. 1. Block diagram of a self- adaptive continuous system. 


The auxiliary networks include a functional device 
for determining the control performance index F, a de- 
vice for determining the extremum of the control per- 
formance index, and elements that ensure the adjust- 
ment of the parameters according to the signals 2,, ..., 
Zp. The form of the auxiliary networks for the extre - 
mal adjustment depend on the method used in determi- 
ning the extremum of the control process performance 
index, This control performance index is a functional 
of the parameters of the objects of control and the com- 
pensating device, and also of the characteristics of the 
controlling signals and the disturbances, 

One of the general methods of obtaining the ex- 
tremal adjustment is to determine the gradient of the 
functional in the space of the parameters xX), X, .... Xp 
and to arrange to have the adjustment of these parame- 
ters in the direction of this gradient, The method of 
estimating the control performance index depends on 
the type of problem and on the form of the disturbances 
experienced by the system, 

When the controlling signal Z(t) and the parame- 
ters of the object y,, ..., Yen Vary slowly, and conse- 
quently the processes in the device that adjust the com- 
pensating network parameters are slow relative to the 
processes in the main network, and when the random 
disturbance X(t) is a stationary random function, wide 
application is found in systems of automatic control of 
estimating the control process performance by the mean 
square of the control error, The functional-character- 
ized network for obtaining an estimate of the perfor- 
mance of the control process consists, in this case, of a 
multiplying block for obtaining the square of the con- 
trol error 4= Z—Y and a device for calculating the 
mathematical expectation of the square of the error. 

It is, however, difficult in practice to construct a device 
that will obtain the mean square error A, It is simpler 
to carry out the averaging operation by using an inertial 






link with time constant T, which guarantees accurate 
averaging. 

The equations for such a functional-characterized 
device are 

U= A, (Tp +1) F =U. (2) 

Further auxiliary networks contain synchronous de- 
tectors which include multipliers for calculating the 
product of the signal F and the test signals (i=1, 2, ...,n) 
and the subsequent filters and inertial links for avera- 
ging the product signals Fw, . At the output of the syn- 
chronous detectors, signals are obtained containing quan- 
tities that are proportional to the gradient components 
of the functional F with respect to the parameters x 
(i=1,..., m). The test signals are also fed directly to 
the adjustable parameter network. For continuous ad- 
justment of the compensating network parameters, the 
rates of change of the parameters Xi (i=1, 2, ..., m) are 
proportional to the measured components of the gradi- 
ent, and the equations of the adjusting network are 


6, = Fw, (3) 
(tp + 1) 4 = — xf, 
px = 2 + py (i = 1, 2,. 0-5 a) 


where the quantities T; are the time constants of the 
filters. 

The test signals w; (i= 1, 2, ..., m) for determining 
the components of the gradient of F can be both sinu- 
soidal harmonics and random, unrelated signals, When 
regular test signals are used, their frequency must be 
chosen sufficiently low, so that the signals will pass 
through the filter of the functional- characterized net- 
work forming the quantity F. If random test signals are 
used, they must possess a sufficiently low frequency 
spectrum, 

It is obvious that the time constants satisfy the re- 
lations T;> T(i=1, 2, ..., n). The system of equations 
(1-3) completely describes the operation of a self-a- 
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daptive system, The problem is to estimate the control 
error 4 in the system with various regimes for the adjust- 
ment of the parameters, and to estimate the adjustment 
network parameters kK, and T, that will ensure the neces- 
sary accuracy of control in steady and transient regimes. 


2. The Approximate Equations for Linear 
Self-Adaptive Systems with Extremal 
Continuous Adjustment of the Compen- 
sating Networks 

The direct use of the nonlinear system of equations 
(1-3) for solving the problems posed above is very dif- 
ficult, These problems can, however, be solved by me- 
thods from random function transformation theory. In 
order to apply these methods, we first transform the 
equations, In the neighborhood of the extreme values 
of the variables x;,(i=1, 2, ..., m), the operator cor- 
responding to the compensating network can be written 














, (%1, +++y In, P) _ (4) 
— Dy (Zo, coe Tro» P) + >? uyb, (p) of 0 0 09 
f=} 
where d®, (z;, p) 
b, (p) = [Pe] (v= 4, 2,..., n)s 
v Xj—Xin 


Xjo are the values of the variables x; for which the 
functional F attains its extremum, and u,, = x,y — Xyo 
are the variations of the parameters from their extre- 
mum values, 

We write the variables occurring in the equations 
(1) and (2) in the form 


Y¥=Yo+y ¥,, (5) 


i=1 


V=Voty Vir A= Ac+ 5} A,, 
i=1 


i=1 


n au n 
U =U,+ > oa Se F=Fot > Fi 

i=1 
where the zero subscript indicates values of the varia- 
bles corresponding to the regime of operation for opti- 
mal values of the parameters x; = xig (i= 1, 2, ..., m). 
The components Y; , Vj, 4;, Fj are the variations of 
the variables for deviations of the ith parameter x; from 
its nominal value, If we use the expressions (4) and (5), 
and limit ourselves to small quantities of the first order 
relative to the variations of the parameters and the. 
variables corresponding to these parameters, we obtain 
the following system of equations for determining the 
variations in the variables: 


Y; = Dp (yi, . - - » ym p) Vis 
(Tp + 1) Fy = Uj, 
V; = ®, (Zio, +++» Eng» p) A, -$ C;U;, 


i=1 


(6) 
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n 
0; = Fw; f p> F yw,, 
k=1 


Ay == Y;, (tp +. 1) Z,--_ %:9;, 


0U 
U,= aA; Ai, pu, = % + pw; — Pip 


(f= 1,2,..., 2), 
where cj (t)=b;(p) M[4y + X)(i=1, 2, ..., n), the sym- 
bol M indicating mathematical expectation. 

A system of equations for determining the varia~ 
bles Yy, Vo, 49 can be obtained from the equations (1) 
for xj = xjg (i = 1, 2, ..., m). In the neighborhood of the 
extremum, the functional F, and consequently, the 
functional U, will have zero linear terms relative to 
the ith variable, because of the conditions for an extre- 
mum. Therefore in determining the value of 6U/d4; 
we must take into account the quadratic term, As a re- 
sult, the value of 6U/6A, will be obtained from the 
equality 

au (1) 


n 
oh = Dy Os i= 4,2,..., 0), 
j=1 


In the general case, the investigation of the accu- 
racy of control and the accuracy of adjustment of the 
compensating network parameters is therefore reduced 
to the integration of the system of equations (1) for 
Xj = Xjg and the integration of the system of equations 
(6) with the conditions (7) and the subsequent calcula- 
tion of the probability characteristics of the control 
error A, It should be noted that the analysis of these 
equations does not present any difficulty in principle 
for arbitrarily rapidly varying parameters, or, therefore, 
for any operating regime of the system. 

In the practically important type of case when the 
processes in the parameter adjustment network are slow 
compared with those in the main network, the system 
of equations (6) can be simplified. In this case, after 
the main network has been operating for a certain time, 
the variations uj = x; — xj9 (i= 1, 2, ..., n) can be consi- 
dered as constant parameters. The network equations 
then take the form 


(Tp + 4) (tip + 1) % = — rawt 2 dys — 
j=1 
(8) 

pu; = 2 + pw; — pr, Ek > Deas 

— Gi = Uw; + n 
+>) b,j (uu; + usw, + uw; + ww) w; — wi >} bjt; 

ky i= j=l 

Cyn Ser 
bj = On9j5 (ky f = 1,.2,.-. 4m) 


(8") 


co 
ay =\ ga, (i, t) Cr (t) dt {k= 1,2,..., #). 
0 








The quantity ga, (t, T) is a weight function for the 
output Ay, where this output satisfies the equations 


Y; = Dp (y1, > + +9 Um p) Vi, (9) 
V; = D, (10, . - . , Sng, p) Ay + 5 (t — 1), 
A, ——" Y;, 


where § (t-T) is the delta function. 

In order to obtain the system of equations (8) from 
the equations (6) we use one of the equations from the 
system (1) for X4 = Xjo : 


(Tp + 1) Fo = Us. 


In the case we are considering, the system of equa- 
tions (8) can be used for the analysis of the accuracy of 
adjustment of the parameters x; (i= 1, 2, ..., nm), and 
for the calculation of the mathematical expectation 
My; and the mean square of the deviations of the quan- 
tities uj. These accuracy estimates for the adjustment 
can be applied in the choice of the optimal values of 
the amplification coefficients Ki and the time cons- 
tants 7% and T, and also for the choice of the test signal 
amplitudes w;. 

The final calculation ofthe mathematical expec- 
tancy of the control accuracy for a self-adaptive sys- 
tem is made by using the formula 


(10) 


(11) 
ma ma, + > a\m,,, 
i=<1 


where m, is the mathematical expectation of the con- 
trol pin. in the optimal values of the parameters 
Xj = Xjo, and my, is the mathematical expectation of 
the parameter adjustment errors uj = Xj — Xjo. 
The standard deviation of the control error is cal- 

culated from the formula 

n nn (12) 
Da-~ Da, + 2 > aiK aa, + > > a;0jK 4 .4;, 


ta= } i=1j=1 
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Fig. 2. Self-adaptive system with continuous extre- 


mal adjustment of the compensating network am- 
plification coefficient. 




















where D A, is the standard deviation of the control error 
A, for optimal values of the parameters x; = xj9; Kad: 
K A; A, are the correlation coefficients of the random 
components of the control error, 


3, An Application of the Theory 

We will take as an example to illustrate the appli- 
cation of the theory, a system with one adjustable pa- 
rameter in the compensating network - its amplifica- 
tion k, and with the gradient determined by the syn- 
chronous detector method [11]. Fig. 2 shows a block 
diagram of such a system, The control error A= Z - Y 
is squared and the square averaged in a filter with time 
constant T. A second filter with time constant T is 
used to average the periodic signal of the multiplier of 
the synchronous detector, 

Let the random disturbance X be a stationary ran- 
dom function, let the control input Z = z») = const, and 
let the parameters of the object of control N and Ty 
vary slowly with respect to time (relative to the pro- 
cesses in the main control network). 

On the basis of Fig. 2, the equations of the main 
control network are 


(Tp + 1)¥ = Nk(A+ Xl, A=Z-Y. 





(13) 


If we note that m, =0, then for the mathematical 
expectation and the standard deviation of the control 
error, we have the expressions 
an sageieilien (14) 
ms = TNE’ 

C Nk 2 
0 

For Gx(W) =G»= const, corresponding to a distur- 
bance X in the form of white sound, formula (15) yields 
the expression 





nN? Gok? (16) 
2 (i + Nk) T> 
for the standard deviation, 

If we recall that the performance of the control 
process is measured by the expression F = mi, +0", we 
find that the extremal value k can be obtained from 
the condition oc . 
ar ma + Salean, = 0, (17) 
where m, and 0), are given by the expressions in (14) 
and (16), 

The equation determining the extremal value ky 
in the given case is 


ANGoky + AN*Go3ky + ANG 2k. — 4T 9% = 9. 
(18) 
Equation (18) is best solved graphically, when the 
functions N(t) and T)(t) are given, The result yields 
the optima! value of the amplification factor ky(t) of 
the compensating network. 
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The complete set of equations for the adjustment 
network in the case we are considering is obtained as 
the special case of the equations (2) and (3) for one 
adjustable parameter, 

Because of the slowness of the processes in the 
adjustment network, the approximate equations of the 
type (8) for the adjusting network are 

(19) 
(Tp + 1)(tp + 1)2 = 

= — xw*a*u — “up, pu=z + pw — pho. 

If we neglect small quantities of the third order, 
then on the basis of the expression (8") the disturbing 
function ¢ satisfies the relation 

7 = Uw, (20) 
where Uy = aj is the square of the control error for k=ko. 

Since we have assumed that the variation of the 
variables N, Tp, kg is slow compared with the processes 
in the main network, the coefficient a is given by the 
expression 


0 Pag [teen 


For the analysis of the system of equations (19) and 
(20) that follows, we must first of all transform the ex- 
pression (20) by applying the method of statistical lin- 
earization to the quadrator [10]. Since the quadrator 
has an even-symmetric characteristic, we must apply 
the statistical linearization method in the form [10] 


U,= AS = fy + hy AS, 
where fy= m'a,+ on is the mathematical expectation 


of the function U,; Af the random component of A, having 
a zero mathematical expectation; and k,=2m,_ is the 
statistical amplification coefficient of this random 
component, 

We now use the system of equations (19), apply 
the method of statistical linearization [10] to the quad- 
rator to obtain the mathematical expectation of the 
adjustment error m,,, and obtain the system of linear 


equations 
(Tp + 1) (tp +1) z= 


= — xw*a’m,, — xw (m4, + 04,), (21) 


pm, = 2+ pw — pkg, 

Since the time constants T and T of the filters are 
such that the variables m,, m Ay’? Ay can be assumed 
to vary slowly with the processes occurring in the fil- 
ters, then the system of equations (21) reduces to the 


following equation for m,: 
pm, + xbym, = — pk, + pw — xb, (mi, + 4,]. 
Here (22) 


bi = \ ge (t, t) w® (t)a® (t) drt, 


ow 


t 
be = \ es (t, t} w(t) dt, 
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where g,(t, T) is a weighting function for the left-hand 
side of the first equation of the system (21). 

The quantities b, and b, are mean values of the 
variables a*w’ and w, and are obtained by using two 
filters in the adjustment network, The quantity b, is 
close to zero, while b, * const. The solution of equation 
(22) can be used to obtain the function \ 


m, (i) = 


=m,, (0) e-*bt + [et—a [w (t) — ke (x)] dt. 
0 
In particular, for the test signal w= py sin wt, we 
obtain 
_ fe isiiin al 
_ + — sin @,/. 
[+3 
® 


1 








For finding the random component of the adjust- 
ment error u° (superscript "0" denotes the centralized 
random function), we obtain from the system (19) the 
equation 


(Tp + 1) (tp + 1) 2 = (24) 


= — xw*a*u® — xw2ma,, Ap, pu® = 2. 


The system of equations (24) is linear with vari- 
able coefficients, and can be solved by using linear- 
theory methods. However, if we take into consideration 
what was said about the system of equations (21), the 
system (24) can also be written as 


pu® + xb,u°® = 
=— 2, (tp + 1) (Tp + 1) & = 2xmawds. 


The function A? is a stationary random function 
with a spectral distribution 


(25) 


kg? 
Ty@* + (1 + hol)? 


The function z§ is obtained from the function A? 
by multiplication by the test signal fy sin w,t in the 
synchronous detector and averaging in two stationary 
filters; zs is a stationary random function with a spec- 
tral density of the form [9] 


Ga, (o) _ Go. 





G,,(@) = 


x2 m2 . 
(T +t? (@ + a) + [1 — Ts (@ + oi)? t 
2?mi, Be 


26 
(T + ct) (@ — oi)? + [1 — 7+ oa}: om 





Ga, (@) { 





+. 


The solution of the first equation of the system (25) 
can be used to determine the quantity u’(t) and the 
standard deviation from the expression 








0 


where 


k, (t) n= 


i t 
D,()) = GohN*? mi’ 5 \ e- xbibidty {em ntatnh, (¢ — t2 + th) die, 






0 


£2, (@) cos@t da, 


i 





gs, (@) = 
” Too? + [1 + Nol? 


+ 1 


T+Po+o+i—meror + 





(T + t? (@ — om) + 


In a steady state the standard deviation D,, can 
be obtained by calculating the integral 


ioe] 


<2 Za73..3. 3 ..3 oe ee 28 
Dy == GokoN x mass | 82, () "ee ioe (28) 


Moreover, by analyzing the initial second moment 
of the adjustment error for various values of the para- 
meters K and T, we can choose their optimal values 








to correspond to the minimum of 
,2 
Mimaxt dD, a aa 4 = -{- (29) 
“b> . 
1+ pr 
oF 
Gok2N tz? eC 
ae 6s, (0) a — de 
(1+ kN) § ' a + @? 
" 


The optimal values of the parameters « and T of 
the adjustment network depend on the level of the 
spectral density Gg of the random disturbance, on the 
size of the useful control signal Z», and also on the 
size of the amplification coefficients of the object N 
and the correlation network ko. 

The parameters of the adjusting network can thus 
only be chosen to be optimal for definite calculated 
conditions, 

It can be seen from the expression (29) that the 
adjustment error depends greatly on the test signal 
amplitude, For fixed values of the parameters k andT, 
the error in the adjustment of the amplification coef- 
ficient of the compensating network for varying external 
conditions is obtained for a steady regime from formulas 
(23) and (28). 

We now obtain the mathematical expectation of the 
deviation in the estimation of the accuracy of control 
in the form ma, = amy. The standard deviation of the 
accuracy of control is calculated from the formula 


Da, zs a*D.. 


According to the formulas (11) and (12), the accu- 
tracy of control in the case of a self-adaptive system is 
calculated from the expressions 








ma = mg, + ma, Da = Da, + Da, + 2Kaya,, 








l. 


{1 — 7 (@ — wm PP 





where KA,A, is the correlation moment for the steady 
regime, equal to 
90 


0 


em 


#74) @toaytit 





[ — Tt (@ + )* 


g™ iw 








+ — Tt (@ — w)* + (T + T) i (@ — a) + 1 ||ao, 
and PN aNGoks 
M ~ ZA+FRAN)T, * 
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A method is analyzed for solving the problem of determining the stability of a certain system in which the 
controlled object has a nonlinear characteristic for large initial perturbations. 


1, Statement of the Problem 

The stability of a controlled system with a nonlin- 
ear actuating mechanism was studied in [1, 2, 3]. In 
these papers the perturbed movement of the controlled 
object was described by linear differential equations on 
the basis of linear object characteristics or by lineari- 
zation, As we know, the form of the equations for per- 
turbed motion of a controlled system depends on the 
conditions of the specific problem, We shall study the 
horizontal flight of a controlled aircraft at a constant 
velocity. In this case the equations for perturbed motion 
are of the following form: 


Tp + Up + K(8)+p=0, 
p—B—nB =0, 
H=/(0), c=ap+ Ep+G% —typ. 





(1) 


Here ¥ is the yawing angle; 6 is the slip angle; p 
is the deviation angle of the steering wheel; f(@) is the 
nonlinear characteristic of the servomotor and belongs 
to functions of class A or A’ [2]. The function K(8), 
which characterizes the lateral stability moment, is 
essentially nonlinear (Figs. 1, 2). 

A problem of this type was investigated for the 
general case of sufficiently small perturbations in a mo- 
nograph by A. M, Letov [2]. Letov proved the theorem 
on the stability (in a first approximation) of an auto- 
matic control system with an essentially nonlinear ac- 
tuator, But if we do not limit our analysis to movement 
near the origin, then it is necessary to take into account 
the effect of a change in the sign of the function K(8). 

«In other words, in that case we cannot linearize the 
first equation of the system (1). 


At present a general method for solving this type 
of problem does not yet exist. In this paper we study 
one particular method for analyzing the stability of the 
system (1) and demonstrate that the known Lyapunov 
methods can be extended to a control system with a 
nonlinear object. 


2. Investigation of the Stability in the 
General Case 

In order to investigate the stability of the system 
(1) we make use of the second Lyapunov method. Be- 
fore formulating the Lyapunov function we shall reduce 
the original equations of system (1) to dimensionless 
form, From Fig, 1 it is evident that the nonlinear func- 
tion K(8) cah be replaced approximately by the func- 
tion k8 + mB*, where k and m have opposite signs, Fig. 
1 corresponds to k> 0, m<0, and Fig. 2 corresponds to 
k<0, m> 0, av 

Introducing the new variables 9,=¥%, 22=T—, 








N,=8 and t=TT, we obtain = 

Th = Na, 

Te = Gaia + Goate + GasMs + MaMa” + M25, 

Ns = a + Gass, re 


S = git + GaN + 9sMs + GaMa%s* + Gs%s° — 16 — Sf (0), 


K(f) KA) 
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whose derivative satisfies the relationship 





P ioglt — alT* ees —T(U +/B) 
=e (a — pe Si gy ap 
dm, 7 ons (42171 + Gao + ests) + (5) 
be Qs3= = ae nT oF 
s— WG) 43 , 
am + 1G + Ons (Ne+ 4381s) = — A (n,* + no” + 0%), 
T? 
m= — Woe , where A is a specified positive constant number. 
From Eq.(5) we can find the unknown coefficients 
mm IT? al (UG? — ET) Bj), which are expressed in terms of a1, a9, agg, agg, 
R= e: 2A TTrc’ and A, 


We shall study the function 





‘in wat ke | U + 1B |(UG*— ET) 

ee T? (T? + 1G? . i 
si - eal V=F(m, N, Ns) + 7% 

nkG? k (UG? — ET?) 








ea ae 














= and shall formulate its total time derivative. In view of 
° " % liao (2), we shall have 
) - — mG? _ nmG? m (UG* — ET®) 
| = —p > = —p- + TT —V =—f*y, — 
1 — 1UGs— ers ™ 
T (T* + 1G) ° “i (2211 + Ge2Meo + Goss + Mes” + N25) — 
T? + 1G? 


s=T Bs, 1)=hs +9(0), 9(s)>0. — Fa (M2 + ears) — 


Rename erm " a 


We shall assume that the characteristic equation ra tin oe 
S[qimi + GaM2 + 9sNs + Y2NMNs* +9303* — 75 — sf (9)]. 
—i 1 0 Denoting nj by 7 and making use of relationship 
D(a) =| 91 %2— A a3 | =0 (3) (5), we rewrite —Vin the following form: 
0 1 33 — A 


. 4 4 | er 
—V == An? + = An’? +(5A—m,B sail 
has roots restricted to negative real parts, Under these 2am + 2 4M + (5 2 231) Ns 


conditions the following inequalities must be satisfied 4 - . 7 
on the basis of the Routh-Hurwitz criterion: =(m_Bi2) M3 — (mM2Be2N) Nes + oy A(,? + No” + N3")— 
eg + Aggy <0, ag:433 > 0, (4) — (91 + 2By2) 6M — (G2 + M2Bo2+ 920) ons — 


| — (qs + n2Bos + 930) ons + ¥5" + sho* + s—p(s) <. 


e's ass & 
(Ges + 3) (@an402 — 421 — 43) + Gnndas < 0. This expression can be treated as the sum of two 


) quadratic forms relative to ;, 1», Ns, and 0, whose 
. Then, according to the Lyapunov theorem there exists = coefficients depend on 7. In order to find the sufficient 














1 conditions for stability of the system (2) we formulate 
a positive-definite quadratic form F — — B , y ys 
~ ’ 2. ry — os following two determinants: 
1 i 
y A 0 oy m2Byay 
1 1 
D,=| 0 zA =F MBrsy : 
i i 1 
7 mBiun > mBarn 7 (A— mBsn) 
1 i 
2 A 0 0 "5 (91 + m2Brz) 
i 4 zs 
i} 0 zA4 0 7) (g2 + m2Bes + gen) 
D, = 4 ‘ 
0 0 z > (93 + 2Bes + 9s) 
1 i o 1 2 
Fi t+ mB) FZ (ga + meBee+ a2n) = (91+ m2Bes + gan) 1+ sh 
































y 
é>é, >€>0 z, ed, ‘ E “o 
é; 
z, 
Fi] 
0 z 0 z 
Fig. 3. Fig. 4, 


Expanding them, we obtain 


D, = A? + Byy + Cyn, 
(6) 
D, = A*(y + sh) + Ag + Ban + Cyn’, 


where the coefficients B,, Cy, Ap, Bz, C, can be found 
easily. 

If the inequalities Dy>0, D,>0 had been valid for 
any value of n, then in accordance with the Silvester 
criterion the form for — V would have been definitely 
negative in the entire space (1), >, 13, %); therefore 
the system (2) would have been asymptotically stable 
as a whole, Since the coefficients of the polynomials 
(6) may be negative, the inequalities D,>0, D,>0 are 
satisfied only for a certain finite set of values 7. 

Basing ourselves on this fact, we use the Lyapunov 
theorem to state that if the inequalities D,>0, D,>0 
are valid for values of » which belong to a certain un- 
broken continuum [0, 9] containing the origin, then the 
system (2) will be asymptotically stable for any initial 
perturbations of 7;, 2, and © and for initial perturba- 
tions of 1; such that its absolute values do not exceed 
4%]. The existence of such an > 0 is guaranteed if 
A*(y +sh) + A,> 0, since in view of the continuity con- 
dition it is possible to find a limiting value 7 < 7 * such 
that over the semi-interval [0, n* } the inequalities 
D,> 0, D> 0 are valid, It is evident that n* = min( nf?) 
n(?)), where n() and n(”) are the first roots of the equa- 
tions D,=0, D,=0 on the semi-axis (0, + ©). 


3, Investigating the Stability Conditions 
We now proceed to a detailed study of the appli- 








=z Atm, Bian 
D, — i 
> MBean 
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5A 0 
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D, an 0 > A 





S(Bul+a+9n) $(Balt+a) +7 





cability of the inequalities D,> 0, D, > 0. For simplicity 
we assume that n=0, G*=0, and f(0) belong to the 











Fig. 5, 


class A’ [2}. Then ,=1 3 and the system (2) will be 
written as 


mh = Ne 
Ta = Gait: + Agathe + mani? + n96, (7) 


o = 9M + GaN + GNI — 15 — sf (0), 
where 

a =—(k+al), ay =— Ot) | 
m,=—m, n, = l, 


h = (b+ al), 





E —E 
@=a—4,(U +E), q=—e", 


— El 
i Serta wet, 


s=+, fc) =ho+9(0), h>0, 9(s)>0. 
In that case the determinants D,, D, are written as 


5 (Bul + 91 + 9) 
+ (Baal + 92) 

















where By», Bz, are expressed in terms of A by the follow- 
ing formulas: 


T(i+k+al)A 


oe 
22 Tk + al) (U + IE) * 





Expanding the determinants D, and D, we obtain 





A 2m m*T? (4 +- k + al)? 
D = 41+ k+al ' ~~ (e+ al*(U + 16p sf 


Dy = 24*(F +1) — 
E (k +-al) Em 3 
5 — Fal — 





—4|~fa- 


(1+-k-+<al)ITA E 
—4 [Fane +4 — 7 U +16) 


We introduce the substitutions 








Saline, 4 4 
ng —1=% mn = , 
atz=X, <(1+2)—y+F-=Y 


and select the as yet undefined constant A in the form 
Aa ee are 
“ae 


Making use of these substitutions, we write the suf- 
ficient conditions D,> 0, D,>0 for stability of the sys- 
tem (7) in the following form: 








4 21 2 
1+2>— EG+ts! x0, (8) 
X—Y!— #50, (9% 


where — €[0, £}, E= mn. 

Based on inequality (8), it is possible to plot the 
stability region on the plane of parameters x, y for dif- 
ferent specified limiting values F = m¥ (Figs. 3, 4). An- 
alogously, inequality (9) determines the boundaries of 
the stability region on the X, Y plane for various values 
of F =m? (Fig. 5). 





Inequalities (8), (9) must be satisfied simultaneous- 
ly for different values of F = mi. In view of continuity 
and for an appropriate selection of the parameters a,j, 
E, and h, such a requirement may be satisfied. 


In conclusion the author expresses his appreciation 
to A. M, Letov for a number of valuable suggestions and 
remarks, 
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The paper studies the problems of applying the harmonic linearization method (Krylov and Bogolyubov [1, 2]) 
to an investigation of periodic oscillations in control systems with lag that are described by linear nth order 


differential equations with a lagging argument. 


In [3] a number of specific examples involving the 
application of the asymptotic methods developed by 
N. M. Krylov and N. N. Bogolyubov for investigating 
periodic oscillations in both linear and nonlinear control 
systems with lag are studied without rigorous mathemat- 
ical foundations, The essence of the asymptotic methods 
was reduced to finding the solution for the equation un- 
der study in the form of the sum of its particular solu- 
tions expanded into formal series in the small param- 
eter €, For the case where a nonlinear equation was in- 
vestigated the so-called “generating” solution was deter- 





mined first (i.e., a periodic solution of an ordinary lin- 
ear differential equation with a lagging argument,which 
is derived from the nonlinear equation for €=0), For 
this purpose it is quite natural to begin with an attempt 
at applying asymptotic methods (the method of har- 
monic linearization) in the investigation of periodic 
oscillations in linear systems with lag. Our present 
paper is devoted to precisely this subject. 

1, Assume that a controlled system with lag is 
described by ordinary linear differential nth order equa- 
tions with constant coefficients and a lagging argument; 





n 
Liy OQ) =y™ & + Dd [ay (d) + boy (t — +] = 0, (1) 
v=1 
where L is a linear differential operator; a,,,b,, are con- more rapidly than c|t|% can be written as the sum 
stant real coefficients, and 1,, are real positive constants m 
which characterize the time lag, y (t) = >) onye (0) (3) 
Such an equation was investigated in detail in k=1 


[4-7] and in other papers. Schmidt [5] studied an equa- 
tion of the type (1) and de :nonstrated that if a) the solu- 
tion y(t) is continuous, together with its first n deriva- 
tives, over the entire real axis, and b) for t*> + the 
function y(t) does not increase more rapidly than c|t{* 
(a, c>0 are constants), then the nature of the solution 
for Eq. (1) depends on the existence of real zeros for the 
transcendental characteristic function defined by the 
equation 


n 

L(@) = (jo)" + >} [a, (jo)"—" + 6, (jo) e*%), 

vel (2) 

Schmidt also demonstrated that the function lw) 

has a finite number of real zeros and proved the follow- 
ing theorem: if 2 (w) contains exactly m real zerost 
(m> 0), then a) the homogeneous differential equation 
(1) has m linearly independent particular solutions, and 
b) any solution y(t) which for t> 4 © does not increase 


where a, are arbitrary constants and y, (t) are those 
particular solutions that are contained in any series 
yi(t), ye(t), . . . » ¥m(t) which includes functions of the 
form exp (jAt), texp(jXt), exp(jr"t), . . . ; here X, 
d', and d” respectively denote all simple, all double, 
and all triple (etc.) real zeros of the function (w), 

2. One of the possible Krylov-Bogolyubov methods, 
which is called the harmonic linearization method in 
automatic control theory, can be reduced to the follow- 
ing when applied to nth order equations of the form (1), 
Equation (1) is replaced by an ordinary linear nth order 
differential equation of the following form, which has 
no lag in its argument; 


* Paper presented November 19, 1959 at the Seminar on 
Differential- Difference Equations Held in the Lomonosov 
Moscow State University, 

T We mean «4, Gb, ... » Wm. 




















y™ (t) +5 +. ony" (= 0, (4) 





where c,, (v =1, 2,... , mn) are the coefficients to be 
ieee determined. In order to determine the coefficients c,, 


§ (t) = — Ao sin of, 
9 (.) = — Aw? cos ot, 





Substituting the values of y(t) and its n first deriv- 
atives from Eq. (5) and (6) into (1) and (4), we set the 
coefficients of the sines, cosines, and corresponding 
powers of w equal to one another. As a result we ob- 
tain 


we assume formally that a=-— a sin @T:, 
y () = A cos ot, (5) co = a1 + bi costs — = sin @T:, 
in accordance with the Krylov-Bogolyubov asymptotic 
methods [2]. b,, 
Differentiating Eq. (5) n times with respect to t, we Cn = Any + bn; COS MTn_, — COS On, (1) 
obtain 


Cnt, = An + Dy COS OTy. 


Thus the linear differential nth order equation, with 





| 4 (0) = Aw’ sin ot, (6) periodic solutions of the form (5) that are common to 

| Searle See PR Pye > the differential-difference equation (1), will be 

ia 

| L’ ly (jl= (4 —-; sin ot) y™ (t) + &. (2 + b, cos at, — +" sin ot,4) (8) 


















where L® is a linear differential operator. 
Applying the Euler method [4], consisting of the 


substitution 

y (t) =e, (9) 
to Eq. (1) and (8), we obtain the characteristic equations 
for (1): 


F (A) = a" + > [aA"—" + byAn—Ye**¥) —0 
=i (10) 
and (8): 


P@=(t1—-2 


® 


sin ots) A” + 


+ 2 [(a. + b, cos @T, — “est sin @T,+1) rn] 


+ an + by cos wt, = 0, (11) 


respectively; here F and F* are characteristic functions 
corresponding to the operators L and L*. 

3, We shall prove a number of theorems by means 
of which we can clarify the possibility of applying the 
Krylov-Bogolyubov asymptotic methods to an investiga- 
tien of per.odic oscillations in controlled systems with 
lag. It is assumed that the controlled system can be 
described by nth order linear differential equations with 
constant coefficients and a lagging argument. 

Theorem 1. Assume the characteristic equations for 
the linear differential equations (1) and (8) have purely 
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i x y(™—") (2) + (an + bp COS Wtp) y (t) = 0, 
| 





imaginary roots, Then the characteristic functions 

F (jw) and F*(jw) have an equal number of identical real 
zerost which determine the periodic soiutions of Eq.(1) 
and (8). 

Proof. We know that for the existence of periodic 
solutions of a linear differential equation (even one with 
a lagging argument [5]) it is necessary and sufficient 
that its corresponding characteristic equation have 
purely imaginary roots, 

Assuming \ = jw, we obtain 


F (jo) = (jo)"+ >) (a, (jo)"—” 


v=] 


+ b, (jo)"— e~*¥*) = 0, (12) 


F* (jo) = (4 — = sin ot) (jo)" + 


n—1 


+ > [(a. + b, cos ot, — met sin ©, 41) (joy—| 


v=] 
+ an + 6, cos ot, = 0. (13) 


Using the Euler formulas exp(-jwt,)=coswt, — 
—jsinwtp, it is easy to see that Expression (13) can be 


transformed to if 
F’ (jo) = 


=(jo)" + Dy (av + bye) (ja)""}] =0. 4) 


v=1 





tWe mean the roots uy, u,... 


From this it follows that for w equal to the zeros of 
the characteristic functions we have F (jw) = FYjw)=w) 

Based on the results cited by Schmidt [5], it canbe 
stated that the characteristic functions F(jw) and FXjw) 
have the same number of identically equal real zeros 
which determine the periodic solutions of Eq. (1) and (8). 
Thus we have proved the theorem. 

Theorem 2, Assume the periodic solutions of Eq{1) 
are determined by the initial functions 9(*)(t) (r= 
=0,1,... , 2m~1) specified over the initial set Et, 
consisting of the points on the time segment tp-T<t<tp; 
assume further that for t=tg=0 these functions have the 
value 9) (t,), Now assume that the characteristic func- 
tions F(jw) and F*(jw) have m real zeros. Then the dif- 
ferential equation (8) has periodic solutions of the form 
y(t) =B8 ay sin w, t + By, cosw, t(k=1, 2,..., m) 
{where 5,-,, 5, are arbitrary constants, uw, are the real 
roots of the characteristic functions F(jw) (or F{jw))) 
which are common to Eq. (1), provided that the arbi- 
trary constants 8,),, 85, of the periodic solutions for 
(8) are identically equal to the arbitrary constants of 
the periodic solutions for (1). 

Proof, In order to determine the arbitrary constants 
Bs, 89, of the periodic solutions for (1) (the quasipol- 
ynomial of these constants has m real zeros) it proves 
sufficient to specify the numerical values of the initial 
function 9(t) and its (2m-1)th derivative at the point 
t= ty. 

From the Schmidt Theorem it follows that any solu- 
tion y(t) of the differential equation (1) (which satisfies 
the conditions assumed in the theorem) can be written 
in the form of the sim 


y (t) = >) ye (2). 


k=1 


(15) 


In accordance with Theorem 2, y;(t)™ yj; (t) (k = 
=1,2,...,m) where yy, (0 are those particular solu- 
tions which are determined by the real zeros of the char- 
acteristic function F{jw). Here the arbitrary constants 
Bok, 8 can easily be determined from the solution of 
algebraic equations which is obtained if Eq. (15) is dif- 
ferentiated (2m-1) times with respect to t,and the result- 
ing expressions are set equal to the numerical values of 
the initial function 9(tg) and its (2m-1) derivatives for 
t=t,=0. Thus Eq. (8) has periodic solutions 


y (t) = Box—1 Sin @gt + Ba, COS@st (kh = 1,..., m) 


(16) 


that are common to Eq. (1). 

Thus the theorem has been proved. 

Theorem 3, The sum of the particular solutions for 
Eq. (1), 


y(t) = > (Box—; COS @yt +B, Sin @,t), (17) 


k=1 









(each of these solutions is found using the Krylov- 
Bogolyubov asymptotic method) satisfies the differential- 
difference equation (1). 

Proof, The statement of the theorem stems from 
the fact that the particular solutions of (1) that are found 
using the asymptotic methods coincide exactly with the 
particular periodic solutions of that equation. The sum 
of the latter is also a solution of the differential-differ- 
ence equation on the basis of the Schmidt Theorem 
cited above. The theorem has been proved. 

Ir 1959 the paper by A. Khalanai [11}was published; 
this paper provided the foundations of the averaging 
method as applied to systems of quasilinear differential 
equations with a lagging argument, It proved the theo- 
rem (analogous to the corresponding theorem formulated 
by N. M. Bogolyubov [12] for systems of differential 
equations without lag) stating that the error in deter- 
mining the periodic solution of the averaged system is 
small compared to error in the periodic solution of the 
original system of differential equations with a lagging 
argument, It should be noted that the theorem formu- 
lated by A. Khalanai was proved for a system of differ- 
ential equations with a lagging argument reduced to 
standard form [2}. However, differential equations with 
a lagging argument can be reduced to standard form by 
substitution of variables only for the case where the 
linear terms in the system do not contain lag. Thus in 
the general case the problem of reducing a system of 
differential equations with a lagging argument to stand- 
ard form should be considered still unsolved, Such re- 
duction is obviously complicated (even when we seek 
only periodic oscillations) by the fact that the second 
order differential equation with a lagging argument 
which describes the periodic oscillatory mode may 
have two proper frequencies (this is in contrast to an 
ordinary differential equation and was demonstrated 
above). This result also follows from [ 13}. 

4. We shall study the equality between the periodic 
solutions for nth order equations of the form (1) with a 
lagging argument and periodic solutions of ordinary dif- 
ferential equations of the 2m-th order. We shall assume 
that an ordinary linear equation of the 2m-th order has 
periodic solutions equal to those of an nth order differ- 
ential equation with a lagging argument if; a) the char- 
acteristic equations corresponding to the differential 
equations with and without lag in the argument have 
equal purely imaginary roots; b) the periodic solutions 
are determined by equal arbitrary constants; c) the origi- 
nal equations are satisfied for any combination of peri- 
odic solutions of the form 


y (t) = B, sin @,t + B, cos wf + B, sin gt + 
+ By COS Wet + - - - + Bam—1 SIN Omit + Bom COS Wml. 





We shall prove a number of theorems which can be 
used to provide the foundation for postulating equality 
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between the periodic solutions of differential equations 
with a lagging argument (1) and the periodic solutions 
of equations without lag. 

Theorem 4. If the quasipolynomial F (jw) of the 
linear differential nth order equation (1) has m real roots, 
then; a) the characteristic equation for the ordinary 
linear differential equation which has purely imaginary 
roots equal to those of Eq. (1) is of degree 2m; b) itscor- 
responding differential equation is of order 2m. 

Proof, ‘From the general theory of linear differen- 
tial equations it is known that the maximum number of 
linearly independent solutions of a linear differential 
equation is equal to its order. It is also known that if 
the characteristic equation has 2m purely imaginary 
roots, then the roots can only be conjugate pairs. There- 
fore the characteristic equation for a linear differential 
equation without lag that has purely imaginary roots 
equal to the characteristic equation for the linear differ- 
ential equation (1) will have the degree 2m; the differ- 
ential equation that corresponds to it will have the order 
2m. The theorem has been proved. 

From Theorem 1 it follows directly that the char- 
acteristic quasipolynomial F(jw) has a finite number of 
real roots, Assume that the number of purely imaginary 
zeros of the quasipolynomial F(jw) is equal to m. Then 
from Theorem 4 it follows that the characteristic equa- 
tion for a linear differential equation with no lag and 
periodic solutions identically equal to the solutions of 
the equation with a lagging argument is of degree 2m. 
We are, confronted with the problem of how to determine 
the coefficients of this characteristic equation, We 
know that between the coefficients of the algebraic 
equation 


VW 4+ art 2. + dam = 0 (18) 


and its roots A;, A», . . 
relationship: 


- » Agm there is the following 


Nyt eg He. hem = Dy hy = — Oe 


i=1 
2m 
Mare + MAs +. . «+ Aem—sAom = 2 AA; = ae, 
i, j=1 
(i<j) 


Naga + Agheha +. - et Aam—2hom—14em = 
“ne > MyAjhe = — ay, 


i, i, k=1 
G<i< k) 


Aye eee Vem = (— 1)?"dom- 


(19) 


Thus if we know the purely imaginary roots of the 
characteristic equation for the linear differential equa- 
tion (1), 
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A, = + joy, Ag = — jo,... 





Kom—1 = + JOm, dam = — J®m, 


then it is always possible to determine the coefficients 
of the characteristic equation in unique fashion; thus it 
is possible to determine the coefficients of an ordinary 
linear 2m-th order differential equation in the form 


y™™) (2) + ary—D (1) + 
+ ary?" —2) (t) +. ..+ dgmy (t) = 0. (20) 


Here it is easy to show that the coefficients ay, ag, 
as, ... for the odd-order derivatives will always be 
equal to zero, 

Theorem 5. Assume the periodic solutions of Eq.() 
are determined by the initial functions ¢(")(t) (r = 
=0,1,... , 2m-1) specified over the initial set Et, 
assume further that for t=t,g=0 these functions have the 
values g(t) (tg), and the periodic solutions of Eq. (20) are 
determined by the initial conditions for t= t,=0, yf”) 
(r=0,1,..., 2m-1). Now assume the characteristic 
equations corresponding to the differential equations (1) 
and (20) have 2m equal purely imaginary roots, Then 
for ¢) (t4)= yf") (r=0, i, ... . 2m-1) the differential 
equations (1) and (20) have equal particular periodic 
solutions y, (t), and any solution y(t) can be written in 
the form of the sum 


™m 
y() = 2 yw (d, 
k=1 
where y,(t) are those particular solutions which corre- 
spond to the purely imaginary roots of the characteristic 
equations, 

Proof, The equality between the particular period- 
ic solutions y, (0 for Eq. (1) and (20) follows directly 
from the equality between the purely imaginary roots of 
the characteristic equations and the equality between 
the arbitrary constants determined from identical initial 
conditions. The second postulate in the theorem for 
Eq. (1) stems from the Schmidt Theorem; for Eq. (20) it 
stems from the determination of its coefficients, The 
theorem is proved. 

5. As an example that illustrates the theorems 
proved above, we shall study the equation* * 


y(t) — 2y (i _ 7) + 3y () = 0, (21) 


whose periodic solutions are determined by the initial 
functions g) (t) (r=0, 1, 2, 3) that have the values 
9€) (14) for t= tg= 0. 





* * Equation (21) was presented to the author for solution 
by A, M. Zverkin in the course of discussing this paper 
at the Seminar on Differential- Difference Equations in 
the Lomonosov Moscow State University. 













The quasipolynomial corresponding to Eq. (21), 


wi 
E(\)=M—2e* 4+3=0; 


has purely imaginary roots \ = +j, A =43j (this can 
easily be determined by direct substitution); the periodic 


solution of Eq, (21) corresponding to these roots is written 
as 


y (t) = Bi sint + Becost + 


+ Bs sin 3¢ + Ba cos 3¢. 
We shall make use of the Krylov-Bogolyubov 
asymptotic method for determining the periodic solutions 


of Eq. (21). For this purpose we transform Eq. (21) to the 
form 


2sin 5 @ ‘ 
(14 cet Hs y (t) — 2 cos 5 oy (t) + 
+ 3y (t) =0. = 
The corresponding characteristic equation will be 
2sin > 
F° ay=( i+ mrged Ww— 2 cos + wh + 

(23) 

fb 3 = 0. 


In order for the differential equation (22) to have 
periodic solutions it is necessary and sufficient that 
A=jw. Substituting \ = jw into (23) and separating real 
and imaginary parts, we obtain 


s 2 sin 4 @ 
2a cos o = 0, —w* i + = 0. 
(24) 
From the simultaneous solution of Eq. (24) we findw =1; 
w=3, Thus \ =+j, \=+3j, and the example does not 
contradict Theorem 1. The periodic solutions of Eq.22 
corresponding to the particular values of w are 


y (t) = Bi sin ¢ + 


(25) 


+ Pecost, y (t) = Bs sin 3¢ + Ba cos 3¢. 


The arbitrary constants in expression (25) are deter- 
mined from the solution of the algebraic equations 


P (toc) = Ba + Ba, (te) = Br + 38s, 
@ (to) = —B2— 98a, @ (to) = —B1 — 278s (26) 


and will be equal to 


g, = P(t) + 9 (te) 
1 3 ’ 








__ 9 (te) + (to) 
8, = = Ly ‘ 


8 = — 





@ (to) + @ (to) By = — ited + @ (te) 
24 , ba Esa 

(27) 

respectively. 





Thus Eq. (21) and (22) have general periodic solu- 
tions of the form (25). The results which we have ob- 
tained do not contradict Theorem 2. It is easy to see 
that any solution y(t) that equals the sum of particular 
periodic solutions y(t)=6,sint + 8,cost + B,sin3t + 

+ 8, cos 3t satisfies Eq. (21). The latter does not contra- 
dict Theorern 3. 

Further, the ordinary linear differential equation 
whose characteristic equation has the purely imaginary 
roots A; 2= +j, Ag 4= +3j and whose coefficients are de- 
termined from formulas (19) is written as follows; 


“Y(t) + 109 (t) + 9y (t) = 0 (28) 


(the form above can easily be verified). The order of 
Eq. (28) is 2m=4, The latter result does not contradict 
Theorem 4, Assume the periodic solutions of Eq. (28) 
are determined by the initial conditions for t=tg=0, 
yl" g(t) (ts) (r= 0, 1, 2, 3). Then the general solution 
of Eq. (28) will be 


y (t) = B, sint + B, cost + B, sin 3¢ + f, cos 3¢, 


in which the arbitrary constants 6,, 8,, 85, 8, are deter- 
mined from the solution of the algebraic equation (26). 
This result does not contradict Theorem 5. 

The author considers it his duty to express his sin- 
cere appreciation to all participants of the Seminar on 
Differential- Difference Equations in the Lomonosov 
Moscow State University conducted under the supervision 
of L. £. E1'sgol'ts, who took an active part in discussing 
this paper and helped substantially in the formulation of 
the paper. 
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The paper defines certain equivalent transformations of sequential machines and the corresponding structures of 


the machines, 


1, The Structure of a Sequential Machine 

We specify the sequential machine A (in abbrevi- 
ated notation s.m.A) which is defined by the system of 
functions 





s=fi(g), 9g = fe (x, 4q) 


0<sgm, 1<2rgqm,, 1 [qa <Kmey). (*) 


The values x are called input values; s are called 
output values; q are called numbers of states (or in ab- 
breviated form x, s, q states), The time is varied in 
discrete fashion, and the values x, s, q vary only at the 
discrete instants of time ty, ts,... The “prime” de- 
notes the state of the s,m, at the next instant, 

By specifying the series of input values x;, Xg,.. . 

. +X, at the instants t,, tg,... tp and the initial state 
do, we shall obtain the series of output values s;, S2,... 
- 8). The two series xy, Xg,... » Xp, and Ss, S,... 
- » 8p are called an experiment of length n. Two ex- 
periments are equal if both their input and output series 
are equal, 

For a fixed input value x the equation q*=f,(x, q) 
defines a certain transformation of states q (i.e., it de- 
fines a certain generalized substitution ax of the ele- 
ments q), We shall denote the operator for this substitu- 
tion by the symbol Ax, and shall write the equation q’ = 
= f(x, q) as q’=Ayq. Such a notation is sometimes 
more convenient, We depict the substitution a, graphi- 
cally as follows; a) we mark mg points that are num- 
bered 1,2,..., Mg; b) the transition from the preced- 
ing state q to the next state q’ is indicated by an arrow 
directed from q to q’ (Fig. 1). 

The graphical representation of the substitution, 
which indicates the corresponding value of s near q, is 
called the structure of the substitution, The structures 
of all substitutions a, form the structure [1] of s.m.A 
(Fig. 2). The structure of an s.m, defines it completely. 

Discrete (digital) automatic systems are specified 
in various ways. But for synthesis of such systems we 


specify them in the form of s.m, (cf., for exemple, [3]). 
Depending on the method used for this transition, vari- 
ous s,m, are obtained, 

Thus the same digital automatic system can be 
described by different pairs of functions (*). Thus we 
are confronted with the problem of choosing the opti- 
mum (in a certain sense) s,m, for synthesizing the speci- 
fied digital automatic system, 


2. Equivalent Transformations of Sequen- 
tial Machines 

Two s,m. A and B are called equivalent if for any 
experiment in machine A it is possible to indicate an 
equal experiment in machine B and vice versa, The 
transition from s.m. A to an equivalent s,m, B is called 
an equivalent transformation of s.m. A to s,m. B; we 
shall write this transformation as follows; B= aA where 
a denotes the transformation operator, The simplest 
operator a is the operator for the enumeration of states 
when each state q is assigned a different number in such 
a way that different numbers q correspond to different 
new numbers, We shall assume that such an enumera- 
tion does not alter the structure of the s.m. However, 
note that a different numbering of states q leads to a 
different realization; therefore a successful choice of 
the numbering for q can lead to a structurally simpler 
automatic system. 

Since all equivalent s,m. perform an identical 
manipulation of the input data we encounter the prob- 
lem of choosing an s.m, which will satisfy the specified 
conditions in optimum fashion (in particular, such a 
condition as the smallest number of different states), In 
order to achieve such a choice there must be a4 reserve 
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of definite basic equivalent transformations, Below we 
shall indicate certain elementary transformations, but 

they are insufficient for transferring each s,m. into any 
equivalent s.m. The problem involving the complete- 
ness of the set of basic equivalent transformations will 

be studied in a different paper. 


3. Terminology and Notation 








The set of states qj having the property that for a 
given q it is possible to indicate a natural number r; 
for which Al q, = is defined as the center of the sub- 
stitution a,, The sub-sets of the center which are in- 
variant relative to A, are called kernels of the substitu- 
tion, If the equation A, q=q, does not have a solution, 
then the element q, is called the extreme element of 
the substitution. The chain qo, 4;=Ay49,.--» dy = 
= Ay dq is called a trunk of the substitution a, with a 
beginning q, (the extreme element) and an end q,, ; the 
set of trunks with one end q,, is called a tree with the 
base q,,. 

Figure 1 shows the graphical representation of the 
substitution. The elements 1, 2, 3, and 9, 10 formtwo 
kernels, The two kernels form a center, The elements 
6, 7, 8 are extreme elements, 

The two trunks 7,5 and 6, 5 form a tree, In turn, 
the tree forms the trunks 7, 5, 4, and 6, 5, 4, 

We shall denote an experiment by the letter 6. If 
we desire to indicate that the experiment ¢ leads tothe 
state q, begins with the state q, or contains q in the 
middle, we write 64, q6, 6,46, respectively. If the 
experiment q; 6q, exists, then q, is called the state 
which precedes qg, and q, is called the state which fol- 
lows q;. 





4. Identification of Equivalent States 

If the experiments q,6 and q,é are equal, then q, 
G2 are called equivalent states relative to the experi- 
ment 6, If q;, dz are equivalent for any &, then q, and 
dz are called equivalent states and are denoted by qy~qg, 
For the experimenter the states q, and dy are indistin- 
guishable, It is evident that 





1° if §=4i ~ Qjsthen 9; ~ 4; 
2° if qj~ qj and gj ~ q, then q ~ qv; 
3° if i ~~ Qiy then for any = 

we have A.g; ~ Axqj- 


Therefore, the set of states can be subdivided into 
classes Q;, Q9,..., of equivalent states. We shall 
not study classes consisting of one element, The struc- 
ture of the s,m, determines those equivalent transforma- 
tions which can be performed on them, Therefore, we 
shall study the structure of a s.m. when there are equiv- 
alent states, Since the structure of a s,m. is formed by 
the structures of a substitution, it is sufficient to study 
the structure of an individual substitution. 

Assume PR, Pp, ... , Py are classes of equivalent 
states located on a given trunk with a kernel, We shall 
introduce the set M=P,+P,+...+P,. Let qy € M denote 
the element preceding all the remaining elements of M 
on the specified trunk. If M is contained only in the 
kernel, then q, can be any element in M. 

Assume q,, = Aya is the first equivalent state fol- 
lowing q;. Then from the property 3° it follows that 
Ay dg™ Ay, etc. There the function s=f,(q) will be a 
periodic function, beginning with q,, on the given trunk, 
Here the number of elements in the kernel is divisible 











by the period. If the equivalent states are located on 
different trunks of the substitution, then f,(q) acquire 
equal values over segments of the equivalent states, By 
identifying equivalent states of s.m.A we obtain (cf. 
Moore [2]) the s.m.B which is equivalent to the original 
one. Such a transition consists of the following: in the 
substitution ay we retain one of the periods f,(q) which 
forms a new kernel; the equivalent portions located on 
different trunks are joined. 

Example (Fig. 2). In machine A we have two 
classes of equivalent states; Qy { 3, 5, 7, 9} » Qe { 6, 8, 
10}. The transition from A to B is achieved as follows. 
In the substitution a, we join the equivalent states 3, 5, 
9 on different trunks with a common base; we leave the 
states which are not included in the classes Q, and Qs, 
and one period later we replace the designations of the 
elements included in the period with the designation of 
the corresponding class. In Fig. 2, b, c, d we have shown 
the stages of the transition, The transition to equivalent 
machines by identifying equivalent states reduces the 
number of states. A machine without equivalent states 
is called a simple machine. In the analysis below we 
shall study only simple machines, 

The indicated structure of an s.m. when equivalent 
states are present can serve as the necessary identifying 
feature for detecting equivalent states, But note that if 
Q;, Qa, . . « » OQ, are classes of equivalent states relative 
to each substitution, then this does not yet mean that 
they will be classes of equivalent states, For example 
(Fig. 3) the states 1, 3 are equivalent with respect to 
each substitution, but 1 is not equivalent to 3. In order 
to verify this it is sufficient to study two experiments: 





















































5. Annulment of Regular Sub-Sets 
Assume the state q of s.m.A has the following prop- 
erty: for any experiment with the initial state q there is 
an equal experiment in A. The question arises as to 
whether this state q can be discarded. But if we discard 
the state q then the structure of the machine is disrupted 
since the states preceding q will not have a next state. 
Therefore the annulment of one state q involves the 
annulment of all preceding states, We shall formulate 
this more precisely, The sub-set R of the set of states 
in s.m.A is called regular if together with q€R all the 








preceding q also belong to R, and if for an experiment 
with its initial state in R there is an equal experiment 
with its initial state not in R. 

Theorem 1, The annulment of regular sub-sets 
transforms an s.m. into an equivalent one (cf, Appendix 

Example. We shall study machines A (Fig. 4a) and 
B (Fig, 4b). 

It is evident that A~B. In order to verify this it is 
sufficient to check state 5, For experiments 6, begin- 
ning with x=1 the state 5~1; for experiments 6, begin- 
ning with x=2 the state 5~4, The transition from A to 
B is achieved by discarding the state 5 rather than by 
identifying equivalent states, 


6. Partial Enumeration 

As we have already noted, enumeration is an equiv- 
alent transformation. We shall now determine the struc- 
ture required for s.m.A in order to ensure that the enu- 
meration of states over less than all substitutions will 
be an equivalent transformation. 

The states qy, dg,... . 4, which are such that from 
the presence of just experiment 6q;,it follows that there 
exist r equal experiments 6q, (v =1, 2,... ,1),are 
called similar states, 

Theorem 2. If the states qy, dg,... , dy are similar, 
then any permutation q), qg,... » 4, in any substitution 
of s.m.A transforms A to the equivalent s.m.B for which 
Gi, Wg,» - - » Gp remain similar (for the proof, cf. Appen- 
dix). 

_ Example, We shall study machine A (Fig. 5a), 
Machine A does not contain equivalent states and regular 
parts, It would seem that it is impossible to reduce the 
number of states. Note, however, that the states 1, 3 
are similar. Interchanging the positions of states 1, 3 in 
the substitution ag, we make the transition to the equiv- 
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alent machine B (Fig. 5b). In machine B the state 3 can 
be discarded, We obtain D~ A with a smaller number 
of states (Fig. 5c). As we can see, the simplification was 
achieved at the expense of the possibility of enumerat- 
ing the states in one substitution. 

The machine is called “simplest” if for each initial 
state there is an experiment which distinguishes this state 
from others. If there ts just one substitution in the 
machine for which there are no equivalent states rela- 
tive to the substitution, then the machine is simplest. 
For the simplest s.m.A the sufficient condition formu- 
lated in Theorem 2 for the equivalence of partial enu- 
meration is also necessary; this flows from Theorem 3 
(for proof, cf. Appendix). 

Theorem 3. We specify the simplest machines Il, 
and IIx, Assume q;>q,~. . . > 4; ~ q; is one of the 
cycles for the operator 8, If the machines A=II, + Il, 
and B= 611, + Il, are equivalent, then the elements 
ys Gg,» » + » Gy are similar, 


7. The Search for Equivalent Transforma- 








tions 

The structural properties of s.m. studied above 
make it possible to find equivalent transformations in- 
volving identification and annulment of regular subsets, 
and partial enumeration. The corollaries of Theorem 4 
and its proof cited in the Appendix can facilitate the 
search for equivalent transformations considerably. 

Theorem 4, If in the simple s.m.A there is just one 
state q, such that for any experiment qo6 there is an 
equal experiment q}& in the simple s.m.B (where qj can 
be different for different experiments q,6), then there 
are identical submachines A® and B® in s.m.A and 
s.m.B. 

The theorem is proven in two stages, First we find 
the states q, €A and qi €B which correspond to one 
another for any experiment 6, Then it is demonstrated 
that the states in A which follow q, and the states in B 
which follow q} form identical structures for the two 
submachines A® and B® (for proof, cf. Appendix). Note 
merely that if q, and qi correspond to one another for 


a Q—© 


z=! 


any 6, then the state A, qy corresponds to the state 
B, 4} for any &. Therefore, such corresponding states 
q, and qi must be sought for in the kernels of the sub- 
stitutions, 

Assume that in s.m.A there is a substitution whose 
center consists of one kernel, In that case any state in 
this kernel can be chosen as q;. All states following 
q; form the submachine A*. Therefore an investigation 
of the possibility of applying equivalent transformations 
to s.m.A must be performed among states which are not 
included in submachine A*. 


Assume now that s.m.A has no substitutions whose 
centers consist of one kernel. Without undertaking a 
more detailed study of the problem, we shall say that 
it is not possible to accurately indicate which kernel 
of the substitution a, contains the state q,. In that 
case we formulate the submachines Aj (v =1,2,..., 
r; rk) (certain of these submachines may _ coincide 
or may be contained in one another) for each of the k 
kernels of some one substitution a,. 


It is necessary to perform investigations of the pos- 
sibility of applying equivalent transformations to s.m.A 
among those states of A which are not included in sub- 
machine Af; then we investigate the states of A which 
are not included in Aj, etc. 


Let us restate the fact that these conclusions apply 
to simple s.m. (i.e., to machines without equivalent 
states), 


APPENDIX 


Proof of Theorem 1. If the experiment ¢ begins 
with q eR, then it does not contain states in R in view 
of the regularity of R. Therefore, in machine B the 
same experiment & corresponds to this experiment. If 
the experiment & begins with the state q, €R, then ac- 
cording to the condition of the theorem there exists 
de ER, and qg~4q, relative to 6, Therefore, in machine 
B the specified experiment with the initial state q, 
must correspond to an equal experiment with the initial 
state qo. 














Proof of Theorem 2. A portion of the substitutions of 
an s.m. forms the structure of a certain machine I,, and 
the remaining portion of the substitutions forms the struc- 
ture of another submachine II,. We shall represent the 
original s,m. with the sum of these submachines, 

Assume A=TI, + Il,, B=8T, + Mz, where 6 is the 
transposition of similar elements q,, qa. We shall not 
only prove that the theorem is satisfied, but we shall also 
demonstrate that for equal corresponding experiments &, 
and 6, the last states are identical. For experiments 6, 
of length 1 the theorem is obvious, since these 6, can 
serve as the experiment 6, Assume the theorem and the 
corollary are valid for experiments of length ?<m, We 
shall prove that for this case the theorem is valid for ex- 
periments of length m, The experiment & is written in 
the form L, Gacy « 

We shall study two cases, 

A, Assume q, coincides with q, or dg, and q,, sepa- 
rates II, from Il,; i,e., the next-to-last input value xq 
belongs to II, and the last value xy belongs to II,. In 
order to be specific we shall assume qq =q;. We shall 
formulate the experiment 6;, which is equal to L493 4y » 
in the following manner, In machine A we take the ex- 
periment 6/'q, equal to 6,q;. For machine B we take the 
experiment 6,4, (which is of length m~1) to be equal to 
6/2. Continuing it by means of the input value x, , we 
obtain the experiment 6"q,qy which is equal to $141 4y : 

B. If the conditions for case "A* are not satisfied, 
then the formulation of 6 is simplified. On the basis of 
6;4q, we formulate an equal 6/*qq in B and continue it 
one step to 6/q,,qy using the input value xy: Thus we 
have proven Theorem 2, 

Proof of Theorem 3, We have the arbitrary experi- 
ment 6,q;. We shall continue it in the machine IT 
using the experiment 6, (i.e., we take the input values 
of 6, from II,) which distinguishes the state q,. For the 
experiment &, q, 6, we can indicate an equal experi- 
ment & q' 6, in the machine B which is equivalent to A. 
But since L, distinguishes the state q, in II,, we have 
q'=q;. | 

Thus in B there is an experiment 6 q, equal to 64). 
We shall now continue the experiment 6q, by means of 
the experiment 6, in the machine II, which distinguishes 
the state q,. The experiment 6,q,L, has an equal experi- 
ment 6,q'6, in the machine B. But & distinguishes q, in 
Il,. Therefore q' must be that state which became a3; 
ie., q'=q,. 

Thus, we have the experiment & q, which is equal 
to &\q;. We shall now take the experiment 6, q;. Rea- 
soning analogously, we find that for 6 q, there is an 
equal experiment 6,q,;. Therefore, having 64), we ob- 
tain an equal experiment 6,q,. Reasoning analogously 
with respect to 6,q., we obtain an equal experiment 














EaQr-y, etc. Therefore we immediately have r equal to 
the experiments 
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Proof of Theorem 4, We shall take an arbitrary ex- 
periment q 9%, in the machine A. For the machine B 
there is a set Q of states qj such that qj} 6, are equal to 
Ao%- When the experiment 6, is extended the set Q 
will be decreased; however, due to the condition in the 
theorem the set will not be vacant. Since the set Q is 
finite it is possible to find a continuation of qg6, such 
that the corresponding set Q, will not be diminished for 
any continuation of it. 

Assume the experiment q96, ends with the state q, 
and that equal experiments in B end with the state qi 
(i=1,2,...,k). Amy experiment q,é in the machine 
A has k equal experiments q,6 in the machine B, In 
order to verify this it is sufficient to study the experi- 
ments q,6 and qi & as continuations of the experiments 
Go%, and Le. Since the experiments qié are equal, we 
find that qj (i=1, 2,... , k) are equivalent states of 
the machine B; this contradicts the simplicity condition 
for machine B. Therefore the set Q, consists of one ele- 
ment qi. Thus, the state q, corresponds not to k states 
but only to one state q} for any experiment ¢. But if q, 
qi correspond to one another for any experiment, then 
the states following them (as parts of equivalent experi- 
ments q,6 and qj}é) will also correspond to one another. 
From the simplicity condition for the machines A and 
B it follows that the correspondence will be mutually 
single-valued. 

Assume that the states following q, form the set 
A*, and that the states following qi form the set B*. 
We shall show that A* and B® are identical sequential 
machines, Assume q, € A® and qi €B* correspond. Then 
from equality of the experiments qg6q3 an? q36q} and 
the simplicity of the machines A and B it follows that 
qs and qj also correspond. Thus we find that the mutu- 
ally single- valued correspondence established between 
the states A* and B® is maintained for any experiment, 

Therefore the machines A* and B® are identical, 

Corollary. If the machine A is closely coupled and 
is equivalent to B, then A is identical to a part of B. In 
fact, for closely coupled machines A* coincides with A, 
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A theory is provided for the operation of a magnetic modulator with crossed fields and a periodically changing 
inductance of the control winding. The basic conclusions arrived at in the work are confirmed experimentally. 


Paper [1) provides the theory and design of a mag- 
netic modulator with mutually perpendicular fields and 
a fundamental frequency output for an ideal case when 
under the effect of the current in the excitation winding 
W, (Fig. 1), the inductance of the control winding W, 
varies periodically between the values Lmax and Lyin= 
=0. Infact, a modulator with a 1000-2000 ferrocart 
core has a ratio of inductance L,,9, to Lyin amounting 
to 4-5 owing to the limited power of the source andnon- 
uniform cross section of the core for the exciting flux, 
thus leaving some of the core sections unsaturated, 
Hence, the operation of a magnetic modulator with fi- 
nite variations of the control winding We inductance is 
of considerable interest, 

Let us examine a steady state condition in a control 
circuit with the signal source connected and a rectangu- 
lar wave modulating current supplied, for instance, by a 
magnetic transitron generator through a large resistor, 

Let us find the relation between the amplitude of 
the first harmonic voltage at the output of the modulator 
and the ratio Lyjax/Lmin =m for a signal voltage E. 

Figure 2a shows the characteristic variations with 
time of the excitation current i,; Fig. 2b, those of the 
inductance L of winding W,; and Fig. 2c, d, those of i, 
flowing through the winding, At instant t=0 the induc- 
tance of the control winding rises instantaneously from 
Lmin *© Linax: It is known that with the change in the 
inductance the value of the flux linkage # cannot 











change immediately (otherwise the control circuit would 
have an infinitely large emf induced in it [2]). Hence, 
if we denote by I,;, the value of the current in winding 
We of inductance Lmin before the instantaneous change, 
the value of current Iq at the end of the instantaneous 
change of the inductance to Lmax will be represented by 


9 (0) = VacLinig = Jol wax = Const, (1) 


i.e., the current in the control winding will decrease in 
an inverse proportion to the inductance of the winding: 
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After current i, drops to the value of Ij, at t= 0 it 
will begin to rise exponentially with a time constant of 


7% t 


i, = Sie + (Io — Ine) (1 —e *)* 
(3) 


Os ot <a, 


where 


0,4nW2 poS 


L 
r=, R= R+M, a 


R Imax = 

c 
and Ig=E/R, is the direct current component in the con- 
trol circuit at the steady state condition, 

At the instant t="/w the inductance of the control 
winding is reduced instantaneously to Lmin; the control 
current,which has by then reached the value of I, in- 
stantaneously rises to Iz, and then falls exponentially 
with a time constant of 1, =1/m: 


ig = I+ (Io — In) (1 —e * (2), 


’ 


(4) 
nm < wt [2x. 

In a steady state at t=2%/w , the control current i,= 
= Ip 

It is obvious that by similarity with equation (2) it 
is possible to represent the ratio of currents I), and Igo 
at the beginning and end of the instantaneous variations 
of the inductance from Lyyjgx tO Lyin by 


In — Emax _ 


Te Toni ) 


By inserting in equations (3) and (4) respectively 
the values of t= "fw and t= 2%/w we obtain for currents 
Iy, and Ip, the expressions 


Ip Ip + o= ha e*), (6) 


Toy = 49 + Uo —In)(1—e **). (7) 


By solving equations (2), (5), (6), and (7) with 
respect to Iy9 and Ig we obtain 


Tio = I,K, (8) 
T20 = IT, Ki, (9) 


where 
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By inserting the values of Iy9 and Ig in equations (3) and 
(4) we have 





t 
hes = Jo — Jee ra — KR), 0< wt < 2x, (10) 


tm 
og = Jo — Jee "(1 —Ki), x< ot <2n. (1) 


I 
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For a voltage at the output of the modulator with a 
small internal impedance of the signal source we have 


Mm = i,,Ri = Toh: — Toh: (1 — K)e * 


0O< ota, (12) 


tm 


tna = fyi = Toi — FoR: (1 — Kiyo * 
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Table 1 











| Core Dext. Dint» | Ho. 

| number mm mm _ | gauss/oe 
; | 29.4 7.6 1200 
| 2 29.4 7.6 2260 














ee Cross section 

6, H..0e | T;,,°C Ho ee a 
udinal flux‘ 

oe a a ta 

2400 0.1 120 0.68 

2000 0.4 80 0.68 














The amplitude of the first harmonic of the output 
voltage is 


Uim = | Ut, + U2, 


where 
Qn 
Via = . \ Uy Sin @t dwt, 
0 
2x 
Uy» = 4 \ Uy Cos wt dat. 


0 


After integration we obtain 


ER a : 
Uw = — Felten (i te “)—K)- 
— ae “Kite %*)). 


ER — 
Je = —r laa! +e *) (4—K) — 


pe. San or (1—K,)\(iten =) 


m?+ wt? 


Voltage Uym rises monotonically with an increas- 
ing m and wr, tending to a definite limit. For wr>10 
it is possible to neglect component Uy, of the output 
voltage and obtain for the amplitude of the first har- 
monic the following expression: 


s 2% 
a ate rearr ( an +e: — 


2,2 _2m , (14) 
Soo ee (1 — Ki) (te =]: 
If wr ©, we shall obtain for the modulator output 
voltage» 
(15) 


If m also tends to infinity we obtain from (15) the 
expression 


Umm = Usama ig =1.27E 3+. ae) 
c 


Let us note that for wt ~ © and m ~ © the coef- 
ficient K > 0 and Kj ~ 2. This limiting operating con- 
dition of the modulator corresponds to a square wave 


current in the control winding (Fig. 2d), whose amplitude 
is twice the direct component Ip = E/R,, of the current. 


In fact, the magnetic energy stored in the longitudinal 
field of the modulator for L = L,,,,, is completely dis- 
sipated in the load in the following interval, when 
L= Lmin. thus doubling the load current compared with 
the direct component of the current in the steady state. 
Figures 3 and 4 show the relations of coefficients K 
and K, to m for wr values of 10, 25, 50, and 100. It 
will be seen that K is relatively independent of wr. 


The relation of the first harmonic amplitude to m, 
calculated from formula (14) for various values of wr , 
is shown in Fig. 5. It will be seen that with a decreasing 
mandwt the amplitude of the first harmonic voltage at 
the output decreases, but even for wr = 25 and m=5 it 
remains equal to half its maximum possible value, and it 
is only 25% below Uj which corresponds to wr, 
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Table 2. 





















































m 
wr | Re: 2 | 3 | 4 | 5 Note 
___ | (ohms) | calc, lexpt. | calc.| expt. | calc. | expt. |calc. lexpt. 
Modulator 
50 ~ A .6 é é 
on 0.4 0.4 0.01 | 0.6 0.11 | 0.7 | with core No. 1 
100 | 1000 | 0.43 | 0.42 | 0.63 | 0.6 | 0.74 | 0.72 | 
60 | 2000 | 0.4 0.4 | 0.64 | 0.6 0.72 | 0.7 | 0.76 | 0.74 Modulator 
ne with core 
120 | 1000 | 0.43 | 0.42 | 0,63 | 0.63 | 0.75 | 0.74 | 0.83 | 0.831 No.2 
! 









































Tests of several modulators with two different types 
of cores, whose characteristics are given in Table 1, 
have fully confirmed the theoretical conclusions ob- 
tained above. 

Figure 6 shows the experimentally obtained charac- 
teristic variations in the inductance of the control wind- 
ing (W, = 2800 turns, Rg= 200) of the modulator with 
cores No, 1 and No. 2 when the excitation circuit (W, = 
= 200 turns) is fed by a direct current, These character- 
istics were used for calculating the transfer constants of 
the modulator when excited with a square-wave current. 

Table 2 gives the computed and experimental 
values of the transfer constant Uyp/E for the above mod- 
ulators at 1000 cps for various values of wr and m, 

It should be noted that not only the transfer con- 
stants but also the curves of the experimental and theo- 
retical current ic variations in the control winding coin- 


cide, For instance,Fig. 7 shows an oscillogram of the ex- 
citation current i, and the control circuit current i, for 
wt 50 and m=5, The points shown on the curve i, = 

= f(t) have been calculated from formulas (10) and (11) 
for the above values of wr and m, for which Iy9=0.3lg, 
Igq= 1. 71g, Ig, = 0.341g, and Ig, = 1.519, It will be seen that 
the shape of the curve for the current in the control wind- 
ing, and hence for the output voltage, does not differ 
greatly from the shape of the excitation current for large 
values of wr, and in the half-period corresponding to L= 


=Lmax it changes very little. 


CONCLUSIONS 

The theory of a magnetic modulator with a periodi- 
cally changing inductance in the input signal circuit has 
been derived, It has been shown that a sufficiently high 
transfer constant can be obtained with a relatively small 
variation of the inductance, 

For this purpose it is necessary to make the maxi- 
mum value of wr in the input circuit of the modulator 
sufficiently high (wr= 25), Under this condition the 
shape of the output voltage wave resembles closely 
that of the excitation current, 

The basic conclusions hold both for modulators with 
mutually perpendicular fields and for those of the same 
type but with the signal and excitation fields connected 
in parallel, since a parallel excitation field produced by 
a half-wave rectifier will also vary the inductance of 


winding We periodically from Lyjax tO Linin: 
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The conditions for obtaining increased efficiency of reversive magnetic amplifiers with dc output are determined. 
New circuits for reversive amplifiers satisfying these conditions are presented; among these are some that have 


increased efficiency without using transformers. 


Introduction 

The use of dynamoelectric amplifiers in automatic 
control systems presents numerous difficulties; the most 
substantial of these difficulties are their lack of dependa- 
bility in operation, the fact that they are complicated to 
adjust, their need of periodic maintenance, etc, There- 
fore, a tendency has lately developed to replace dyna- 
moelectric amplifiers by magnetic amplifiers, The use 
of magnetic amplifiers makes it possible to create con- 
trol systems which require almost no adjustment, and 
which are distinguished by high technical and economic 
merits, In the creation of such systems, a considerable 
part is played by the right choice of the magnetic am- 
plifier, by the correct approach to the design of the sys- 
tem, and by the using of advantageous design data, 

Any reversive full-wave magnetic amplifier is com- 
posed of two, or, in the general case, of several half- 
cycle (half-wave) amplifiers. 

There are two basic principles which may be used 
in the design of full-wave, two-semicycle amplifiers: 
the principle of summation of voltages obtained at the 
outputs of the semicycle amplifiers and the principle of 
summation of the output currents, A scheme based on 
the first principle is shown diagrammatically in Fig. 1. 
In this scheme, separate resistors Rj, serve as load resist- 
ances for each of the two amplifiers. To obtain a maxi- 
mum of power in the output when this scheme is used, 
the total common load resistance R,, must be equal to 
YZ R,. Amplifiers built according to this principle have 
very low efficiency; therefore, their application is 
limited to the area of low output power [1]. 

The diagram of a scheme based on algebraic sum- 
mation of currents is shown in Fig. 2. When the signal is 
applied to the control winding, the current in the output 
of one of the component amplifiers is increased; in the 
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output of the other, decreased. To obtain maximum 
power output, it is necessary that the output current of 
the second component amplifier have a minimum value 
when the output current of the first is at a maximum. 
To achieve this, the two component half-wave ampli- 
fiers must be made to work in the combined full-wave 
amplifier in the same way as when they operate sepa- 
rately, 

Direct subtraction of the two direct currents ob- 
tained from the output of the half-wave amplifiers is 
impossible because of the mutual interaction between 
the two amplifiers when they are connected with each 
other as shown, Usually this interaction is decreased by 
insertion of ballast resistances R, (shown in Fig. 2 in dot- 
ted lines), But these ballast resistances greatly decrease 
the efficiency. Instead of ballast resistances, control 
valves can be used inthis scheme [2]. In our paper revers-- 
ive amplifiers are reviewed which do not contain such 
valves, 


Basic Principles for Design of More Effi- 
cient Reversive Magnetic Amplifiers 

Let us discuss the character of the interaction of the 
two amplifiers in Fig. 2 in more detail, To present a 
clearer picture of the conditions, we shall assume that 








* In the original Russian and in the author's abstract in 
English, the term “reversive" is used, Inasmuch as this 
term is unfamiliar in the USA, it has been replaced in 
the title by “duodirectional,” which is the shortest equiv- 
alent term for this concept accepted in authoritative 
texts of the English-speaking world; see, e.g., Magnetic- 
Amplifier Circuits by William A, Geyger, 2d edit., 
McGraw-Hill, New York, 1957, p.77. In the text, how- 
ever, the word "reversive” has been kept as the shorter 
of the two terms [Publisher‘s note]. 
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for the half-wave amplifications single-core magnetic 
amplifiers are used; a circuit diagram of these is shown 
in Fig. 3. As will be made evident in the forthcoming 
discussion, this amplification scheme serves as the basic 
element for construction of more complex schemes of 
magnetic amplifiers, Before going to the examination 
of the basic circuits for connection of half-wave ampli- 
fiers into reversive combinations for dc amplification, 
let us point out some specific features in the operation 
of the scheme of Fig. 3 which must be taken into consid- 
eration in the design of such reversive amplifiers. 


‘ | 
fy Ae | Ie 
uh, 
Fig. 3. 


Let us assume that the core in the device of Fig. 3 
has a rectangular hysteresis loop (Fig. 4). In order to ob- 
tain the lowest possible value of the load current in this 
scheme—Iioaq min ~ it is necessary as a rule to have the 
following condition fulfilled; 

Tie 
\ udt < (B, — By) w,S, (1) 


0 


where u is the voltage applied during the active (output 
contributing) half-period to the operational winding Wpi 
By, the value of the magnetic induction in the core at 
the end of the control half-period (at the start of the 
active half-period); B,, the induction at saturation; and 
T, the period of the supply voltage. In the specific case 
when By= -B;, condition (1) assumes the following form: 


T/2 
\ udt < 2wpSB,. (2) 





When condition (1) is fulfilled, the core either does 
not reach the saturation state at all, or reaches it only 
at the end of the active half-period when wt=", when 
condition (1) is not fulfilled, then the induction of the 
core becomes equal to the saturation value B, at some 
moment during the active half-period which is deter- 
mined by the expression 


B,—B 
ot = a = are cos | Rie. |<Sxbe wp | < %. (3) 


U,, 
During the time interval when a<wt=<m there will be a 
current i= u/R going through the load, where R is the 
complete resistance of the load circuit, In formula (3) 
Um is the amplitude value of the supply voltage, and 
w=2"f, its angular frequency, In the half-period 
scheme of Fig. 3 by itself there is no difficulty whatso- 
ever in fulfilling conditions (1) or (2), and in obtaining 
a high ratio of high to low load current values, How- 
ever, when two or more half-wave amplifiers of Fig. 3 
are combined into a reversive scheme, condition (1) is 
usually violated, because of an increase of the left side 
of the inequality, and/or because of a decrease of the 
right side, 

An increase of the left part of inequality (1) takes 
place as a result of a redistribution of voltage in the 
windings of the reversive scheme which are in the active 
state during a given half-wave period, A decrease of 
the right side is caused by a redistribution of voitage in 
the active circuit while it is in the control exercising 
half- period. 




















Fig. 4. 


Let us illustrate the nature of the voltage redistribu- 
tions and the consequences thereof by discussing, for 
example, these phenomena in the most simple full-wave 
magnetic amplifiers (Fig. 5). Let us examine first the 
scheme of Fig. 5a, which represents a differential con- 
nection of two half-wave amplifiers the active half- 
periods of which are shifted in phase by a half-period. 
The instantaneous polarities of the supply voltage for 
the separate amplifiers are shown in the figure by aster- 
isks, When no signal current is present, the following 
current relationship is in effect; 


i) = —i(t+D) (4) 








































































































Fig. 5, 


and the dc component is absent from the output load. 
Let us assume now that to the amplifier input a current 
signal I, is applied which is of such magnitude and po- 
larity that core II is saturated during the entire active 
half-period, andthe current i,=u/R is flowing. In order 
to obtain the maximum intensity of current in the out- 
put from the amplifier, the current i, must have the 
least possible value, that is,the core I must not get satu- 
rated before the end of its working half-period. This 
can be achieved only if, during the control half-period, 
an induction By will be established which will satisfy in- 
equality (1). However, in the scheme of Fig, 5a, this 
condition cannot be fulfilled. Indeed, if core II, during 
its active half-period, is saturated, then only a small 
part of the supply voltage will be applied to rectifier 
valve V, which is equal to 


Urey = u— i,Ry =, (Ry + R,), (5) 


where Rg and R, are the resistances of the winding 
and of the rectifier valve in the forward direction. 

This voltage determines the highest rate of change 
of the magnetic induction in core I that can possibly be 
produced by the signal Iy during the control half- periods 
indeed, when the condition 


wS dB, i (6) 


a “rev 


is reached, valve V, becomes conducting, and the cir- 
cuit loop formed by the valves and operational (active) 
Wp which has a rather low total resistance, practically 
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shortcircuits the control winding wy of core I, because of 
the transformer coupling between the windings Wy and 
Wp of core I. 

The largest possible change of induction in core I 
during the control half-period is 


AB = B, — By= =~ \ weydt = 
7 - 


: (7) 





Ry 
“wv (!— FFRER, | 
2fw S 
where Ugy is the average value of the supply voltage 
and f—its frequency, 

Thus, By for core I will differ but little from B,, 
and during the next half-period the mean value of cur- 
rent i, will differ only very little from the mean value 
of current i,, regardless of the magnitude of the signal. 
Because of this, the maximum value of the constant 
component of the load current in the scheme of Fig. 5a 
will be much smaller than for the nonreversive half- 
wave scheme of Fig, 3 by itself. 

Let us now consider the scheme of Fig. 5b. Here, 
in distinction from scheme 5a, the operational periods 
for the two cores are not shifted in time, but coincide, 
Let us assume that the magnitude and the polarity of the 





Signals are such,that at the end of the control half-period 


we have induction By] = —B, in core I, and in core I 
the induction = +B,. During the entire operational 
half=period core II is saturated, and to the opera- 











tional winding wp of core I the following voltage is ap- 
plied: 


Ug =u + [u— i,(R, + Ry)] = 2u. (8) 


Therefore, to comply with condition (1) or (2), for 
the scheme of Fig. 5b a supply voltage must be chosen 
which is one-half of the one used for the scheme of 
Fig.3. Otherwise, core I also would reach the saturation 
stage approximately in the middle of its operational half- 
period, and this would lead to appearance of a consider- 
able i, current, which, going through the valves and the 
Wp windings, can be much larger than the load current 
ipy 

Thus, we come to the conclusion that in order to 
keep the working conditions of a half-wave amplifier 
unaltered when connecting it with other similar ampli- 
fiers into a reversive scheme, we must remove the effect 
of the load voltage upon the operation of the half-wave 
amplifier which must have a minimum of output current 
at the time. This effect can be avoided by the use of 
the well-known bridge scheme (Fig. 6) of subtraction of 


ci iia 


~ 
‘ + + 8, - 










































































the output currents from two half-wave amplifiers. The 
device is built up from four controlled rectifying ampli- 
fier elements. As such, single-core half-wave magnetic 
amplifiers such as shown in Fig. 3 may be used, or any 
known combinations of part cycle amplifiers, such as 
two-half-wave or multiphase amplifiers, Let the assem- 
bly be so arranged that the rectified current in the out- 
put of elements B, and B; increases, and the rectified 





current in the output of elements B, and B, decreases 
when the control signal is in action. The presence of 
ballast resistance Rj, which is equal in value to the load 
resistance R,,; prevents redistribution of voltage when 
the signal is on, Indeed, the voltage drops across Ry, and 
Rj; are applied to Bz and Bz, with opposite polarities, 
This results in the fact that the operating conditions of 
elements B, and B3, when included in this combined 
operation system, remain practically the same as wher 
they are used separately. Therefore, the current in the 
output of these elements can actually reach minimum 
values at the proper time, and thus assure maximum 
output in the load circuit of the assembiy, in distinction 
from the scheme of Fig, 2, the system of Fig. 6 has a 
much higher efficiency, which reaches 50% in the ab- 
sence of losses in the windings and in the rectifiers, 
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If in the system of Fig. 6 each of the four elements 
By, By, Bg, and Bs consists of a single-core amplifier 
(Fig. 3), then the assembly is the well known bridge 
scheme of reversive full-wave amplification by means of 
half-wave components (Fig. 7a). The scheme of Fig. 7a 
may easily be simplified by using a common source of 
power supply (Fig. 7b) and employing two amplifiers in- 
stead of four; the latter is possible since the cores I and 
Il and the cores III and IV work in identical regimes, 
and it is therefore possible to have only two cores with 
two operational windings on each core, and connect the 
four windings in the bridge scheme (Fig. 7b). This is pos- 
sible to do in every full- wave scheme built up from any 
half-wave elements, 

When, for instance, the component magnetic ampli- 
fiers (controlled rectifying elements) in Fig. 6 consist of 
twin half- period devices with center tap connections, we 
have the amplifier scheme of Fig. 8a. And this scheme 
also can be simplified by the use of common cores for 
pairs of two operational windings through which the same 
currents must pass (as e.g., the windings w,, and w,,, and 

and w,,), and also by decreasing the number of sec- 
ondary windings of the supply transformer; the latter is 
achieved by means of connecting together the points in 
the circuit which are at equal potentials. For this pur- 
pose, the scheme of Fig, 8a is redrawn as showa in Fig.8b, 
wherein those ends of the secondary windings which at a 
given moment are at a positive potential are shown by 
asterisks, From Fig. 8b it is evident that the two second- 
ary windings w, and w,, and the two secondary windings 
w, and w; can be replaced each by a single winding. As 
a result, the scheme of Fig. 8a is transformed into Fig.8ct 

Other magnetic amplifiers with reversive or duodi- 
rectional operation may be constructed in a similar man- 
ner, 


Dependence of Efficiency upon Various 
Factors 

We shall examine the basic factors which affect the 
efficiency of magnetic amplifiers built according to the 
scheme of Fig.6. The dc power output in the load is 
equal to 





Py = Ry = (Lio ee T4)* Ru, (9) 


where Iyp and Igg are the dc components of the output 
currents from the corresponding half-wave amplifiers. 
The effective (root-mean-square) values of these cur- 
rents are respectively 1, = ke; Iy9, Ip =kfgIgq, and Ip; o¢¢= 
= kg Typ 

Let us assume for the sake of simplicity that kg, = 
=kg=ky. Then the entire powerconsumed in the cir- 
cuit is approximately equal to 


P, = (hh + /2)? Rb + 
+(h — In)* Ry + 22R + 12K, 


(10) 








where R is the resistance of one operational winding to- 
gether with the valve connected in series with it. From 
(9) and (10) we find the following expression for the ef- 
ficiency of the amplifier: 





Py 
y P, = 
(7; — J)? Ry (11) 





KB UC + 1)? Ry + (i — a? Ry + (B+ 72) 2R} 


We introduce the following denotations; 1) the ratio 
of the currents in the windings of the two half-wave am- 
plifiers k=1,/,, and 2) the ratio of the resistance of one 
operational winding and the valve associated with it to 
the load resistance m=R/Ry;. After introduction of these 
denotations into (11) and assuming that Ry =R,;, we have 


0-4) 


2K (1 + m) (1 aa =r) ) 





— (12) 


Values of efficiencies calculated by (12) for various 
reversive amplifiers built according to the scheme of 
Fig. 6 are given in Tablel below. The efficiency data 
in this table are all computed for the same resistances 
ratio m=0,1, but for a number of different current ratios 
k of the half-wave component amplifiers. The data in 
Table 1 are arranged in three columns; a) for the case 
when each arm of the bridge in Fig. 6 contains a single- 
half-period single-phase magnetic amplifier according 
to Fig. 7; b) for the case when the bridge arms are made 
up of twin half-period single-phase amplifiers according 
to Fig. 8; and c) for the case of three-phase half- period 
amplifiers according to Figs. 9 and 10 in the bridge arms, 


TABLE 1. Efficiency values of reversive mag- 
netic amplifiers with resistive loading and 














m = 0,1 
Amplifier in bridge arms 

Single | Two half- Two half- 
half- | periods periods 

k period jsingle-phase | three-phase 
k¢=157| kg =1,11 k¢ #1.00 

5 0,11 0,23 0.28 

10 0.15 0.30 0.36 

20 0,17 0.33 0.41 

a 0,19 0.37 0.46 














Tit should be emphasized that the output leads from the 
transformers in the schemes of Fig. 8 are connected to 
center taps of the secondary windings, and that in both 
halves of the secondary windings on either side of the 
center taps there are always equal potentials U.. 
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Fig. 9. 








The values for k¢ given in the table all correspond to the 
maximum value of the load current and to a currents 
ratio k=, 

Reversive schemes of magnetic amplifiers using 
three-phase elements can be built on the same principles 
as those described above, In Fig. 9a a scheme of a re- 
versive amplifier is shown which is composed of a three- 
core magnetic amplifier with internal feedback. Com- 
paring Figs. 9a and 9b, it is evident that by combining 
each two secondary transformer windings of Fig. 9a desig- 
nated by the same symbols Wy Wo, etc., a simplified 
version of the scheme (Fig. 9b) is obtained, For correct 
operation of the amplifier, the voltages across the trans- 
former windings w+4, Wy,, and w-y, must be shifted in 
phase by 180° with respect to the voltages of the corre- 
sponding transformer windings w'y,, W'y,, and wy. The 
entire amplifier is built on six cores. If the three-core 
amplifier in the scheme of Fig. 9 is replaced by a 
three-phase bridge magnetic amplifier composed of six 
cores, the result is the scheme shown in Fig. 10 [3]. 




















Fig. 10. 


Two-half-period Schemes of Reversive 
Magnetic Amplifiers of High Efficiency 
Without Transformers 

In the schemes of Figs, 8, 9, and 10, a supply trans- 
former with a power rating exceeding 2-3 times the load 
power (see Table 2) is an inevitable part of the device. 

However, utilizing the principles for designing re- 
vers.ve amplifier schemes stated above, it is possible to 
construct schemes of two-half-period reversive ampli- 
fiers without the use of supply transformers, It would 
appear that such an amplifier could be constructed by 
combining two single-half-period reversive amplifiers 
built according to the scheme of Fig. 7b and delivering 
during alternating half-periods their outputs to a com- 
mon load, One of the feasible schemes of such a com~- 
bination of single-half- period amplifiers is shown in 
Fig. 11a. In this scheme a common ballast resistance 
Ro= Ry is connected in series with the power supply. 
However, this scheme is inoperative, for this addition 
of resistance Rg eliminates the redistribution of the sup- 
ply voltage during the active half-period only in the 
windings of an unsaturated core after saturation of one of 
the other cores, However, if one of the cores (for in- 
stance, I) is saturated during the entire working half- 
period, then the reverse supply voltage on the valves 
connected in series with the operational windings of an- 
other core (IV),which is at the time in the control half- 
period, will be close to zero; and that will be so, specifi- 
cally, also for the core which, for efficient operation, 
should have a minimum of current in its windings at that 
time. As was pointed out above, this lowering of the re- 
verse voltage on the valves leads to a considerable in- 
crease of the current through the operational windings in 
the next active half-period of operation; and thus the 
current in the output load is considerably lowered, 

Better operational characteristics are obtained by a 
two-half-period reversive scheme suggested by the au- 
thors which is shown in Fig. 11b. 

As pointed out above, to be sure of high effic 
in any reversive magnetic amplifier, the following two 
conditions must be fulfilled; 1) saturation of one or sev- 
eral cores during a half-period must not lead to an in- 
crease of the voltage applied to the terminals of the 
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operational windings of the unsaturated cores during that 
half-period; 2) presence of a signal at the input of the 
amplifier, and, hence, of a voltage in the output load 
must not lead to a decrease of the reverse voltage on the 
valves in the regulatory half- period. 

Satisfaction of both above conditions is achieved in 
the scheme of Fig, 11b by adding, prior to the entrance 
of the supply voltage to each of the two apexes of the 
two bridges of Fig, 11a, a resistance R,/2 and a capaci- 
tance C in parallel with each of these resistors, The 
value of Ry is chosen to be equal to one-half of the load 
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resistance Rj. As a result of this, when, for instance, 
cores I and II are saturated, the voltages on the opera- 
tional windings of cores III and IV remain unaltered, be- 
cause during each half-period one-half of the supply 
voltage is absorbed in the Rg resistances connected at 
opposite ends of the bridge. The shunting ofthe resistors 
by the capacitors C produces the result that during the 
regulatory half-periods the reverse voltages on the valves, 
connected, for instancc, in series with windings wf, 

Wily “fy and wy, also are practically not decreased at 
all, evenwhen cores I and II are fully saturated during 
the working half-period. The magnitude of capacitance 
C is not critical; it may be chosen to equal C»(I-2)/R,f), 
where f is the frequency of the supply voltage, Electro- 
lytic condensers may be used for this purpose, since, first, 
the voltages across the resistors Re/2 are of constant po- 
larity as shown in Fig. 11b, and, second, leakage currents 
in the condensers have no substantial effect on the opera- 
tion of the amplifier. 

A warning, however, must be given against the 
choice of condensers of capacitance values so high that 
the voltage across them during the entire period of the 
supply voltage will practically be equal to the full am- 
plitude variation of voltage across Rg at C= 0: for, since 
the voltage across the condensers is subtracted from the 
supply voltage, with such very high capacitances an out- 
put voltage from the amplifier is observed which is some- 
what lower than the maximum possible, 

In place of the scheme of Fig, 11b, the scheme of 
Fig. 1lc can also be used, which makes it possible to 
satisfy the above indicated condition with fewer resistors 
and condensers, This scheme, howéver, has the draw- 
back that the reverse voltage on the valves may reach 
higher values than in scheme 11b. Let us assume, for 
example, that the magnitude and polarity of the signals 
is such that cores I and II are completely saturated dur- 
ing the respective working half-periods, The peak in- 
verse voltages on the valves connected in series will then 
be 4 U;,, for the scheme of Fig. 11b, and U,, for the 
scheme of Fig. 1lc, where Upp, is the value of the supply 
voltage amplitude, 
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Scheme Uy T, ore _4 e U, I, P, # Py Pe n 
a | ae fem | a | a | ee ma 
of: H H " 4 " H n u u , 
144d | 0.781 1.74] 6.94 | 3.44] 0.78 | 2.45 | 19.6 | 2.96 | 0.406 
| Fig. 8 
eae <-* was = eles Eu 
0.855] 0.59 | 4 6 2.4 | 0.59 | 41.24 | 14.9 | 2.74 | 0,5 
Fig. 9 
10.37 | 0.82 [0.43 | 5.45 | 1.05 | 0.58 | 0.61 | 44.6 | 2.4 | 0.5 
Fig. 10 










































A consideration of the actual working conditions 
will decide which of the above schemes is the most ad- 
vantageous ia every specific case, 

In Fig. 12, characteristic curves are given for ampli- 
fiers constructed according to the scheme of Fig, lla 
(curve 1), according to the scheme of 11b with ripple- 
smoothening condensers (curve 2), and according to the 
same scheme without ripple-smoothening (curve 3), The 
differences in characteristics evident in Fig, 12 will be- 
come still larger when there is an opposing emf or an in- 
ductance in the amplifier output load, 

Experimental comparisons of the proposed schemes 
with all other known schemes of reversive half-period 
amplifiers with dc output have shown that the proposed 
schemes not only make it possible to operate without a 
power transformer, but also have the highest efficiency, 
which approaches 0,4 when losses in windings and valves 
are small. 


Comparison of Various Schemes 

Selection of one or the other scheme of a reversive 
magnetic amplifier for application in a specific control 
system must be made after all the technical and eco- 
nomic aspects of the feasible modifications have been 
discussed and the most advantageous"decided upon. 

We present above in Table 2 a number of data, 
characteristic of the operation of our schemes shown in 
Figs.8, 9, and 10, The testing of these schemes was 
carried out with purely ohmic resistances in the loads; 
the remanent inductance of the operational windings and 
the leakage inductances of the transformer as well as 
losses in the windings and valves were neglected in these 
experiments, All data given in this table are relative— 
with reference to the average values of the amplifier 
output characteristics I,; (current), Uj; (voltage), and 
Py; (power). 

The rated core power Poore, given in the table, is 
the product of the maximum rms values of the voltage 
in the operating winding Up: the current through it Ip 
and the number of operating windings n on the core, 
The rated rectifier power P, is understood to equal the 



































product of the peak inverse voltage Up at the valve with 
the maximum rms value of the current Ip through the 
valve, , 

In the above table ng designates the number of 
cores, and np the number of valves used in the scheme. 
Prr is the rated power of the transformer. 


CONCLUSIONS 

On the basis of our investigation the following con- 
clusions can be drawn: 

1) For construction of reversive (duodirectional) 
schemes of magnetic amplifiers with dc output operating 
with higher efficiencies, it is advantageous to use a 
scheme based on subtraction of currents of nonreversive 
amplifiers, | 

2) To obtain increased efficiency in these schemes 
the following conditions must be fulfilled; 

a) saturation of one or of several cores must not 
lead to increase of the voltage at the terminals of the 
working windings of any of the nonsaturated cores during 
the working half-period; 

b) the voltage on the load must not produce a de- 
crease of the inverse (rectifier) valve voltage during the 
control half- period, 

3) Fulfillment of condition 2) can be assured by con- ; 
nection of half-wave amplifiers into the bridge scheme 
of Fig. 6. 

4) To obtain higher efficiency, it is necessary in ad- 
dition to fulfillment of conditions 2) to assure a high 
ratio of high to low current in the working windings of 
the nonreversive amplifiers and small pulsations in the 
rectified current. Therefore, all other conditions being 
equal, reversive amplifiers using three-phase current 
have the highest efficiency, 

5) In this paper new schemes for reversive magnetic 
amplifiers with dc outputs without the use of transformers 
have been proposed; they satisfy conditions 2) and have 
high efficiencies, 
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Electrical circuits are reviewed in which transistors and two position vibrators are used for voltage regulation of 
dc generators, Arrangements are presented which provide increase of the armatures frequency of vibration , as 
well as schemes using additional transformers or choke coils to secure better performance of the regulator at 
elevated ambient temperatures, Results of experimental investigations are given, 


Vibratory electromagnetic relays and regulators 
have found wide application in various technical fields 
[1]; for instance, voltage regulation of modern dc and ac 
generators of automobiles is carried out exclusively by 
means of vibrating regulators [2]. 

A serious shortcoming of vibratory regulators is the 
presence of a continuously vibrating switching system 
which all the time breaks and makes contacts to open 
and close an inductive current-carrying electrical cir- 
cuit. 

It is known that the wearing out of contacts is due 
to mechanical, chemical, and electrical phenomena. 
The greatest wear is caused by electrical phenomena oc- 
curring when an arc is produced during the breaking of 
contacts, However, when the current in the circuit of 
the contacts is lower than 0,05-0,06 amp, and the volt- 
age applied to them is less than the breakdown voltage 
of the air gap at the contact, then there is practically no 
erosion taking place when contacts are broken [1]. Such 
conditions for operation of contacts can be obtained by 
various schemes using semiconductors, 

In the present paper, electrical regulation schemes 
are reviewed which use transistors and two- position vi- 
brators for control of dc generators, 


1, Simplest Circuits for Voltage Regula- 
tion 

A number of schemes for voltage regulation of dc 
generators in which transistors and vibratory elements 
are employed were reviewed in [3], The vibrators re- 
viewed there consist of nonsymmetrical polarized two- 
position relays (one of the two contact positions being 
“predominant”), 

In actual operations two-position vibrational regu- 
lators having one pair of normally closed contacts proved 
the most dependable; our manufacturing industry has be- 
come quite familiar with the production techniques of 
such regulators, Therefore, it is most advantageous in 








schemes employing transistors also to use vibrational reg- 
ulators with one pair of contacts, 

One of the feasible schemes for voltage regulation 
by means of a vibratory device in combination with a 
transistor is shown in Fig, 1. In this figure, OB is the ex- 
citation winding of a generator, and OO and K are the 














aD z : 


Fig. 1. 





ic, amp 
20 7 























12 =a 



























































4 fj 
7 he ‘ 
a oe a ae oe 
0 a2 a4 a6 QU... Vv 





energizing coil and the contacts of the vibrator. Rep is 
a resistance connected as shown between the emitter 
and base electrodes of transistor T which serves the pur- 
pose of decreasing the potential difference between 
these electrodes when contacts K are broken. Semicon- 
ductor diode D serves to eliminate the effect of the self- 
induction emf in the generator's excitation winding upon 
the operation of the transistor; see [4] and [5]. The gen- 
erator excitation winding is in our figure connected in 
the collector loop of transistor T. However, it may be 
connected instead in the emitter circuit. The power 
being broken by the separation of the contacts is the 
same in both cases; therefore, in general, both of these 
methods of connection of the excitation winding are of 
equal value. 

However, in single conductor systems it is more ad- 
vantageous to have the generator excitation winding in 
the emitter circuit of the transistor: in this case, if the 
generator minus pole is connected to the housing, it 
will not be necessary to insulate the transistor housing 
from the generator housing. In this way the generator 
housing may be utilized as a thermal sink to carry away 
the heat from the transistor. Such schemes will be dis- 
cussed below. 

The scheme of Fig. 1 operates in the following man- 
ner. 

When the generator voltage is lower than the one 
to be maintained by the regulator, the contacts K are 
closed, and the transistor base is connected to the col- 
lector. We then have uch=0, where uch is the voltage 
between the collector and base. 

In Fig, 2 the volt-ampere curve of power transistor 
P 208 at uop=0 is presented; this curve can be expressed 
with an accuracy sufficient for practical purposes by 
the formula 


Yec * Ung +Rric, (1) 


where uec is the voltage between the emitter and col- 
lector of the transistor; i,, the collector current; up, 
the projected value of the emitter-collector voltage at 
ic = 0 (e-c voltage with collector circuit open). R, is 
the average dynamic resistance of the transistor, for 
which we have the expression 


Ry = tga = a ° 


For the particular characteristic curve of Fig. 2, ex- 
pression (1), after substitutions of the numerical values, 
assumes the following form: 


Use = 0,25 + 0,029 ig? 


Using approximation formula (1) it becomes possi- 
ble to obtain the simplest possible expressions for the 
voltages and currents which are characteristic for the 
processes taking place in the schemes under discussion. 
For example, when self-excitation takes place in the 
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generator, the voltage generated uc, is balanced by the 
sum of the voltage drops in the transistor and in the ex- 
citation winding; that is, we have 


UG = Uec+ int, , (2) 
where tp and i, are respectively the resistance of and 
the current in the excitation winding. 

Considering that ip =i, (Fig. 1), and substituting (1) 


and (2), we obtain the following expression for the gen- 
erator excitation current: , 


YUG~"ToO (3) 
th + Ry 


When the generator voltage u, reaches the standard 
output value, Contacts K are opened, Then we have 
Ueb * 0, since the base is usually connected with the 
emitter by a resistance Rep of low ohmage. In order to 
obtain simple relationships, let us assume that the in- 
ductances of the excitation winding of the generator and 
of the vibrator energizing coil are constant, and that the 
resistance of the transistor in the reverse direction is in- 
finite, Then for our scheme in its present condition the 
following equation will be valid; 


ip = 


Lg “8 + iptg+up = 0, (4) 
where uy, is the potential drop across the rectifier diode 
D, Lg~the inductance of the generator excitation wind- 
ing. 

Since in this state the resistance of the transistor is 
considerably higher than that of the diode which is now 
in the forward conductance state, it is possible to con- 
sider that as a first approximation the extra current 
caused by the closing of passage through the transistor is 
going through the diode loop, and that 


ip * ip: (5) 
Vibrating voltage regulators usually have high ratios 
of release currents to holding current reaching up to 095, 
By the time. the next contact make is due, the current in 
the diode usually has not had time to go down to zero. 
Therefore, the volt-ampere curve of the shunting diode 
may be expressed, with accuracy sufficient for practical 


purposes, by an experimental formula similar to (1), 
which is as follows: 


Up =U), + IpRp- (6) 
Solving equation (4) for the generator excitation 
current ip, taking into account (1), (5), and (6), we ob- 
tain the following formula for the excitation current 
during the time when the regulator contacts are open: 


- “ ret ?p e 
in = ing? eigeetlie tno thc gil 
r, + Rp , 





where ipo is the magnitude of the generator excitation 
current at the moment of opening of the contacts, 


























Expression (7) makes it possible in a relatively sim- 
ple way to obtain an approximate estimate of the aver- 
age value of the diode current, and to choose the proper 
type of diode. To do this it is sufficient to substitute 
into (7) the value of ip from (5), Also the time interval, 
must be known during which the contacts remain open, 
and the vibration frequency, Then for the average value 
of the diode current I, the following formula is obtained: 


fs fly 
D r, + Rp 
ratRkp 
; “po \ > as, “no 
| (em | 4 F Rp} (t—e B prone ty (8) 


where f . is the frequency of vibration of the contacts, 
and t,, is the time the contacts are open during a cycle, 

y the use of a shunting semiconductor diode not 
only is the proper operation of the transistor under the 
prevailing conditions of back and forth switching secured, 
but also the appearance of excessive voltage between the 
emitter and collector is excluded, 

When the contacts are open, the generator voltage 
u and the voltage across the diode are of the same po- 
larity; consequently 


Solving equations (6), (7), and (9) simultaneously, 
we obtain an expression for the voltage uec between the 
emitter and collector of the transistor at the time when 
contacts K are open: 


Uec = Ug+ Upo+ 


oo JTetRp | a rgtFD 
+- Rd ie bp — (4 -— € bp )] . 
(10) 








“DO 
a+ Rp 


The product of Rp with the expression in the brackets 
can amount to only a few percent of the voltage to be 
regulated, The same is true for the voltage u,. There- 
fore, one may consider that in schemes where the exci- 
tation windings are shunted by a semiconductor diode, 
the voltage between emitter and collector is approxi- 
mately equal to the voltage of the generator which is 
being regulated, The larger the resistance of the diode 
circuit loop, the higher is the emitter to collector volt- 
age, 

Formulas (9) and (10) make it possible to estimate 
the maximum voltage that can be expected on the tran- 
sistor in this application, and accordingly to choose the 
proper transistor for it, 

When a shunting diode is not present at the moment 
the contacts are opened, the voltage across the transistor 
may exceed the regulation voltage of the generator u, 
by a factor of several times [3]. 

It is the resistance Re} connected between the emit- 
ter and the base (Fig, 1) which determines the potential 






YecV 
Fig. 3. 

difference between the emitter and base when the con- 

tacts K are open, Thus, the resistance of the transistor 

in the nonconducting state is determined by the value 

of resistance Rep. 

In Fig. 3 are presented volt-ampere curves of tran- 
sistor P 208 for various values of resistance R,}, to illus- 
trate the effect of the shunting resistance upon the tran- 
sistor resistance in its nonconducting phase. [,; in Fig,3 
is the initial collector current. The characteristic curves 
of Fig. 3 were taken at an ambient temperature of 50°C, 
The measurements were made after a two-hour warmup 
of the transistor. The surface area of the heat sink, which 
was made of aluminum, was 750 cm’, 


Voltage regulator schemes such as those shown in 
Fig. 1 have two substantial shortcomings: 1) Relatively 
low speed of the regulation process, and large fluctua- 
tions in the voltage being regulated, both of these being 
due to the low frequency of armature vibrations; and 
2) the current through the contacts and the resistance of 
the transistor in its closed (nonconducting) condition are 
both determined by the magnitude of the shunting resist- 
ance Rep; the lower Rep, the larger the current that is 
interrupted by breaking of the contacts; but if Rep is in- 
creased, then the resistance of the transistor in its 
blocked condition is lowered, However, these short- 
comings can be easily eliminated by the use of special 
schemes, 

2. Methods for Raising the Workin 
quency in Voltage Regulators by 
of Transistors 

Increased vibration frequencies in voltage regulators 
are obtained by the use of additional windings and resist- 
ances (Fig. 4). 

In the scheme of Fig. 4a the frequency of the arma~ 
ture vibrations is raised by means of an acceleration 
winding YO which is connected in parallel with the 
generator excitation winding OB; in this and in the fol- 
lowing schemes the beginnings of windings are indicated 
in the figures by points, In order to compensate for the 
lowering of the resultant magnetic force with increasing 
speed of rotation of the generator armature, caused by 
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the presence of the acceleration winding, a compensat- 
ing winding KO is added which is wound in a sense to 
produce a magnetization counteracting the magnetiza- 
tion by the basic acceleration winding. 

When the contacts of the vibratory regulator are 
closed, the acceleration winding YO has applied to it 
the generator excitation voltage with the voltage on KO 
added; thus we have 


Uy = Up + UC= (Fa +7 odin, (11) 


where uy is the voltage on the acceleration winding; u,, 
the voltage on the generator excitation winding; uc, the 
voltage on the compensating winding; and r_,, the resist- 
ance of the compensating winding. 

When contacts K are open, the excitation winding is 
shortcircuited by the diode; the voltage on the accelera- 
tion winding then is equal to the potential drop across 
the diode. When the contacts are broken, the voltage on 
winding YO is decreased abruptly by 


Auy = (rp + 1c) in — Uy. (12) 


The change of current through winding YO can be 
expressed by the following simple formula (neglecting 
the mutual effects of windings OO and OY): 


. — fa t'obe Ty _»( Ly ) 
ee ry y * (13) 





























Fig. 4, Scheme of voltage regulation with acceleration of vibration 


of the regulator armature. 





where ly and Ly are the resistance and inductance of the 











acceleration winding. 
It follows from expression (13) that at the instant 






ty = *y In E gE ma] 
ry “un 
the current in the acceleration winding changes its direc- 
tion: when 0=<t=<t, the resultant magnetizing force is 
equal to the sum of the magnetizing forces produced by 
the principal and acceleration windings, and whent>t, 
—to their difference. 

Consequently, when the contacts are opened, the 
resulting magnetizing force of the regulator windings is 
substantially decreased, which leads to a considerable 
demagnetization of the vibrator core, and this in turn 
brings about a faster closure of the contacts [remember 
that the vibrator being of the “normally closed" type, 
the contact “make” is mostly caused by spring action on 
the armature opposed to the magnetic action= Translator], 


(14) 


If the "make" of the contacts occurs at time t, after 
the preceding break, acceleration winding YO will have 
at that moment the voltage uy applied to it, and the 
current in this winding will begin to rise; this change 
will take place according to the following expression: 


j — 2 tpt) 
patetoe(y_ Be) 4 


r. 

y 
ge 
e y 


ry 


*Y (4t) 
wg oo (15) 






































The increasing current in the acceleration winding 
will speed up the occurrence of the “break” of the vibra- 
tor contacts, 

In the circuit of Fig. 4b acceleration is achieved by 
a winding YO connected in series with the emitter of the 
transistor. When the contacts are closed, the current in 
this winding is determined by the expression 


— Fe u 2% t) 
ly = ge t'_“w(4_ " (16) 
where ip is the excitation current when the contacts are 


closed, 


R=tg +Ryp+ly +l: and L= La + by + boi 

here ty,lo, Ly, and LL, are the resistances and induct- 
ances of the acceleration winding and of the compensat- 
ing winding, connected now in series, 

The volt-ampere curve of the transistor in this case 
is again represented by formula (1), 

When the contacts K are opened, the generator ex- 
citation current is substantially going only by way of the 
diode circuit, while the winding YO is practically cur- 
rentless, Thus, the magnetizing force of the regulator 
winding is diminished, and thereby conditions for faster 
closing of the contacts are created, 

The frequency of the regulator armature's vibrations 
can be raised also by connecting an acceleration wind- 
ing in series with the shunting diode (Fig.4c). The ac- 
celeration winding YO is connected in opposition to the 
principal regulator winding OO. In the circuit of Fig.4a 
the winding YO must be of high resistance, in the cir- 
cuits of Figs. 4b and 4c—of low resistance. 


’ 



































In Fig. 4d a circuit for voltage regulation of a gener- 
ator is presented, in which a frequency increase of the 
vibrator armature is obtained by means of an additional 
“acceleration resistor” Ry with another resistor Ry con- 
nected in series with it. 

The entire resistance (Ry +Rp) is connected in shunt 
to the transistor, and the principal winding of the vibra- 
tor is now connected to the tap between Ry and Ro: in- 
stead of to the minus pole of the generator armature as 
before. The magnitude of the shunting resistance (Ry + 
+Rp) is determined by the excursions of the excitation 
current in the process of the voltage regulation of the 
generator. In voltage regulators used in automobiles, 
the Ry part of the auxiliary resistance serves not only for 
purposes of acceleration but also as a heat sink for the 
principal regulator winding OO. 

The voltage applied to the principal vibrator wind- 
ing when the contacts are closed is almost equal to the 
generator voltage, The voltage across the shunting resist- 
ance (R,,+ when the contacts are broken is deter- 
mined by (9) and (10) above, The voltage on the. vibra- 
tor winding OO is then equalto 


tee ye 
Ry + Rp : 


= (17) 


Obviously, the smaller Rp, the smaller the voltage 
on the OO winding at the moment the contacts are 
broken, and the faster the falling off of the current in it 
takes place, When Rp=0 the accelerating action of the 
OO winding in the scheme of Fig. 4d is similar to that of 
Fig. 4a. 



































Fig. 5. Scheme of voltage regulation with an additional blocking 


circuit. 





However, these schemes are not altogether equiva- 
lent to each other: in the scheme of Fig. 4d closing and 
opening of the contacts is affected by changes in the 
magnetizing force of only winding OO, while in the 
scheme of Fig. 4a both windings OO and YO are having 
an effect. Also, the magnetizing force produced in the 
OO winding varies only within relatively small limits in 
the process of regulation, The latter property is also 
shared by the schemes of Figs. 4b and 4c, 

The scheme using an additional acceleration resist- 
ance has,in comparison with the schemes using addition- 
al windings for acceleration, also the advantage that less 
space for accommodation of the windings is needed. This 
makes it possible to decrease the size of the regulator 
magnetic system and to save on copper, 

An analysis of this kind of scheme for raising the 
vibration frequency of the regulator armature is given 
in [2}. 

In the schemes of Fig. 4 discussed above the gener- 
ator excitation winding was connected in the emitter 
circuit of the transistor. However, in all these schemes 
the excitation winding may be just as well connected in 
the collector circuit, as was shown in Fig. 1. 


3. Voltage Regulation Schemes with an 
Additional Blocking Circuit 

In the voltage regulation schemes of Figs, 1 and 
4, the current load on the vibrator contacts is much 
lower than in the ordinary regulators used at the present 
time. However, depending upon the magnitude of the 
Rep resistance, the currents being interrupted by the con- 
tacts can still be of considerable magnitude. 

Using additional blocking circuits it is possible 
practically to eliminate the currents depending upon the 
presence of the Reb resistance, and, consequently, to un- 
load the contacts still more and practically exclude 
erosion, Such schemes are presented on Fig. 5; in these 
schemes instead of the resistor R.p a transformer Tf or a 
choke coil Ch is used in suitable places, 


In Fig. 5a the primary winding of the transformer 
connects the emitter and collector of the transistor, and 
the secondary winding, the emitter and the base. 

If the primary and secondary transformer windings 
Ww, and we are traversed by the same changing magnetic 
flux, the following emf's are induced in them: 








d® d® 
Ge Dyar » = — Sa 


from this it follows that 
bf 


ee We © 


Consequently, if the effect of the voltage drops in 
the windings is neglected, the voltage variations across 
the emitter and base when the contacts are opened and 
closed will be of the same character as those across the 
emitter and collector, When contacts K are closed the 
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voltage on the primary winding is determined by the 
transistor volt-ampere curve of Fig. 2; it is equal to the 
emitter-collector voltage corresponding to the collector 
voltage u,.=0, When the contacts are open the voltage 
on the primary transformer winding is determined by - 
formulas (9) and (10), The voltages applied between the 
emitter and base will be similar in shape. When con- 
tacts K are open, the positive pole of the potential, pro- 
duced by the emf induced in the secondary transformer 
winding, will be applied to the base, the negative—to 
the emitter, In this way passage of current through the 
transistor is reliably blocked. In the scheme of Fig. 5a 
the primary and secondary transformer windings must be 
of high ohmage. 

The scheme of Fig. 5b works in a similar manner, 
Here the primary transformer winding is connected in 
parallel to the generator excitation winding. In Fig. 5c 
a voltage regulation scheme is shown having a low ohm- 
age primary transformer winding connected in series 
with the diode. If the variation of the voltage drop 
across the diode with current is neglected, the shape of 
the voltage induced in the secondary winding will be the 
same as that induced in the scheme of Fig. 5a, In the 
scheme of Fig. 5d, a positive potentiai on the transistor 
base relative to the emitter is secured by the use of 
choke coil Ch connected between the plus pole of the 
generator and the transistor base. 

In order for the induced emf to last longer, the time 
constant of the choke circuit must be increased, This is 
easily accomplished by connecting another semiconduc- 
tor dicde D, in parallel with the choke. In the scheme 
of Fig. 5d, the diode D, serves also to limit the magni- 
tude of the applied signal, 

To decrease the power handied by the contacts when 
a break occurs, the secondary windings of the transform- 
ers and the choke winding must be of high ohmage: if 
the currents from these windings are made sufficiently 
small so they can be neglected, then there will be only 
the load of the transistor base current on the contacts, 

The voltage regulation schemes with a supplemen- 
tary blocking circuit shown in Figs, 5a, 5b, and 5c can 
be realized also by connecting the generator excitation 
winding in the collector circuit loop of the transistor. 

The reviewed voltage regulation schemes with sup- 
plementary blocking circuit also insure dependable opera- 
tion of the entire regulator at elevated ambient temper- 
atures, 
In Fig. 6* the relationships between the emitter- 
collector voltage of transistor P 208 and the collector out- 
put current are given for different emitter-base voltages, 

The heat sink was made of aluminum, the active 
heat absorbing area was 750 cm? at an ambient temper- 
ature of 85°C; this is the highest permissible ambient 


* The curves of Fig, 6 are for the same transistor 
specimen as those of Fig. 3. 
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temperature when the collector of the transistor is loaded 
to the maximum rated power output, The voltages be- 
tween the emitter and base were applied from a source 
which had a constant output resistance, 

With increase of positive potentiat on the base rela- 
tive to the emitter, the collector current first decreases, 
At a potential between 0.2 and 1,2 v the collector cur- 
rent reaches a minimum. With further increase of posi- 
tive potential on the base, the collector current some- 
what increases, Hence, when P 208 type transistors are 
used, the transformer or choke must be so chosen, that 
positive potentials within the indicated limits are pro- 
vided on the transistor base when the contacts are open. 
Base emitter potentials close to those required can be 


obtained by including in all these schemes an additional 
germanium or silicon diode in the emitter-base circuit 
in a manner similar to Fig. 5d, 


4, Some Results of Experimental Investi- 
gations 

To illustrate the processes taking place in some of 
the schemes for voltage regulation described in this 
paper and to confirm some of the mathematical rela- 
tionships given, a number of voltage and current oscillo- 
grams are given below in Figs. 7 and 8; these oscillo- 
grams were taken on a de generator with a power output 
of 1500 w at a rated voltage of 28 v, controlled by a 
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Fig. 7. Oscillograms of voltage 
of the circuit of Fig. 5a. 


regulator consisting of a vibrator in combination with a 
transistor, The minimum speed of revolutions of the 
generator armature at the rated voltage and at full load 
was 1550 rpm with the generator still cold. A transistor 
of the P4A type was connected in series with the excita- 
tion winding which had a resistance of 8 ohms; a VG-10 
valve was in parallel with the excitation winding; the 
load on the generator consisted of a storage battery 
which was being charged at a current of 15 amp, To in- 
crease the speed of the regulator operation an accelera- 
tion resistance Ry ¥ 0,14 Rp was connected in the circuit. 

Figures 7a and 7b represent oscillograms of the gen- 
erator voltage uc, of the emitter-collector voltage ug,, 
of the emitter-base voltage u,}, and of the collector 
current i, of the transistor, The potential difference 
Ugh when the transistor is passing current is shown above 
the zero line, the blocking u,} potential—below the 
zero line, 

These oscillograms illustrate the operation of the 
scheme shown in Fig, 5a at n=1500 rpm. The average 
value of the generator excitation current was 1,43 amp, 
the generator voltage was 28.2 v. To obtain the block- 
ing voltage u,}, an interstage transformer from a 
Rodina-47 receiving set was used, The scales of the up 
and Ue- voltages on the oscillograms were the same. 
The u,,, voltages are shown on a magnified scale, 

The emitter-collector voltage ug, is here approxi- 
mately equal to the generator voltage u,,; thus relation- 
ships (9) and (10) are confirmed. The shape of the Usb 
voltage variations is approximately the same as that of 
the Ug, variations; the slight differences observed be- 
tween these two are due to deviations of the vibrator 
natural frequencies. 

The other schemes work in a similar manner, 

Figure 8 shows oscillograms of current ip, and also 
of voltages uc and Ue, illustrating the blocking of the 
transistor by means of a choke coil according to the 
scheme of Fig.5d, The primary winding of an inter- 
stage transformer from a Rodina-47 receiver-set was 
used as a choke coil. The emitter-base part of the tran- 
sistor was shunted by a semiconductor diode of type 
DG- Ts24 in order to limit the blocking signal. 


CONCLUSIONS 
At present the industry of our country has mastered 
the production of germanium transistors for currents up 
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Fig. 8. Oscillograms of u, 
illustrating the operation of the 
circuit of Fig, 5d. 


to25 amp [6]. By means of the above described schemes 
it is possible when regulating the voltage of generators 
with excitation currents to 25 ampto have at the vibra- 


tor contacts being broken currents of only 0.3 to 0.8amp. 


The great majority of automobile generators have 
excitation currents of 1.5-3,0amp, When vibratory regu- 
lators in combination with transistors are used on these 
generators, the current broken by the vibrator contacts 
will be from 0.03 to 0.09 amp. 

By combined use of vibratory electromagnetic regu- 
lators and of transistors it will be possible to avoid ero- 
sion of contacts, 

The schemes reviewed here by no means exhaust 
the many feasible combinations of transistors with vibra- 





tory electromagnetic relays and regulators. However, 
by the use of these schemes it is already now possible to 
lower the currents being interrupted at the contacts to a 
considerable extent, and thus to considerably increase 
the dependability of vibratory electromagnetic voltage 
regulators. 
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A method is presented for computing the parameters of circuits for stabilizing the speed of compound de and 
ac drives using a special three-winding transformer, A graphical method for plotting the circuit characteristics 
of the drive stabilizer is given; the non-linearity of the magnetization curve of the three-winding transformer 


core is taken into account, 


The methods presented are illustrated by an example in which the parameters and characteristics of a 
stabilization circuit are given; the circuit controls a dc drive powered from the ac line through a rectifier. 


Circuits using a special three-winding transformer 
for stabilizing the speed of dc and ac motors have 
found contemporary use [1]. Such circuits have a num- 
ber of important advantages, In particular, they permit 
one to realize feedback compensating for perturbations 
with a single device, to decrease the weight of active 
materials in the controlling device, to simplify the 
adjustment of the system, and to decrease the variety 
and number of circuit elements, The advantages of 
such circuits permit them to be widely used in a num- 
ber of dc and ac electric drive systems. The correct 
choice of the stabilizing circuit parameters is impor- 
tant in designing such systems, At present, problems re- 
lated to the computation of such stabilizing systems 
have not been dealt with in periodical literature, This 
article presents a method for computing the parameters 
of circuits for stabilizing the speed of dc motors with a 
choke regulator connected to the ac line through a rec- 
tifier, and three-phase motors with a Shenfer rotor; a 
method for calculating the characteristics of the three- 
winding transformer used in the speed stabilizer of the 
drive is also presented, 


1. Designing a Three-winding Transfor- 
mer for a DC Motor Speed Stabilizer 

A three-winding transformer is used as the com- 
pounding element in the circuit for controlling and 
stabilizing the speed of a dc motor shown in Fig. 1. It 
can be seen from this circuit that the load current I~ 
flows through the current winding of the transformer; 
the voltage winding in series with an additional resis- 
tance Rq is connected to the supply voltage U~, and 
the secondary winding feeds the stabilizing winding of 
the controlling device through a rectifier, It was shown 








in [1] that the expression for the current in the second- 
ary circuit of such a transformer, I», is 

I, =k(U. —kl-), (1) 
where 


k ah WwW, WR, 


Er) = Wr 
Wy y/ + nant (Ge) s 


if the voltage and current windings are connected op- 
posing. 

Here W; is the number of turns in the current win- 
ding of the transformer; W,, the number of turns in the 
voltage winding; W, the number of turns in the second- 
ary winding; Rj the resistance of the voltage winding 
referred to the secondary winding; x’, the inductive 
reactance of the magnetizing circuit referred to the 
secondary; and R,, the resistance of the secondary cir- 
cuit Ra=r,+R,. 
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The expression for the rotary shaft speed of the 
motor connected to the supply through a rectifier is 
U_—I1_Re (2) 


c 


n= - 
n 


where 
Re = bs (Rk, + Rat Rp), cn = | 
ky ky, 

Here R, is the rectif’ resistance, R, and R, are 
the resistances of the motor armature and brushes, re- 
spectively, kj and k,, are rectifier coefficients for cur- 
rent and voltage, > is the motor field winding flux, 
and c is a constant, 

Comparing equations (1) and (2), we can conclude 
that the three- winding transformer is a static speed 
transducer when the condition 


. (2a) 
W,, 


- Ke = const 
is met, 

The initial data for designing the compounding 
transformer are the parameters of the rectifier, motor, 
and transformer load, It will be assumed that the cur- 
rents and voltages in the transformer are nearly sinu- 
soidal in all operating modes of the stabilizing system, 
and that the rectifier resistance and coefficient are 
constant, Such assumptions are fully justified, since 
the transformer operates over a linear portion of the 
core magnetization curve. As was noted in [1], k,, va- 
ries over a small range. Therefore, by setting k,, at 
some value within this range, it can be considered 
constant with sufficient accuracy, It will be further 
assumed that the load resistance R, is significantly 
greater than the secondary winding resistance, and that 
the effect of the magnetizing circuit is negligibly 
smalli(X,,,* ©). The assumptions mentioned are justi- 
fied by the good agreement between the computed and 
experimentally found transformer parameters, It is 
recommended that a toroid core made of a high-per- 
meability material (type NP-50, NP-65, etc.) be used 
for the flux path of the transformers in order to dimin- 
ish the leakage flux and eddy current losses in the 
steel. 

We will keep the operating range of the transfor- 
mer within the limits of the linear portion of the mag- 
netization curve so that the signal will be undistorted. 
Considering the assumptions made, equation (1) can be 
written in the following form: 





I we KyE,, 
where Ra + KR: ’ 
En = Can, (4) 
_ Ws (5) 

_— Ww, ° 

Considering (3), the power in the load, P, , can be 

expressed in the following way: 
p K,°E,,"Rs (6) 

L 





(a +K,2 Ry 
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Let us find the condition of maximum power trans- 
fer to the load by the transformer. To do this, we will 
find the greatest maximum of the function P=f(Ky. Ry): 





Sok = (Ry + Ku?Ry)* — 2K (Ry + Ku?) = 0, 
2 
whence 
R, _ R,K,,’, (7) 
and oP, = 
a(K,") 


= (Ry + Ky2R;)? — 2Ky? (Ra + Ky®R,) R, = 0, 


whence 2 As 
a (8) 
It can thus be seen from equations (6) -(8) that the 
greatest maximum condition for function P, = f( « Re) 
is expressed by the equation 
E 2 
n 
R, = Pi, . (9) 
It can be shown, by using the initial set of trans- 
former equations given in [1], that the expression for 
the current in the voltage winding is 


U_’+1_'R: 

Ry’ + Ry . (10) 
in which the apostrophes denote that the corresponding 
quantities are referred to the secondary winding. 

Solving equations (5) and (10) together in terms of 

I, Ux ., we get a final equation for the current in the 
voltage winding: U_+1_Re 
a, a 

The voltage across the voltage winding of the 

transformer can be expressed by the following equation: 

U,, = 4.44/W B mabkegs 10°, (12) 


where B_,, is the effective flux density in the transfor- 
mer core; a and b are the dimensions of the sides of the 
core cross section; and kgf is the core fill- factor. 
Taking (8) into account, the solution of (12) for 
the product of the sides of the flux path cross section 


assumes the form 
V BR: 108 


LAAB, keg Wi, 





I,! = - 


e*= 


(1) 


ab = 





(13) 


Keeping equalities (3) and (8) in mind, the expres- 
sion for the ampere-turns of the secondary winding can 
be written in the following manner: 

E,W, 
IWa= $y. 

If q, 4, and kg are the wire cross section, current 
density, and fill- factor of the flux path window for the 
various windings, and Dg, is the diametez of this window 
the joint solution of equation (14) and the window fill 
equation: 


(14) 





n(Dagy— a)? “> 





15 
—T_ = (15) 
_ gWs0-* = q,W,,10-* rs qW,10-* xD? 
ke key kf 4 














gives the expression for the average diameter of the 
toroidal flux path: 








Day =a+ 
TW IVT ATW AO*  2E,WAO : 
/ aS ~ RARE, nA ke, Arkenks, + Def. 
(16) 


Hence the average line length is 


lay =x (a -t- 
a1,W i0-* 
+/ mA ke; 
(17) 


The amount of additional resistance in the voltage 
winding circuit can be found from the equation 





4], W, 10-2 
TAL KE, 





2E,,W A0-* 
2 


Ww 
Ry = Ry — 2 407, (18) 


u 
where 
( rey 9 tu q,W ,,10-* - ol a 
a . 
~lavhiy 

Ay is the thickness of the voltage winding cross section 
and p is the resistivity of the winding wire. 

The expression for the resistance of the secondary 
winding is ply Wa 


3 > q2 10°°, (19) 


Ly = 2a + 2b + 4hy = 





where 


ly, = 2a + 2b + 8Ay + 40g = 


‘= u . qWs ) 
= 10° + 2,——- 10° 
2(a+b+ 4—GA- 10% + 2 ee 
(A, is the thickness of the secondary winding cross sec- 
tion). 
The weight of the copper in the transformer is 
Yu fier . JytRy Je?r 
G= Bart ast] eo 
(here y ,, is the specific gravity a the copper). 
The total weight of the transformer is determined 
by the sum of the weights of copper and steel: 


Gror = Gu + Vee Yors (21) 


where yg is the specific gravity of the core material, 
and V,, is the volume of steel in the core: 


Ver = ablavirs- (22) 


This method for designing a transformer is illus- 
trated by an example given at the end of this article. 


2. Designing a Three-winding Transfor- 
mer for Stabilizing the Speed of an AC 
Motor 

Three- winding transformers are used as compoun- 
ding elements in circuits for stabilizing the speed of a 
three-phase motor with a Shenfer rotor, and a two- 
phase motor with a thin-walled rotor, which are shown 
in Figs. 2a and 2b. These transformers are built like 





































Fig. 2. CW - motor control winding, LW ~ motor line 
winding, Ug, - feedback voltage « 


the transformer shown in the circuit of Fig. 1, except 
that an additional inductive reactance X, is connected 
in the circuit of the voltage windings; X; is a choke 
with an air gap, 

It has been shown in [1] that if the current and 
voltage windings are connected opposing, and X; = 
(Ri sing)/cosg, the expression for the current in the 
secondary circuit is 

T, = kgU. — I ky (23) 

Using the designations mentioned above, we get 
the oS expressions for the coefficients k, and kg 
when X_, 





~ V (RK sR 


RW, 
W008 @ 
hs = VK + RP Fa hy 

Here ¢ is the phase angle between the current in 
winding W, and the voltage across Wy. 

For the circuit shown in Fig. 2a 9=¢,~30° (see 
Fig. 3), and for the circuit shown in Fig. 2b 9, = ¢, 
where ¢, is the phase angle between the load current 
and the voltage. 1_RwWw 

We will designate the expression // = ay 

or u 


by the letter A and write equation (23) in the following 
form : 





KA 
fa an <p eel a 
" (AKL? + Ae)? + APR, ? (24 . 











Considering (24), the expression for the power in 
the transformer load can be written in the following 
form: 
K,," A*R; 

= 2 = 4 
PL = Isha = TERS RS ROR Me * (25) 





Let us find the condition of maximum power tran- 
sfer to the load by the transformer.To do this, we will 
find the greatest maximum of functionP, =f(K,, Ry): 





oP 
ane = RitKu' (1 + te*@) — Ry* = 0, 
whence R, = RK, =a ' (25a) 
and = ap 
Ky == R,* — R,*Ky* (1 + tg*p) = 0, 








4 
whence R, = R,K,,? 4 (25b) 





It can thus be seen from equations (25a) and (25b) 
that the greatest maximum of P, = f(KE,, R,)is expressed 


by the equation degpeepmcmate 
Ky = Y= cos 
we (26) 


The joint solution of equations (25) and (26) 
gives the following expression for the resistance of the 
voltage winding: 


R, = 





A* cos @ (27) 
2P, (1+ cos@) ° 


Using the initial set of transformer equations given 
in [1), it can be shown that the expression for the cur- 
rent in the voltage winding circuit is 





hs 


f Ww W, \2 
Ky \/ UntKel + 2U.qy2lRe op 0089 +12 Ret (+) 








The voltage across the voltage winding can be 
found from the following dependence: 
Uy = 4.44 SW. Brmabks 10°°. (29) 
Considering (26), the solution of (29) for the pro- 
duct of the sides of the flux path cross section is 
| ab = G44jB, kip W, 
If Nyg= 0, then I, must also be 0, Consequently, 
equation (24) can be written thus: 








(30) 


U_=cl., (31) 
where pa RW, 
— Wycos@ * (32) 











cos ¢, = 0.5; Py = 1w; Ry = 1000 ohm; I, = 0.033 amp. 


Vira + AiK,,)* + Ai? tg*pXK,! ; (28) 





Hence, the transformer design consists of the fol- 
lowing series of steps: 

1, Using equations (27) and (31), we find c and Ry 
under starting conditions, The quantity A is determined 
from conditions when the motor is running steadily at 
its nominal rpm, 

2, We find K,, by using equation ( 26). 

3. Having been given a series of values for W;, we 
find the corresponding values of W, from (32), and 
knowing K,,, we then find Wp». 

4, Using (24) and (28), we find the currents in the 
secondary and voltage winding circuits, 

5. Using (30), we find the transformer flux path 
cross section for the values of W,, obtained in step 3. 

6. Choosing a flux path with a square cross section 
and expressing the wire cross section in terms of a per- 
missible current and current density, we find the dia- 
meter and average line of the flux path from (15). 

7, Using (18), we find the value of the additional 
resistance, 

8, We find the copper weight and total weight of 
the transformer from (20) and (21) for the various val- 
ues of Wj assumed in step 3, 

9, The last step is to construct a graph of the total 
transformer weight as a function of the number of turns 
in the current winding. We then find the number of 
turns W; which gives the minimum total weight, and 
we take all the parameters computed above which 
correspond to this W;. 

The graph Grote = f( Wj ) is shown in Fig. 4 for a 
transformer operating as in Fig. 2b, with the following 
initial data: I~ = 0.25 amp; UW = 50 v; U,, = 25 vy 














The transformer core material is 50NP. 
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3. Computing the Transformer Charac- 





teristics 


The change in the permeability of the iron in the 
core must be accounted for if the accuracy of the cal- 
culations of the three- winding transformer characteris- 
tics is to be increased. Using the equivalent circuit of 
a transformer operating in a dc motor speed stabilizing 
circuit (Fig, 5), we find the load current by the equi- 

















&s|}-— \ : ‘ valent generator method, The linear inductance in the 
g. . | voltage winding circuit must be taken into account tn 
y the calculation. 
oe If it is assumed that / _’ — I.’ (cos@_— / sin 1); 
mg the expression 
U_’ +1_'Rs (cos @, — 7 sin 
| l= ee ory tae 
0 50 100 180 W.., Turns r+ My 
Fig. 4. 
Ry ft G: 
—— —_ 
AS 8 
Ux, | An k 
| [Pies ic 
(1 
a 
Fig. 5. 
is justified for the current 1, (see the circuit of Fig.5b), f Ff 
whence hes = Roo . 
Uan=U.'—1,R,’ he (34) AB t+ IAM (37) 
If the transformer parameters are known, and a 
; ; nee number of values I~ are given, the voltages correspon- 
wi [U’ + (1 cos oy — Jl’ sin pi) Ra} Ri ding to them, Uap, can be found, We will solve equa- 
fi’ + Rs tion (37) graphically. Both axes will be given identical 
When cos ¢ = 1, equation (34)takes the form voltage scales (Fig. 7), We will plot the voltage Ey~ 
: . the voltage across the non-linear reactance X,,~along 
U., 2 UL UL +4’ Rs . (3s the ordinate, and the voltage across Rap will be plot- 
Ry’ + Rs ted along the abscissa. We draw a number of circles 
The resistance of the equivalent circuit with with radius U AB’ centered at the origin. Each such 
respect to terminals AB is equal to 
. 
RR, 
= 3 
Rap Ri’ + Ry 3% 


Fig. 6 shows an equivalent diagram of the trans- 
former, consisting of a source Uap, linear resistance 
Rap: and non-linear reactance X,, Connected in series. 
The current in this circuit is found from the formula 


‘i 


Fig. 6, 




















circle corresponds to definite values of the voltage U~ 
and current LJ. If the magnetization curve of the steel 
core is also plotted on these coordinates, using the 
equations 

Ey’ = 4.44 [W,BwSke¢ 10°, (38) 


Hylay= V 2-0.42Ites'W 2, (39) 

the points where this curve intersects the circles will 

be operating points for the transformer, Having a family 
of curves E', = f(UresR ap) for U_, = const , and 


E'y = f(IpesRap) for LL = const , we find the required 
dependencies 1, = f(U_,) and I,=f(I~). The function 
I, = f(U_) for 1. = const is analogous to the function 
n= f(U_) for a drive with a constant load torque on the 
motor shaft. The function I, = f(I.) for U. = const 

is analogous to the mechanical characteristic of the 
drive, n= f(M), for a constant supply voltage. 

Example 1, Design of a stabilizing three- winding 
transformer for a “supply-rectifier — dc motor™ system. 
The drive speed as a function of load current is expressed 
by the equation 

—Rel U_—viAdiL 
re = ~ 37-105 


on 





(ex. 1) 


The nominal load current L~ = 20 amp;nominal shaft 
speed n= 3000 rpm; transformer power rating P; =2 watts 
(U,= 80 v, I,= 0.025 amp). Core material is 50 NP, 


Considering the size of the load current (I.=20 amp), 


it is best, for technical reasons, to make the transformer 
with a one-turn current winding(W ; = 1), The area oc- 
cupied by the current winding must be chosen on the 
basis of an appropriate value of Dp. 

1, The supply voltage is found from expression (2), 
using the nominal values of I~ and n: 


VU = cn +1 Re= 


== 3.7-10-*.3000 + 20-0,15 = 14.1 v. 


(ex,2) 
2. The circuit resistance of the voltage winding is 
found from (9): 


(c,")® —_(3,7-10-8-3000)? 





Ry = iA, = a = 15.2 ohms. (ex, 3) 
3, The number of turns in the voltage winding, ac- 
cording to equation (2a), is 
Ww, ~i = te = 101 tums. (ex. 4) 
4, The transformation ratio K,, equals 
Ky= V/ a= Va =14.5  (€X- 5) 


according to (8), 
5, We compute the number of turns in the secon- 
dary winding from (5): 


Wa= KW, = 14.5°101 = 1470 turns, = (€X. 6) 








6, The current in the voltage winding circuit is 
found from (11): 


U_+7_Re 14.144 20-045 
2 3-452 








= 0.562 amp. (ex. 7) 


7, Choosing the current density in the secondary 
and voltage windings to be A= 3 amp/mm’, we find the 
wire cross-sectional area of these windings: 


I, 0.562 2 


= 0.1842 mm”. (ex, 8) 


Tae es 
8. We choose B,,, = 12000 gauss from the magneti- 
zation curve of NP= 50 steel, and we make a=b in or- 
der to optimize the core geometry, 
9. In correspondence with (30), the side of the to- 
roid flux path cross section is 


Ya. VAR 
ein’ T44jB, keg W. 


/ 











=~ aaa 108 ¥ 2.15.2 
= 4.44-400-12000-0.43- 104 





= 0,77 cm. (ex, 9) 


10, According to (32), the average diameter of the 
flux path, for D)= 10 mm, is 


Day = 4+ 





41,W, 10°? 2E&,W,10-* 





+ Do? = 9.77 + 








mA kfy ~ *ReAskfe 
4-0,562-101-10-? 2-3.7-10-$- 3000-10? , 4 = 
3.14-0.3-3 3.14-0.15-0.3-3 

= 2.3 cm. (ex, 10) 

11, The length of the average toroid line, 7 ay’ is 

found from (33): 
lay= wD. .= 3.14-2.3 = 7.24 cm. (ex, 11) 
12, The volume of the flux path is 

Vet = @hykge = (0-75)*-7.24-0.43 = 1.85 cm®. (ex, 12) 


13, The average length of a turn in the voltage 
winding is found from the equation 


9,,W,,10-* 
‘wu = 4 («+ kfylay | 








: 0.1842-101-10-2 
=4 (0.77 +357 ) =3.34cm, (eX 13) 
14, The additional resistance is 
p 
Rg = Ri— A 10-1215.2 — 
qu 
_ 3.34401 
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3) 


4) 


],, ma 






































Fig. 8, 
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16 
I=20amp 
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Fig. 9, 


15, The average length of the secondary winding 
is found from the equation 
WuW 10% g,W210-* 4 
=4 2 <a 
Va (: bn lav key , laykf, } 


0.1842-101-40-* 0,00814-4470-40-2 
7.24-0.8 + 7.24-0.3 





— (0,77 42 





) = 3,98 cm. 


(ex, 15) 
16, The secondary winding resistance equals 


Plas 2 3-98-1470 
q2 = 57-0,0081 10°= 128 ohms. (ex, 16) 





r= 





17, We use the well-known formula 


Tw TOR, 12. 
oe | oad ba 


(0.562)?-0.318 — (0,025)?.428 
= 8.9.57 [ces an) = 101 g 


(ex, 17) 


to find the weight of the copper, 


18, The total transformer weight is equal to 


Gy = Gy + tstVst = 10,4 + 8.2-1.85 = 25.2 g. (ex. 18) 


Example 2, Graphical computation of the charac- 
teristic U. = f(1,) and I,=f(I~) for the three- winding 


transformer whose parameters were found in Example 1, 


The circuit parameters are R, = 15,2 ohms; R,= 3200 
ohms; W; = 1; Wy= 101; W, = 1470;r, = 128 ohms. 

1, Let us determine the equivalent values of resis- 
tances, currents, and voltages: 


Ww 1470\2 
Ry'= Ry (iv : i= 15 15.2( or) = 3240 ohm, 


1’= 1. W, = Ti’ 
W, 1470 
0! = Uw =U. oH 


(ex, 19) 

2. Having been given a series of values for L~, we 
will find the equivalent currents, The results are shown 
in Table 1, 


Table 1 





amy 0 | 


Pry ma| 0 1.36] 2.72 | 4.07 | 6.8 10,2| 13.6 | 17 


4 | 6 10/15 | 20 | 25 




















3, Having been given a series of values for U., we 
find the equivalent voltages, The results are shown in 
Table 2, 























Table 2 
| U_.v 2 2 | 12 “ 16 
| ol. 29.1 87.2 | 116.4 | 174.5 | 232.5 
Table 3 

¢7 aS 4 | 10 | 45 25 





~ [eon] a | 48.2 | 43.4 | 





4, Using equation (36), we find Rap: 


Ri Ra 3240 - 3338 
AB= Ry +h™ 3240 + 3338 = 1625 ohm. (ex. 20) 





R 


5. Using equations (38) and (39), we plot the mag- 
netization curve of the 50NP core material on the co- 
ordinates E,,' and I’,..Rag . The curve obtained for the 
transformer core is shown in Fig. 7. 

6, We find the U,, corresponding to the different 
L." found in step 2, for each U~ assumed in step 3, 

The results of these computations are given in Table 
3 for U., =8 volts, 

Fig. 7 shows the construction for finding the trans- 
former operating points on the plane E,'= f(T’ es" ap) 
Using the expression 

I, = E,'/ Re, (ex. 21) 
where 
R,= RK a re (ex, 22) 


we find the dependence I, = f(L_ (cf. Fig. 8). The func- 
tion. I,= f(1~), observed experimentally with a trans- 
former made according to the design data of Example 
1,is shown also in Fig. 8 (circles denote the calculated 
curve, and crosses show the experimental one), We 
find the dependence U,, = f(I,) for L. = canst from anal- 
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ogous constructions for other values of U~, Fig. 9shows 
U~ =f(1,) forl.=20amp. This figure also shows the 
function U~ = f(1,) observed experimentally with the 
transformer previously mentioned. 

The good agreement between the calculated and 
experimentally observed characteristicsshows that the 
assumptions made above are justified in the design of 
a three- winding transformer. 


LITERA TURE CITED 
1, O, B. Rosenbauli and R. N. Rodin, "Motor speed 
stabilizing circuits using a three-winding trans- 
former,” Avtomat, i Telemekh, 20, No, 3( 1959)? 


* See English translation, 























AN APPROXIMATE DETERMINATION OF JET REACTION 
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A method is presented for computing the approximate magnitude of the force with which a jet of working fluid, 
issuing from a nozzle, affects the flapper. The results obtained are compared to experimental data. 


In nozzle-flapper hydraulic amplifiers [1] the jet 
of working fluid issuing from the nozzle acts as a force 
upon the flapper. The magnitude of this hydro-dynamic 
force (let us call it the jet reaction) is an important 
factor in choosing the control element for the amplifier, 
and in analyzing the joint operation of the control ele- 
ment and amplifier. 

In recent years, a number of experimental and 
theoretical works concerned with the study of the inter- 
action between the jet and flapper have been completed 
[2-4]. This paper presents a method for calculating 
approximately the jet reaction; it agrees well with the 
experimental results given in [4]. 

One can consider that the jet reaction is composed 
of three forces: 

A. Force N,, which varies with the mass flow rate 
of working fluid issuing from the nozzle. This component 
can be found from the expression 


N, =m (V. oe V.), (1) 


where V,, is the projection of the working fluid velocity 
at the point of exit from the nozzle along an axis per- 
pendicular to the plane of the flapper, Vj, is the pro- 
jection of the working fluid velocity at discharge [over 
the area of the nozzle end-face of external radius 
Thy = dy;/2 (see figure)] along an axis perpendicular to 
the plane of the flapper, and m is the mass flow per 
second issuing from the nozzle. 

Setting the velocity of the discharge equal to zero 
(Vj = 0), and assuming that 

Q nd.* 


m= Qp, Vo= “F? F = * (2) 


4 , 





where Q is the volume rate of working fluid discharged 
through the nozzle, p is the density of the working 
fluid, F is the sectional area of the nozzle end-face, 
and d, is the nozzle diameter, we get 
2 
N,= a . (3) 

B. The force generated by the pressure of the wor- 

king fluid on the nozzle end-face area; it equals 


Newt 
7 ee (4) 





where p,, is the working fluid pressure at the nozzle 
end -face (in atmospheres, gauge). 

The pressure p,, is assumed to be identical over 
the entire area of the nozzle end-face. 

C. The force exerted by the working fluid pressure 
in the gap between the end of the nozzle and the flap- 
per. This third component can be represented in the 
form 


N, =2n \ p,rdr, (5) 


where r is the flow radius of the nozzle end, r, is the 
nozzle radius, tr; is the external radius of the nozzle 
end and py is the working fluid pressure in gap between 
the nozzle end and the flapper, over a radius r. 
If we assume that the pressure p, over the end -face 
radius varies linearly from p,, to 0, that is, 
Pr = Pc (4 = =) ’ 


ee 


(6) 


then we get 





d.*—d2 d,3—d,? 
4 6d 


Ny = ape (—* c aa c 


an 
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h,mm |p. kg/cm? rm N,G Ne-G = | [N—Ne].G 
0.15 10 0.717 96 95 f 
0.20 10 0.667 i101 100 4 
0.25 10 0.670 109 115 6 
0.15 20 0.700 191 190 1 
0.20 20 0.647 200 220 20 
0.20 30 0.675 303 300 3 
0.25 30 0.638 321 320 | 
0.15 | 40 0.637 375 375 — 
0.20 40 0.664 400 400 if 
0.25 | 40 0.644 429 420 9 





The jet reaction, being the sum of the three com- 
ponents mentioned above, will equal 


4pQ? n / @,° ds 
N= +4 (4- +“) Pe 


The formula for volume discharge rate through the 
nozzle is 


(8) 


nd? agi (9) 
Q = be V2 (Pi— Pc); 
where p, is the pressure in the receiver, and y ,, is the 
inherent discharge rate factor of the nozzle (without 
the flapper); we find from this that 


Mee nea 8pQ? (10) 
Pe = Pi Wd a? 
In this case, 
40% 11 





xn [4° q,* 8pQ? 
+5( 2 + dy )(m— wd.* p.? ): 
If the discharge is into the atmosphere, Q can be 
written in the form 


Q=padh VW = p,, (12) 


where yp is the discharge rate factor of the nozzle and 
flapper combined, and h is the gap between the nozzle 
and the flapper. 

Expressions (11) and (12) enable one to find the jet 
reaction N, if the parameters dc, dj;, .. wp» h and p; 
are known, The supply pressure p, can be substituted for 
pressure p,, since they are statically related. 

The calculated values N of the jet reaction were 
compared with experimental data N, for a nozzle of 
the first type with diameters dc = 1 mm, and d),=1.35 
mm [4], working in AMG~=10 oil at a temperature of 
50- 55°C. Expressions (11) and (12) give 

13 
N = [8.03 (uh)? -|- 0.00865) p, ies 
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for a nozzle with these dimensions, and an experimen- 
tally determined p= 0.9. 

The values of discharge rate factor » for the noz- 
zle and flapper combination which were used for com- 
puting the jet reaction from formula (13) were taken 
from reference [5]. The calculated results and experi- 
mental data are given in the table. 

A comparison of the calculated and experimental 
values of jet reaction given in the table shows satis- 
factory agreement, Hence the method being considered 
can be recommended as an approximate determination 
of the force on a flapper, if the amplifier parameters 
approximate those given in [4]. 


CONCLUSIONS 
1, Computational formulas are given which per- 
mit one to calculate the approximate force on a flapper 
as liquid flows around it, 
2. It is shown that calculated and experimentally 
found values of the force on a flapper agree. 
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rhe paper studies the effect of the dynamics of a deterministic function on the transfer coefficient for fluctuations 
when regular and random functions pass simultaneously through a nonlinearity. 


1, Statement of the Problem 

In analyzing the operational precision of both lin- 
ear and nonlinear dynamic systems, many engineering 
computations treat only the first probability moments 
of random functions, The approximate statistical line- 
arization theory [1, 2] which is based on these engin- 
eering requirements has made it possible to perform 
fruitful investigations of nonlinear systems; these inves- 
tigations have exposed very interesting system proper- 
ties which are both harmful and useful from the point 
of view of engineering practice, According to statisti- 
cal linearization theory, two functions are considered 
statistically equivalent if they have identical first and 
second moments when a specified distribution law governs 
their arguments. Here, the first or second statistical 
linearization method is used to determine the transfer 
coefficients of a nonlinear element; Kg is the transfer 
coefficient for the useful signal, and K, is the transfer 
coefficient for fluctuations. 

In many engineering problems, the value of the 

mathematical expectation for the random perturbation 
varies insignificantly, whereas the useful signal varies 
relatively rapidly, Therefore, it is natural to consider 
the dynamics of the useful signal in computing.the 
coefficients Ky and K,. In this regard it is of interest 
to study [3] which provides a method for taking into 
account the dynamics of the useful signal, Using this 
method, we derive analytical expressions which make 
it possible to determine how the transfer coefficient 
for the fluctuations is affected both by the parameters 
of the distribution law for the random process and by 
certain parameters which characterize the dynamics 
of the regular component, 
2. Determining the Transfer Coefficient 
of a Nonlinear Element for Fluctuations 
while taking into account the Dynamics 
of the Regular Signal 

















We shall assume that the nonlinearity described by 
the functional relationship v= f(v9),where vo is the in- 





put and y is the output, is subjected to three compon- 
ents vo= 8 + x + y, where B is the constant component 
which contains the practically constant mathematical 
expectation of the random perturbation; x is the ran- 
dom component which has a zero mathematical expec- 
tation, and y is any regular signal, If K is the equival- 
ent transfer coefficient of the nonlinear element for 
fluctuations, and k is the transfer coefficient for the 
useful signal, then the instantaneous value of the dif- 
ference 


rg=f PB+2+y)—(Ke+hky) () 


is caused by the neglect of the higher harmonics and a 
certain constant component A at the output of the non- 
linearity, The instantaneous value of the difference x, 
is a function of x and y whose properties are charac- 
terized by the distribution laws q(x) and p(y) respec- 
tively. The criterion which we shall choose for selec- 
ting the optimal value of the equivalent transfer coef- 
ficient will be the minimum mean-square error for the 
instantaneous difference x,. Then the problem of find- 
ing the equivalent transfer coefficient reduces to find- 
ing the minimum for the expression 


M= { { Ge (2) p Waedy. (2) 


With respect to the distribution function p(y) we 
must note the following, Assume that the useful signal 
varies according to a sinusoidal law y=asin€, and that 
each of the possible values of y in the limits -a to +a 
corresponds to an innumerable set of values € : 


= mh , (1 y «= 
g. 5 +——oresin = (k = 0, +1, +2...) (9) 
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From this we treat € in accoraance with [4] as a certain random quantity; we find 


| : 





Py) = [* Vey 


Assume now that the quantity € is uniformly dis- 
tributed over the interval - = € = +1; then in the in- 
finite sum (4) only two terms corresponding to k=0 
and k= 1 will be nonzero, Then p(y) can be written in 
simpler form: 





( — » lyl<a, 
| xa V 1—(4) (5) 


ply) = a 
0, ly| >a. 


In this case the distribution function p(y) will 
characterize the dynamic properties of the regular sig- 
nal y which is studied jointly with the random process 
x that is characterized by its distribution function; thus 
p(y) will, so to speak, itself become random, since the 
phase of the regular signal is actually distributed uni- 
formly over the interval - s = € = +m when this sig- 
nal is treated jointly with the random signal. 

Taking the above into account and performing the 
minimization of expression (2) with respect to K, we 
obtain the following result after simple computations 
if we assume that the random signal is distributed nor- 
mally: 


k= | | GB +2+v) a(z) ply) dzdy. (6) 
—oo —00 


where X* is the mean-square value of the input signal x. 
For the purposes of further computation we shall 
modify the expression for K somewhat, Assume that a 
constant component y acts at the input of the nonlin- 
ear element instead of the random signal x; then the 
constant component of the output signal is written as 


sw= (+B +uP Wav. 
- (1) 


Now expression (7) shall be treated as a certain 
fictitious nonlinearity whose input is subjected to the 
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k 
o (6 +o arc sin “|, ly| <a, 
w a (4) 


ly|>a. 


random signal x; then, 
co 
Ky, = “ \ xg (x) q (x) dz. (8) 
—oo 
Substituting Eq. (8) into Eq. (7), we obtain an ex- 
pression for the equivalent transfer coefficient which 
coincides with (6), Such a treatment permits us not only 
to simplify the analytical computations but also exper- 
imentally to determine the fictitious nonlinearities by 
studying the effect of the regular signal on the original 
nonlinearity, It is convenient to perform the analytical 
computation of the fictitious nonlinearities and the 
equivalent transfer coefficients K, by representing the 
nonlinearities in terrns of “choppers” [5]; the expression 
for the nonlinearities are shown in Table 1, The results 
obtained by computing the fictitious nonlinearities and 
the values of K, are cited in Table 2, Figure 1 shows 
the normalized "wansfer coefficients K, and K, , for 
fluctuations in the cases of a relay element ant” an el- 
ement with “negative overlaps” (cf, third characteris- 
tic from top in Table 1) when a constant component, 
noise and a sinusoidal signal act on these elements, 
The graph shows the strong éffect of the normalized 
regular signal proper and the constant component nor- 
malized to the amplitude of the regular signal; the 
constant component shifts the value of the transfer co- 
efficient into the negative region for 1 = 1 and €* 1,2, 
For a relay element with a dead band that is sub- 
jected to the same three coriponents the transfer coef- 
ficient K, also passes into the negative region (Fig. 2). 
Figure 3 shows the values of the transfer coefficient 
K,, for fluctuations in the case of a linear element 
wit a dead band; Fig. 4 shows the transfer coefficient 
for fluctuations in the case of a limiter that is 
pi hat to a constant component, a sinusoidal signal, 
and noise with a Gaussian distribution, From the figures 
the strong effect of the dynamics of the regular compo- 
nent on the equivalent transfer coefficient is evident, 
3. On the Possibility of Generalizing 





the Effect of the Dynamics of an Arbi- 





trary Signal on the Transfer Coefficient 





for Fluctuations 





In closed-loop control systems the form of the pro- 
cesses can be described with sufficient accuracy by the 
expression y= Ce #sin Qx, In order to estimate the 
effect of a signa! having such a form we expand this 
expression into a Fourier series: 


























Table 1, Analytical Expression for Inertialess Nonlinearities 





Shape of the characteristic] Mathematical expression in terms of “choppers” 
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The expression for the coefficients contains still sum of harmonic waves with random phases. For our 
another parameter T = 2% w which represents the period case we introduce the substitution a,,,),=b, and obtain 


over which the regular signal is analyzed, 


the expression for the distribution function of a regular 


In [4] the distribution function is computed for the signal with a zero average value: 
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Then the transfer coefficient for fluctuations can be written as 


m| _ 
Pe 


K,= 


—oOO 


co oo 
\ | (e+ $+ B + ¥)P (2) Won (W) ded, 





Thus when necessary (although it involves great 
difficulties in computation) it is possible to determine 


the effect produced by the dynamics of such a rather 3. 
complex signal on the behavior of the transfer coeffi- 
cients for fluctuations, 

In conclusion it is necessary to note the substantial 
effect produced by the constant signal and the dynam- 
ics of the regular signal on the behavior of the transfer 4, 
coefficient for fluctuations; this underscores the neces- 
sity of exercising great care when using statistical lin- 
earization methods, 5, 
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The paper studies the substitution of an equivalent servosystem for a servosystem with combined control; the 
equivalent system has control only with respect to deviation, and its transfer function has a form convenient 
for application of the logarithm frequency response method. Nomograms are cited which facilitate the plotting 
of logarithmic frequency responses for the elements in the equivalent servosystem. 





A substantial improvement in the dynamic proper- 
ties of automatic control systems is possible when com- 
bined control methods are used; in such methods con- 
trol relative to the basic perturbations acting on the 
system is combined with control relative to deviation 
[1-5]. 

In [6] the possibility of computing servosystems 
with combined control by the logarithmic frequency 
response method was demonstrated, This paper devel- 
oped a convenient method for plotting the logarithmic 
frequency responses for a servosystem with combined 
control while taking the perturbation signal into ac- 
count only in the low-frequency region, 

In this paper we develop a method for plotting the 
logarithmic frequency responses of servosystems with 
combined control over any frequency range, This me- 
thod is based, as in [6], on replacing the servosystem 
with combined control by an equivalent servosystem 
with control solely relative to deviation. However, the 
transfer function of the equivalent servosystem is given 
a form which makes it easy to plot the logarithmic 
frequency responses of the system over any frequency 
range when nomograms are used, 


1, The Transfer Function for the Equiva- 
lent Servosystem 

Figure 1 shows one of the possible block diagrams 
for a servosystem with combined control, This system 
incorporates control relative to the signal g(t) in addi- 
tion to control relative to the deviation of the control- 
led quantity x(t) from the magnitude of the control 
signal g(t). 

To establish the relationship of the quantities ¢(t), 








y(t), and x(t) in the system to the signal g(t), the trans- 


forms of the differential equations can be written in the 
following form: 

& (p) = g (p) — x (p), 

y (p) = © (p) Wi (p) + ¢ (p) W, (p), 


z (p) = y (p) We (p). 


The transfer functions of the system for the error 
€ (p) and for the controlled quantity x(p) relative tothe 
signal g(t) are given by the following expressions on the 
basis of (1): 
1 — We (p) W, (p) 
1 + Wi (p) W2 (p) * 


Ws (p) [Ws (p) + W, (p)l 
1 +.Wi (p) We (p) 
Formulas (2) can be written as 
Ws (p) (Wi (p) + Wy (p)) 9-1 
1 — We (p) W, (p) 8 





®, (pr) = 
(2) 





®,., (p) = 





®.. (p) anes E + 


Ws (p) (Ws (p) + Wz (pr) 
;,4 (7) = ~T— Wp) W, wr) 
Ws (p) (Wi (p) + Wey (p)]) —* 
x [* T— Wy (p) W, @)- : 
From this it follows that a servosystem with com- 
bined control can be replaced by an equivalent servo- 
system with control solely relative to the deviation, 
The transfer function for the equivalent open-loop ser- 
vosystem, Weg, ‘Phis equal to 
Ws (p) (Ws (p) + We (p)I 
Veq)=—T— WwW, £4 
as is evident from (3), 
It is evident that the equivalent servosystem is 
formed from the same basic elements W,(p), W2(p); 
and We (P) which make up the system with combined 


% (3) 




















Wf) 


o}) x(t) 
Mh _@§ m(p) “(p) 




















Fig, 1 






































control. The element W g(t) in the combined system which 
performs control relative to the signal forms a direct 
positive coupling circuit with the element W,(p), and 
forms a positive feedback loop with the element W,(p). 
Figure 2 shows the block diagram obtained in accord- 
ance with (4) for a servosystem with combined control 
that is equivalent to the block diagram in Fig. 1. 

The transfer function of the equivalent servosystem 
(4) can be written differently: 


Ws (p) W, (p) Wi (p) 5) 
Weg\P) = 7 We! Wey |: 

This amounts to representing the equivalent servo- 
system in the form of a serial connection of the elem- 
ent with the transfer function W,(p)/ W_(p) that is con- 
nected by positive direct unity coupling, and the elem- 
ent with the transfer function W,(p) W_(p) that is con- 
nected by unity positive feedback, The block diagram 
of the equivalent servosystem (5) is shown in Fig, 3, 

We denote the transfer function of elements with 
positive direct unity coupling by F, gir (p) and the 
transfer function of elements with unity positive feed- 
back by Fiq,(p); i.e., we write 





Ww 
Pidiz.(P) = 1 + Woy 


Ws (p) W, (Pp) (6) 
ifo ©”) T— Wi) W, 0)" 





Then 
Veq'?) = Fydir(p) afb (P)- (7) 


From the above it is evident that the systems in 
Fig. 1 and Fig, 3 are dynamically equivalent. The sub- 
sequent analysis therefore consists of studying the dyna- 
mic properties of the system (7). For this purpose it is 
necessary to plot the frequency response of the system 
in Fig. 3 for open-loop conditions after first plotting 
the frequency response of the elements with transfer 
functions of the form 


— {p)....¢8) 


Fudirip) = 1+ Wie) Fife (P) = 7 WG): 


2. Plotting the Logarithmic Frequency 
Responses for an Element with the Trans- 
fer Function Fy,, = W'(p) + W'(p) 

We shall study the more general problem of plot- 


ting logarithmic frequency responses for an element 
with the transfer function 
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Fig, 2 
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Fai, (P) = W" (pr) + W" (p). (9) 
To do this the amplitude-phase characteristics for 
elements with transfer functions W*(p) and W''(p) are 
written as 
W' (ja) = H’ (@) e*”™, 


,” ° = HY’ ei?" ), 
W” (jo) (w) (10) 


Therefore 
W (jo) = H (w) &) = Hw) + HY (wo) 
where W(jw)=Fy;, (jw), or (11) 
ary el%o)—0'(@)} 4g 4 Fee cos {0” (@) — 9’ (@)) + 


4 12 
+ jy 2 sin (0° (w) — 6° (w)). — 

Equation (12) makes it possible to use the known 
[H"* (w) — H'(w)}db and [0°"(w) -9"(w)]* to plot the no- 
mograms for the geometric loci of the amplitudes 
(H- H") db and phases [9(w) - 6°(w)]°; these, in turn, 
are used to plot the logarithmic amplitude (H, db) and 
phase (9°) responses for the element with the transfer 
function (9), 

In [7] we plotted the nomograms for the geometric 
loci of the values for the amplitude A(w) and phase 
(Ww) characteristics for a closed-loop system: 

® (jo) = A (w) 2 (13) 
from the values of the amplitude H,)(w) and phase 
6 ,)(w) characteristics of the open-loop system (Fig. 4): 


Woy (i@) = Hey(w) e #° (14) 
The relationship between the amplitude- phase re- 
sponses of the closed-loop ®( jw) and open-loop 
WwW ol (jw) systems: 


Wo (/@) 
® Ge) = TF Wei Ga) as 
can be written as 
t o 
Aw or" =1+ 
(16) 


4 ? i 
+ Fg) cos )(@) — j Fete) 8p} (@) 


on the basis of (13) and (14), 
Comparing Eq. (12) with Eq..(16), we find that for 
the condition 


[H (w) — H’ (@)|db = — A (wpb, 
+ [6 (@) — 0’ (@)I° = #9 @)’, m8 
[H” (@) — H’ (w)ldb = — (Hopdb, 


+ [6” (@) — 9° (@)I° = F &)(o)° 
the nomograms plotted for determining the values of 
the amplitude and phase responses of the closed-loop 
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system from the values of the amplitude and phase re- 
sponses of the open-loop system can be used for deter- 
mining the amplitude difference (H** - H')db and the 

phase difference (9 — 6")°, 

Figure 5 shows the nomograms for the geometric 
loci used to determine the difference between the val- 
ues of the amplitude (H ~ H')db and phase (6 ~ 6")° 
characteristics for the elements with the amplitude- 
phase responses W(jw) and W‘(jw) from the difference 
between the values of the amplitude (H** — H')db and 
phase (8** — 6")° characteristics for the elements with 
the amplitude-phase responses W"(jw) and W"(jw); 
this is done on the basis of the nomogram in Fig. 4 and 
conditions (17), If, as is the case in the example under 
study, W'(jw) = 1 (i.e., H"(w)=0 db and 6*(w}=0°), 
then Eq. (17) is written as 


H (o) = —A (@), +8 (o) = ¥ 9 (@), 


(18) 
H” (@) = — Hoo), +8” (@) = F 4): 


i.e., the nomogram shown in Fig, 5 makes it possible 
immediately to determine the amplitude- frequency 
and phase- frequency responses for the system element 
that is connected by direct positive unity coupling. 


3. Plotting the Logarithmic Frequency 





Responses for an Element with the Trans- 





-20 08, 


fer Function Fy¢,(p) = W(p)/[1-W(p)) 





Assuming F yp, (jw) = B(w)e)8(w), 

W (jo) = H (jo) 2, (19) 
and the amplitude-phase response of the element with 
the transfer function F gfp(p) = W(p)/[1+ W(p)] is 

i 
B(@) 
i a> 
+ H to) © 9 (@) — fj Fy" 0 (w). 

Equation (20) makes it possible to plot the nomo- 
grams for the geometric loci of the amplitude (B, db) 
and phase (6°) from known H, db, and 98°; then it is 
possible to plot the logarithmic amplitude and phase 
responses for the element with the transfer function 
Figh (Pp) = W(p)/[1+ W(p)]. 

Comparing Eq. (16) and Eq. (20), we verify the 
fact that for the conditions 

B (@) = A (@), B° = (180 -} @)° 

H (w) = Hoyo), 8° = (180 + @)° 
the nomograms plotted from Eq. (16) and (20) are iden- 
tical, 

On the basis of the nomogram in Fig, 4 and condi- 
tions (21) Figure 6 shows the nomograms for the geome- 
tric loci of the values for the amplitude (B, db) and 
phase (8°) responses of the element with the amplitude- 
phase response Fig, (jw) as a function of the values for 
the amplitude (H, dh) and phase (9°) responses of the 
element with the amplitude-phase response W (jw), 





elo) = «= { + (20) 


(21) 


Conclusions 
1, In order to plot the logarithmic frequency re- 

sponses of a servosystem with combined control that is 
achieved by introducing coupling relative to the con- 
trol signal (Fig. 1), the block diagram of the equivalent 
servosystem with control only with respect to deviation 
must be represented as shown in Fig, 3, In that case the 
transfer function of the equivale1t open-loop servosys- 
tem is 





Veg?) =F aie €P) - ‘i 
i 
Fidic, (°) = 1 + Weg 
Ws (p) W, (p) 

Fu (”)= T— Wily) W, 7)" 
W,,(p) is the transfer function for control signal coup- 
ling; W,(p) is the transfer function for the servosystem 
elements encompassed by the control-signal coup- 
ling; W,(p) is the transfer function for servosystem ele- 
ments which are not encompassed by the control-signal 
coupling. 

2. In order to plot the logarithmic frequency res- 
ponses for the elements of an equivalent servosystem 
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with the transfer function F yqir(p) it is possible to use 
the nomograms given in Fig, 5; these make it possible 
to use the values of the amplitude and phase responses 
for elements with the transfer function W,(p)/ W,(p) 

to obtain the values of the amplitude and phase respon- 
ses for elements with the transfer function F, gir(P). 

3, In order to plot the logarithmic frequency res- 
ponses for elements of the equivalent servosystem with 
the transfer function F; q,(p) it is possible to use the no- 
mograms shown in Fig, 6; these make it possible to use 
the values of the amplitude and phase responses for el- 
ements with the transfer function W,(p)W.(p) to ob- 
tain the values of the amplitude and phase responses 
for elements with the transfer function F; q,(p). 

4, The nomograms in Figs. 5 and 6 can be used 
not only for plotting the logarithmic frequency respon- 
ses of servosystems with combined control, but also for 
plotting the logarithmic frequency responses of elem- 
ents forming a circuit with direct positive coupling 
(the nomograms in Fig, 5) and with positive feedback 
(the nomograms in Fig, 6), 
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The paper studies the use of the specified transfer function of the closed-loop system to synthesize the control 
program of a digital machine which operates in the closed-loop of an automatic control system. 


At present automatic control systems which include 
a digital machine (DM) are widely used in industry. 
Such systems have a number of advantages over con- 
ventional systems [1]; in particular, they make it pos- 
sible to design control systems with new properties, 
such as the selection of the optimum mode of operation, 
the ability to operate in self-adaptive fashion, etc, 
Moreover, systems including DM have a number of fea- 
tures which require special synthesis methods, 

By treating DM as a special kind of pulse element 
it is possible to use the Ya, Z, Tsypkin theory of pulse 
systems-[2] in order to compute systems that include 
DM, However, the problems of synthesizing a control 
program while taking into account the realization of 
that program in the DM have as yet been treated only 
cursorily in literature, ; 

Experience in applying DM to automatic control 
systems shows that in realizing the control program the 
program volume (in the sense of the number of arith- 
metic operations and “memory*® cells in the machine) 
is of great significance, This fact must certainly be ta- 
ken into account for synthesis, 

Another important problem in synthesizing the 
control program is the program stability, Notwithstan- 
ding the fact that the stability requirement for the con- 
trol program in a DM operating in the closed loop of a 
system has not been considered mandatory by certain 
authors, experiment shows that an unstable control pro- 
gram is inadequate, 

This paper studies the problem of synthesizing the 
control program for specified performance of the 
closed-loop system while taking into account the pro- 
gram realization features indicated above, 


Synthesis of the Control Program 

We shall study the synthesis of the control program 
for a system whose block diagram is shown in Fig. 1, 
The devices for conversion of analog quantities to digi- 
tal quantities and vice versa are referred to the control- 
led object in the closed loop, 
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Assume the object is characterized by the transfer 
function (in the sense of a Laplace transform) 


N (p) 
Wy (p) = M (p)’ 
where N(p) and M(p) are polynomials, M(p)=0 may 
have one zero root*, and the degree of N(p) is lower 
than the degree of M(p). 


Assuming that the DM computes the control signal 
during a time T= const (T is the time required for com-~- 
putation cycle) and assuming that a mandatory condi- 
tion is the presence of a converter at the output of the 
DM fixing the output quantity of the DM over the 
time required for the cycle, we obtain 


k 
W, (P)= do, —z5 


—Tp, -Tp’ 
v=1 e —e 


(1) 
where 


N (T'p,) (1— e7?») 
vy M’(Tp,) Tp, 
: Tp accordance with [3] we use the substitution 


e ‘=z and write the transformation (1) as 
a,2*+ @,_4 a see + ays 

b, 2*4 b_, 2! 4... + bas by” 

where a; and b; are coefficients that depend on C,, and 

4° e7TPy, 





d 
M’ (p) at Sad (p). 


Wis (2) = 





(2) 








4p 


- 








DM Object oe 


























Fig, 1 





* This is the case which is most often encountered in 
practice, 
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Note that (in contrast tothe well-known z transform 
tables) expression (2) takes into account the fixing of 
the DM output (i.e., it takes into account the term 
(1-e" Tpy/ p); therefore the use of (2) simplifies the 
computations which are required for synthesis of the 
control program, 

Note that the zeros and poles of the rational frac- 


tion - Wy (z) are located outside of the unit circle 
on the z plane, provided only that the zeros and poles 
of the transfer function Wy (P) for the object belong to 
the left half-plane of the plane p, This fact, as we shall 
see below, facilitates synthesis of a stable control 
program for specified performance of the closed-loop 
system, 

As we know [3], the control program of the DM 
can be represented by the transfer function 
Ags" + An-12""! +... + Ass + Ay »> & 
By2" + Baz™' +... + Biz + Be 
where A, and B, are the coefficients to be selected in 
the synthesis of the contro] program, 

The synthesis of a stable control program (3) is 
performed for a specified transfer function W,,, (z) of 
the closed-loop system; this transfer function is selec- 
ted on the basis of the requirements imposed on the 
system and the nature of the specified input Xp (its 
rate of change, the presence of random noise), It may 
be written as 


W,, = (2) 





= (4) 
= Wop (2) = Taz + Pez? + Msz? +... + Pnz™ 
giz + go2® + gaz? +... + gy 2! 
Bot 82 + By2* 4---- + itm ik (5) 
Note that in the first case (4) the system has a so- 
called infinite degree of stability and a finite transient 
response time equal to Tm, It is evident that the condi- 
tion for first order astatism in the system will be 


Vp (2) mes 





(6) 
lr +%+0s+...Pm=1. 


For higher order astatism it is necessary that (4) 
contain (1— z) to the appropriate power as a multiplier, 





In the second case, for which W_ (z) is written in 
the form (5), it is necessary that the foots of the denom- 
inator in (5) have a modulus greater than unity in order 
to assure stability of the closed-loop system, 


In the analysis we shall assume that the transfer 
function of the closed-loop system is specified in the 
form (4) or (5), 


By analogy with [4] we subdivide the control pro- 
gram into a series of elements each of which is defined 
by its transfer function (Fig. 2), However, in contrast 
to [4], we choose the transfer functions of the elements 
D,(z) and D3(z). in the form 


W, 
Di (2) = Ws (2), Ds (2) = oe (7) 


Then it follows from a study of the block diagram in 
Fig. 2 that the transfer function for the closed-loop 
system is equal to the specified transfer function; {.e., 
we obtain 


X* = Wop (2) Xo- 

The transfer function D,(z) may be taken in simp- 
lest form, and its coefficients are chosen from the con- 
ditions of proper stability and stability for the closed 
loop of the system (Fig, 2), In other words, the transfer 
function is chosen so that the roots of the equation 

1 +D, (z) W, (z) = 0 
will have a modulus exceeding unity, For example, 
sometimes choosing D,(z) in the form 


Dy (s) = pet he 








buz + dio’ (8) 
where oe > 1, makes it possible to satisfy these 
il 
conditions, 


The resulting control program (Fig. <) is stable; to 
verify this it is sufficient to demonstrate the stability 
of the program corresponding to the element D3(z). 
This is very simple to do; in fact, substituting expression 
(2) into (7), we obtain 


Ds (z) = (Ti + Ta 2 + Moz? +... +3 s™”) (byz*® + byazk—t4, .. bas + be) 





ki 
eye + Op 
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The stability derives from the fact that the zeros 
in the denominator coincide with the zeros of W,,(z) z 
and lie outside the unit circle, It is evident that ‘the 


stable control program for the system in Fig, 2 can be 
written as 


e] = D1 (2) Xq &ey= Ds (2) Xp, (9) 


e, = Dz (ze, —X), e* =e, + ey. 
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The transfer function of the control program ( 3) 
will then consist of two transfer functions: 


e° = Wyy (2) X5 — Dz (2) X*, 


where 





a“ (4) = TyXoq (¢ — T) + TeX (¢ — 27) 4+... 


ey (t) = ee ben (0) — X(N + 5 Les (¢ — 7) - X¢—T) — gee (¢ — 7), 


Way (2) = Di (2) Dz (2) + D, (2). 


Taking (2), (4), (7), and (8) into account, we write 
the control program (9) in the form of the difference 
equations 


+ TmX, (t — mT), 


bu 


19 (10) 


es (t) = CoXe (t) + CiXo (t — T) + CaXo (t— 27) +... + CiXolt— (m + k — 1)T]— 


— Fes (t — T) — Fes (t — 27) —... — Pest — (k— 1) TI, 


@ (t) = &, (¢) + €s (¢). 


Here the values of the quantities in each current 
computation cycle are represented in terms of a func- 
tion of the argument (t), and the values of the quanti- 


ties in the preceding nth cycle are given in terms of 
functions of the argument (t- nT): 


7 ae z Pibo, 
Cr = = (Fabs + Pabo), 


Cy = = (Cabs + Fhr + Tobe), 


It is obvious that it is not difficult to formulate the 
program for computing the control signal in the DM 
using Eq. (10); under these conditions the program vol- 
ume is small, 


Experimental Results 





In order to illustrate the above we shall study an 
example of synthesizing the control program for the 
DM of a system (Fig. 1) for the case where the object 
has the transfer function: 


W, @) = k (Tip + 1) . 
M p(Tsp + 1)(T$p* +2ET7'sp + 1) (11) 





where 
k= 0,0091; T1 = 1.085 sec, T; = 287s = 0,05; 
T3 = 0.0205 sec*. 








Then for a computing-cycle time T= 0,1 sec we 
use (1) and (11) to obtain 





W,, (2) = agz* + a + aez* + ayz e 


(1 — aye T*— 2)(2® + krz the) 


— _ 4424 + asz® + agz* + az 
baz* +- baz® + baz? + biz + by 


(12) 





where a= 0.0625; 4g=—0.045; ay= 0.00571; a, = 0.00771; 
bp = 9.9; by =—24, 75; by = 24.39; by =— 10.54; b= 1; 
Ko = 1.277; ky =— 1.1747. 


We shall choose the transfer function for the closed- 
loop system as 
(13) 
Weph2) s= Ty + Ig2* + Iee® + Maz‘ + 2’, 


where 


mh=T=—lrs=— ly = Fy = 0.2. (14) 
A check of the stability for the closed loop of the 
system (Fig, 2) shows that the transfer function D,(z) 
can be written as the fraction (8), Here the roots of the 
equation 


i +- Dg (2) W, (z) = 0 (15) 


have a modulus greater than unity for a4= 1; ayze"?/ Ty 
= 1.8; by = 0,123 ; by, =— 0,053, 

Substituting the values of the coefficients in (12) 
and (13) into (10), we obtain the control-program equa- 
tion for computing the control signal in the DM: 


es (t) = 0.2 [Xp (t — T) + Xo (t — 27) + Xo (t — 3T) + Xo (t — 47) + Xo (t — ST)I, 

&, (t) = 8.43 fer (t) — X (t)] + 63.5 [e, (¢ —T) — X (t —T)] + 0.42 e& (t —7), 

&s (t) = 31.7 Xp (t) — 47.5 Xo (t — T) 4+ 30.5 Xo (t — 27) — 3.2 Xp (t — 37) — (16) 
—31.7 Xo (t — 57) + 47.5 Xo (¢ — 67) — 30.5 Xy (t — 77) + 3.2 Xo (t — BT) + 
+ 0.72es (t — T) — 0,091 es (¢ — 27) — 0,123es (¢ — 37), 


© (t) = &y (t) + €s (!). 
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Figure 3 shows the graph for the control process in 
the closed-loop system, The control signals were com- 
puted from Eq. (16), In performing the experiment we 
used a specialized high-speed DM which performed the 
computations during a time equal to a portion of the 


cycle T=0,1 sec, We ensured a constant cycle by using . 


a special control network, 

From the graph in Fig. 3 it is evident that in accor 
dance with (15) the system has first order astatism and 
has a velocity gain of approximately 4 sec” *, The seg- 
ment in Fig. 3 corresponding to the initial operation 
of the system after it has been energized is especially 
well illustrated on the typical control process oscillo- 
gram shown in Fig. 4, 

The experimental results not only illustrate the 
described method for synthesizing a control program; 
they also are of special interest in clarifying the speci- 
fics involved in the operation of a system with a DM, 

In fact, the initial segment of the control process 
is usually not taken into account in theoretical papers 
dealing with systems incorporating a DM, This was evi- 
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Fig, 5 


dently the reason for previous underestimation of the 
stability conditions for the control program, 

The indicated segment is not a transient response 
in the general sense of automatic control theory; it a- 
rises due to an absence of information in the “memory* 
of the machine concerning the values of the input and 
the output quantities over the preceding cycles X¢t-hT), 
X(tehT ), oo. (where h=1, 2, so, n),etc,,during the ini- 
tial period of operation for the system, (Incidentally, 
experiment shows that the process corresponding to the 
initial segment may also occur when information is 
lost in the machine due to random failures.) As a result, 
the input of the object will be subjected to a perturba- 
tion (Fig. 3) whose form is determined by the “transfer 
function of the control program” [in the generai case 
(5)] or by the difference equations (12), 

If the control program is stable, then the magni- 
tude of the indicated perturbation decreases; after n 
cycles, where n is the number of preceding cycles which 
are used in the program (in our case eight, Fig. 4), it 
becomes equal to zero, After that the conventional 
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transient response begins in the system in accordance 
with (6), 

For an unstable program the magnitude of the per- 
turbation increases with each subsequent cycle, and this 
leads to an overtaxing of the DM digital capacity and 
cessation of operation, 


The above is verified by the oscillograms for the 
control process (cf, Fig. 5) obtained for an unstable 
control program even though the stability conditions 
for the closed-loop system proper were satisfied. 
Conclusions 

We have analyzed a simple method for synthesi- 
zing a control program for specified performance of a 
closed-loop system when the specifics of the control- 
program realization are taken into account, 

The necessity of satisfying the stability conditions 
for the control program irrespective of the stability 
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conditions for the closed-loop system was verified ex- 
perimentally. 
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In servodrives two channels are often used for 
transmitting signals, The coarse and fine control signals 
are summed at the input of an amplifier after a defin- 
ite matching of the receiver relative to the transmitter 
has been carried out, There are a number of networks 
for performing this summation by means of neon lamps, 
transformers, etc, 

Figure 1 suggests a simpler and more reliable net- 
work for summing such signals, 

The principle governing the operation of the net- 
work is the following: a bias voltage corresponding to 
the operating point on the transfer characteristic is ap- 
plied to one grid of a dual triode, and a negative bias 
voltage with a magnitude equal to or greater than the 
cutoff voltage for the tube is applied to the other grid, 

We shall derive the expression for the summed ac 
component of plate voltage as a function of the mis- 
match angle of the receiver relative to the transmitter, 

The receiver produces modulated fine and coarse 
control signals which can be written as 


Uc = U,, sin 6 sin ot, 


Ur = U,, sin k68 sin wt, 


where U,, Uf are the coarse and fine signals produced 
by the detectors; U_,, is the amplitude of the voltage 














Fig, 1 


tapped off from the detectors; 9 is the mismatch angle 
between the transmitter and the receiver; w is the car- 
rier frequency in sec™*; k is the transfer number be- 
tween the coarse and fine control detectors. 

The grid g, of the tube is subjected to a voltage 
which is phase-shifted through an angle 


Ze 


re + Fin -+ 74 


and has a magnitude equal to 
U.= U,,m sin § sin (ot — ~) — Upjass’ a 





@~ = arc tg 


where 


a! a+ 
V (rin + v5 + may + 2? 
tjz is the internal resistance of the coarse-control de- 
tector; Upias » is the bias voltage at the second grid 
and is determined by the voltage drop across the resis- 
tors fs and ry 
The first grid g, is subjected to a voltage equal to 


Us = Umax "Sin k 6 sin@t — Upjass’ si 





m 





where 


ri . 
~ ri tat ra’ 
Ubjas 1 18 the bias voltage at the first grid and is de- 
termined by the voltage drop across the resistor rg. 
After the signals have been amplified, the plate 
voltage will be 


n 


J, =m (U_,nsinkOsinot— Ui.) + om 

+ ps (U,,m sin 6 sin (ot — 9) — Upjass x 
where py , 22 are the gains for the two arms of the tube, 

Figure 2 shows the waveshapes of the currents in 
the left and right halves of the tubes, as well as the 
waveshape of the over-all plate current for a constant 
mismatch angle; the figure corresponds to the case 
where the voltages U ‘; and U “, are in phase opposition, 
The curve abcd bounds the ordinates for the instanta- 
neous value of the current in the right half of the tube; 
the curve e bounds the ordinates of the current in the 
left half of the tube, and the curve abfd bounds the or- 
dinates of the resultant current, 
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If we approximate the transfer characteristics of 
the tube by a straight line, then we shall have a con- 
stant identical gain for both halves of the tube, In that 
case the ac component of plate voltage (its effective 
value)AU, will be given by the following formula when 
the phase shift angle ¢ is neglected: 


Au, =p (U,, sin 6 r sin 6), (4) 


Um? 
where U‘_= ae is the effective value of the voltage 


at the on 8 arriving from the fine-control detector; 
U "' mn is the effective value of the ac voltage compo- 
nent at the grid g,, corresponding to operation of the 
left half of the tube, 

The graph of the ac component of the plate cur- 
rent as a function of the mismatch angle 9 is shown 
in Fig. 3. 

If we design a network with a deep dip 4u, (this 
can be achieved by increasing the coefficient n), it 
follows that when the drive compensates for a large 
mismatch angle at the controller output, a smaller 
voltage will correspond to this mismatch angle range 
(k® = 360° to 180°); this will cause braking of the con- 
trolled object, 

This fact can be used in certain servodrives (drives 
which have a high overthrust speed, etc.) in order to 
reduce the compensation of large mismatch angles; 
i.e., the rate at which the controlled object arrives at 
the matched position is reduced, 
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Fig. 2 


The bias voltage at the first grid g, of the tube and 
the plate load resistance are chosen from the operating 
conditions of the amplifier based on the fine-control 
voltage; their computation is performed on the basis 
of the requirements imposed on the amplifier, The bias 
voltage at the second grid g, is determined by the re- 
quired range over which operation must be based solely 
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on the fine-control voltage and the tube cutoff voltage. 
The range over which the amplifier (drive) is operated 
by the fine-control signal only (9 ,) is somewhat greater 
than the required operational precision of the drive: 

Ubias2 = Ucot U1, Ui = U,,msin 6; sin (wt — 9), 
where U,._ is the cutoff voltage of the tube and is de- 
termined from the transfer characteristic, 

The resistors rz and r4(Fig. 1) are related by the 

following expressions: 


rs Vbiast 
ra-+rs Vbias2 Ia dite e i dd ra) : Va 





Here I, is the dc component of the plate current, 
and r; is the plate resistance of the tube. 
From these relationships we obtain 


n= (Glon— re) (1 gee) 
a ty Qodgivn 8 5ise Vriase /’ 


= 
rs = 


Ybiasi;, Ubiass 
" Tlasa\ ~ Fines’) 

The capacitor C (Fig. 1) which blocks dc current 
from the coarse-control detector coil must have a max- 
imum capacitance in order to reduce the phase shift 
angle ¢ and increase the coefficient m; however, its 
magnitude is limited by physical size (in practice it 
can be equal to 2 p f). The resistance rs is determined 
from the grid current conditions, In choosing the para- 
meters of the network it is important to know the data 
for the voltage divider consisting of the resistors r, and 
Ip. 


MU 
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Fig. 3 


Using k®9,(* 270°) to denote the angle for which 
the ac component of plate voltage has its maximum dip 
(Fig. 3) due to phase opposition of the coarse and fine 
control voltages, it is possible to write 


m ( U7,sin 0, — U 758in KO, ) = Au, 


m V: 
where Au is the voltage at the input of the last stage 
and corresponds to the beginning of saturation in the 
amplifier or in the drive as a whole, 

From this we determine n: 


V2 
n= (Au+pU,, sin 03) 7, sin kOe 














In designing and checking the network we use the 
following parameters: T is a dual “6N7S" triode;* r, = 
= 27 kohm; r,=1 meg; ts= 1 kohm; r4= 4 kohm; rs= 
= 27 kohm; r,, = 20 kohm; C = 2pnf, 

The signals at the grids of the tube were supplied 
from "SS-405" selsyns, and the transfer number between 
selsyns was 20, 


Notwithstanding its simplicity, the network anal- 
yzed above not only sums signals but also amplifies 
them, 

* The specifications for the "6N'7S" tube are: E¢ = 6.3 v, 
I¢= 0.8 amp, E, = + 300 v, I, = 3.5 ma, E.,=~6 v, 
&m = 1600 ymho, r,, = 22 kh 16 k, Pp = 6 Ww, Pout * 
= 4.2 w [Publisher's note]. 
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